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Abstract 


Over the last years, the increase in energy consumption coupled with ever more 
stringent regulations on pollutants emissions and the massive advent of renewables in 
the energy market, have promoted the development of distributed energy systems and 
thus of an increasing interest towards small and micro power generation systems. In 
this context, the ORC progressively became the leading technology in the field of low 
size energy conversion systems («100 kW) and low temperature applications 
(<150°C). Nonetheless, this technology still deserves further developments, espe- 
cially regarding the design of specific components, which should grant features of 
reliability, acceptable performance level and, often even more important, affordable 
price in order to ensure the attractiveness of the whole energy system. It is the case of 
the small and micro expanders (tens to few kW scale). A possible solution for micro— 
size expanders is the Tesla expander, which is a viscous bladeless turbine that holds 
the desired characteristics of low cost and reliability. This expander was first devel- 
oped by N. Tesla at the beginning of the 20" century, but it did not stir up much 
attention due to the strong drive towards large centralized power plants, where this 
technology becomes no longer competitive against those belonging to bladed expand- 
ers. In the recent years, due to the increasing appeal towards micro power generation 
and energy recovery from wasted flows, this cost effective expander technology rose 
a renovated interest. 

In the present study, a 2D numerical model is realized and a design procedure of 
a Tesla turbine for ORC applications is proposed. A throughout optimization method 
is developed by evaluating the losses of each component and by introducing an inno- 
vative rotor model. The main optimizing parameters of the turbine, such as the rotor 
inlet/outlet diameter ratio, channel width-rotor diameter ratio and tangential velocity— 
rotational speed ratio at rotor inlet are highlighted and assessed. 

The 2D model results are further exploited through the development of 3D com- 
putational investigation, which allows an accurate comprehension of the flow charac- 
teristics, which are difficult to depict with a 2D code. 

Finally, two prototypes are designed, realized and tested. The former one is de- 
signed to work with air as working fluid, with the stator made in ABS with additive 
manufacturing technique, in order to show a possible cost effective way of realization. 
The obtained experimental results of this prototype well match the numerical predic- 
tions. A 94 W net power output with 11.2% efficiency are measured. 

The second prototype is designed to work with organic fluids (specifically with 
R404A), and it is ultimately tested with R1233zd(E). A standard metal manufacturing 


Lorenzo Talluri, Micro turbo expander design for small scale ORC. Tesla turbine, € 2020 Author(s), content CC BY 4.0 International, me- 
tadata CCO 1.0 Universal, published by Firenze University Press (www.fupress.com), ISSN 2612-8020 (online), ISBN 978-88-5518-061-0 
(PDF), DOI 10.36253/978-88-5518-061-0 


Micro turbo expander design for small scale ORC 


is followed for this prototype. The achieved experimental results confirmed the valid- 
ity and the large potential applicative chances of this emerging technology, especially 
in the field of micro sizes, low inlet temperature and low expansion ratios. 371 W net 
power output at 10% shaft efficiency are obtained. 

The experimental results allowed the validation of numerical models, which was 
among the main objectives of this work. In this way, the numerical procedure may be 
reliably employed as the tool for the accurate and optimised design of Tesla turbines 
for organic Rankine cycles but also for applications with gas like air. 

As a final remark, it can be affirmed that the operability of the Tesla expander 
was demonstrated in this work. Thus, it may be considered as a suitable and realizable 
solution to tackle one of the present issues related to micro expanders, namely high 
costs and low reliability, which, moreover, suffers off design conditions only to a 
limited extent. 

The realization of a reliable design tool is another fundamental outcome of the 
present work. 


Chapter 1 
Introduction 


The world scenario recently experienced a strong increase in energy consumption 
demand, associated with a series of issues related to the exhaustion, environmental 
impact and cost of the resources, especially for fossil fuels. This framework encour- 
ages the search of alternative energy solutions for power generation, as well as the 
improvement of already existing conversion systems. 

Over the last years, research on energy systems has focused on small, distributed 
systems for cogeneration, which cover the requirements of heat and power generation 
both in domestic buildings and industrial facilities, with an emphasis on smart grid 
solutions which can effectively deal with problems of load/generation mismatch and 
integration of energy storage. 

When applied to intermediate and low-temperature resources, a modern popular 
technology is the Organic Rankine Cycle, whose applications are being extended to 
small size (5—50 kW.). This technology substitutes water with organic-based com- 
pounds as working fluid. The main advantage of these fluids is that they are suitable 
for low temperature applications, as they allow moderate saturation temperatures and 
pressures and high molecular mass. Indeed, several studies were performed on ORCs 
applied to low-medium temperature thermal resources. Such applications range from 
recovery of heat from gas turbine discharge, internal combustion engines or industrial 
waste heat, energy conversion from biomass, solar or geothermal resources represent 
another common field of application. 

Nonetheless, when micro applications are taken into account, one of the main is- 
sues with Organic Rankine Cycles is linked to the expander, as this component often 
involves high manufacturing costs and offers low reliability. The Tesla turbine, with 
its relatively simple structure, appears to be a potentially reliable and low-cost ex- 
pander, which could find its market in the low-power range. 


1. Motivation 


The application of Tesla turbines to small and mini ORC cycles could allow the 
opening-up of this new niche market, where ORCS have been hindered mainly by the 
high initial investment cost, by delivering an affordable expander technology with 
minimal maintenance requirements. The application to low enthalpy systems will al- 
low the spread of ORC cycles at capillary level, similarly to smart grids, with an EU 
application potential of thousands of units. Therefore, it is of great interest to conduct 
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a research analysis on such an innovative component, which could potentially become 
a breakthrough technology for energy harvesting from industrial wastes of heat and 
low pressure flows, due to its low cost and reliability characteristics. 


2. Objectives and structure 


2.1 Objectives 


This research project aims to the thermo-fluid dynamic assessment of an innova- 
tive boundary layer bladeless expander (Tesla type turbine) for mini and micro energy 
conversion systems, which could become a strong competitor of the actual employed 
micro expanders thanks to its very attractive compromise between efficiency and 
costs. 

The main objectives of the present research can be resumed in the following: 

Development of a numerical 2D model which allows the prediction of the perfor- 
mance of a Tesla turbine for different working fluids, applying real gas assumption 
and introducing sudden expansion and contraction pressure losses; 

Definition of a comprehensive scheme for thermo—fluid dynamic and mechanical 
design and optimization of the expander; 

Development of computational fluid dynamics analysis to depict the flow behav- 
iour inside a Tesla turbine; 

Validation of 2D built numerical model with experimental campaign both on air 
and organic working fluids. 


2.2 Structure 


The manuscript is comprised of five chapters, including the initial introduction 
chapter. 

Chapter 2 is dedicated to the literature review, where the "state of the art" of the 
Tesla turbine researches is assessed. Furthermore, a brief introduction on ORC tech- 
nologies is presented, with a particular focus on micro expanders. At the end of the 
chapter a statistic summary of the available literature on Tesla turbines is reported. 

Chapter 3 is dedicated to the methodology and models utilized in this thesis. Par- 
ticularly the 2D in house EES code is accurately described, presenting each compo- 
nent model. The prototypes design procedure is assessed and the mechanical verifica- 
tion scheme and the computational analysis settings are presented. 

Chapter 4 is dedicated to the analysis of results. The obtained results are divided 
in three main Sections. The first Section dealing with 2D in house code simulation, 
second Section depicting CFD analysis results and last Section displaying the 
achieved experimental data. 

Chapter 5 is dedicated to the conclusions of this research and recommendations 
for future work. 
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Chapter 2 
Literature Review 


1. Word Energy Scenario 


In 2017, the net electricity production grew by 0.8% compared to 2016. A signif- 
icant increase of the power production share (16.7%) was given by renewable ener- 
gies, with a consequent reduction (even if small) of fossil fuels share (by 1%). 

In OECD countries power production by renewable energies accounted for 23.7% 
of the global generation; fossil fuels contribution was of 58.7% and the remaining part 
was filled by nuclear power (17.6%), as displayed in (Fig. 2.1) (IEA, 201 7a). 


Jan-Feb 2018 


10.1% 


ш Combustible Fuels 


ш Nuclear 


8 Hydro 


Geoth./Wind/ 
Solar/Other 


Total = 1856.4 TWh 


Fig. 2.1 OECD Electricity Production by Fuel Type (IEA, 2017) 


The increase in renewable energy share is certainly due by the strong concern 
given by climate change. Particular attention is given to the energy use and green- 
house gases production. Indeed, among human activities that produce greenhouse 
gases, the energy sector is by far the main contributor (68% share (ТЕА, 2017b)). 

The European Union is strongly committed to tackle climate change and it has set 
a comprehensive package of policy measures to reduce greenhouse gas emissions. 
Particularly, H2020 directives on climate change targets the 20—20—20 policy that is 
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of a reduction in 20% of greenhouse emissions (from 1990 levels), a total share of 
energy production by renewable energy of 20% and a 20% improvement in energy 
efficiency. In compliance to these strongholds, the EU policy pushes towards a tran- 
sition to decentralised energy system production, that is through the employment of 
distributed power generation and storage devices in households, as well as to the max- 
imisation of energy recovery from industrial process, which actually waste precious 
resources, such as heat/cold and pneumatic energy (EU commission, 2018). 

This framework encourages the search of alternative energy solutions for power 
generation, as well as the improvement of already existing conversion systems, par- 
ticularly in the field of small and medium power range, which is also the basis to move 
towards the direction of distributed energy systems. Particularly, in recent years, en- 
ergy research focused on small, distributed systems for cogeneration, which cover the 
requirements of heat and power generation both in domestic buildings and industrial 
facilities. Specifically, the affected market ranges from big industrial energy sectors, 
such as textile, food, steel, glass industries to small domestic cogeneration of heat and 
power unites or to inverse cycles (like domestic compression chillers or heat pumps). 

In order to efficiently exploit the waste heat from industrial processes, as well as 
to develop small efficient cogeneration systems, which could also be connected to 
renewable technologies, conventional power generation systems (open cycle gas tur- 
bine and steam cycle) do not seem the most appropriate. Indeed, in the last few dec- 
ades a new technology, based on organic fluids compounds, which are characterized 
by lower saturation temperature and pressure and higher molecular mass when com- 
pared to steam, has taken lead for a wide range of applications where heat and/or 
temperature from the energy sources are limited, such as waste heat recovery appli- 
cations (WHR) or power generation from renewable energies (Fig. 2.2). This technol- 
ogy is known as ORC (Organic Rankine Cycle). 
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Fig. 2.2 Organic Rankine cycle fields of applications (Macchi, 2017) 
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2. Overview of Organic Rankine Cycle 


In order to have of a clear understanding of the reasons of the rising of interest of 
the ORC technology, a comparison with traditional power generation systems both 
from a thermodynamic and a Turbomachinery points of views needs to be carried out. 

First of all, a distinction of power production technologies is given by the archi- 
tecture of the cycle, “open” or “closed”. In open cycles the working fluid experiences 
material exchanges with the environment, both at inlet and outlet of the cycle; an ex- 
ample is the gas turbine cycle, utilizing air as working fluid. Closed cycles, on the 
other hand, are characterised by a working fluid that consecutively operates a cyclic 
series of thermodynamic transformations; an example is the Rankine (or the Hirn) 
cycle, which uses water (steam) as working fluid. Another important aspect to remark 
is the possible transformations that can take place in a power generation system with 
external heat sources (excluding therefore internal combustion engines), which are: 
nearly adiabatic transformation (typically, in pumps, compressors and turbines/ex- 
panders) and nearly isobaric transformations (typically, in heat exchangers). 

Cherishing the above-mentioned difference in cycle architectures and the possi- 
ble thermodynamic transformations, the open-air cycle will be first analysed through 
second law efficiency assessment. Assuming a fixed constant temperature for the heat 
source and a fixed ambient temperature of a simple open cycle, the cycle efficiency 
can be expressed as shown in Eq. (2.1). 


п= (1-22) - (HONS) Q.1) 


1 
max Qin 


To 


e 1— 
max 


any traditional thermodynamic cycle can achieve; 


is known as the “Carnot” efficiency, which is the upper limit that 


AS; . з n" 
• TEN o is the sum of the losses related to each cause of irreversibility. 
in 


Particularly, the second term of Eq. (2.1) can be decomposed in 8 main losses, as 
suggested in (Macchi, 2017) and shown in Fig. 2.3: 


e pressure losses; 

e  fluid-dynamic losses in compressor; 

heat transfer losses in the heat introduction process; 
e  fluid-dynamic losses in expansion process; 

losses due to mixing of hot air to atmosphere; 

heat losses to the environment; 

e  mechanical/electrical losses; 

heat transfer losses in the recuperator (if present). 
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Fig. 2.3 Second Law efficiency and efficiency losses at various heat source temperature. Ney 
is Carnot efficiency, (a) consider simple cycle and (b) recuperative cycle. Optimal cycle pres- 
sure ratio at each temperature is considered (Macchi, 2017) 


As can be noted from Fig. 2.3, the second law efficiency for open cycle architec- 
ture decrease drastically for lower temperatures, due to the increasing of the various 
losses. 

Comparing gas cycles to closed—loop Organic Rankine cycles for temperature 
values below 400°C, the advantages of the ORC solution are quite relevant. First, a 
better coupling of both high and low temperature heat transfer processes can be real- 
ized more easily; in subcritical Rankine cycle, evaporation and condensation pro- 
cesses take place, allowing for large parts of transformations a constant temperature 
heat exchange. This feature is particularly appreciated for heat transfer with the envi- 
ronment, which often requires a relevant heat capacity, and it ensures a major lower- 
ing of the irreversibility in the process of heat transfer. Furthermore, pressurization of 
the cycle can be obtained using pumps (liquid conditions) instead of compressors (gas 
conditions), reducing greatly the amount of work required (and the irreversibility in 
the process). 

Taking as reference the analysis conducted in (Macchi, 2017), where three different 
fluids (water, benzene and MDM) were utilized in order to estimate the efficiency of 
a Rankine cycle with an upper resource temperature of 240°C, it can be claimed that 
Rankine cycles can reach efficiency which are closer to the upper Carnot limit when 
compared to gas cycles. Particularly, as shown in Fig 2.4, the reachable efficiencies 
by a Rankine cycle are in the range of 70-85% of the maximum achievable efficiency 
(compared to the 30% in the case of the gas cycle). Particularly, it is seen, that even 
if the three fluids have very different molecular structures, the achievable cycle effi- 
ciency (when recuperated architecture is utilized) is very close between one and an- 
other. The assumed conditions for the analysis conducted in (Macchi, 2017) are re- 
sumed in Tab. 2.1. 
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Table 2.1 Assumed variables of analysis conducted in (Macchi, 2017) for comparison of gas 
cycles and Rankine cycles 


Variable Assumed value 
Ambient Temperature 15°C 
Condensation Temperature 30°C 
Evaporation Temperature 240°C 
Pump efficiency 0.85 
Turbine efficiency 0.85 
Pressure losses 10% of evaporation pressure 
Thermal losses 1% of heat input 
Mechanical/electrical efficiency 95% 
0.45 
0.40 
0.35 
0.30 
0.25 
(si 
0.20 
0.15 
0.10 
0.05 
0.00 
Ideal Real Ideal Real Ideal Real 
Benzene Water mdm 


Fig. 2.4 Second law efficiency for three different saturated Rankine (ideal and real) cycles with 
assumed condition resumed in Tab. 2.1 The cycle losses represented consider: fluid—dynamic 
losses in pump, fluid—dynamic losses in turbine, pressure losses, heat transfer losses in the 
liquid preheating, heat transfer losses in the evaporation process, heat transfer losses in the heat 
rejection to environment, mechanical/electrical losses, heat losses to the environment, heat 
transfer losses in the recuperator (Macchi, 2017) 


After the comparison between gas cycles and Rankine cycle for low temperature 
heat sources, the reasons why organic fluids are preferable to water for low—tempera- 
ture energy resources are highlighted. The first issue when dealing with the steam 
Rankine cycle for low temperature application is the wet expansion process. Indeed, 
as displayed in Fig. 2.5, the expansion of a saturated cycle is within the liquid—vapour 
dome, on the other hand, for organic compounds, with higher molecular complexity 
(increasing molecular complexity modify the inclination of the vapour curve, known 
also as backward vapour line) the expansion can be dry, which will guarantee that no 
blade erosion issue will present. Furthermore, in the steam Rankine cycle, in order to 
have high turbine performances, the expander design is very costly, as a correct design 
will involve multi-stage turbine, with variable speed shafts. Indeed, for low power 
ranges the construction of an efficient steam expander becomes very difficult, as the 
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steam flow rate would be drastically small with conversely relatively a high expansion 
ratio. Also the development of steam volumetric expanders is subject to many nega- 
tive issues, especially regarding the complexity of the expander (appropriate lubrica- 
tion system, high friction losses, difficulty to realize an adiabatic expander). 


R1233zd(E) 


Т PC] 


Fig. 2.5 Temperature — Entropy diagrams of saturated Rankine cycles for Water and 
R1233zd(E); Evaporator temperature of 150°C, Condensing temperature of 30°C 


Once the thermodynamic (when compared to gas cycle) and turbomachinery 
(when compared to steam Rankine cycle) advantages of ORC for low temperature 
heat resources have been assessed, a comprehensive review on the applications, work- 
ing fluid selection and expanders utilized is required in order to fully grasp the strong 
and weak points of this technology, as well as its maturity level. 


3. ORC State of the art 


3.1 Overview of ORC market 


The first concept of the ORC technology was developed by T. Howard in 1826 
(KCORC, 2018), developing a system to produce 18 kW of energy with ether as work- 
ing fluid. After the first spark, ORC technology started being extensively investigated, 
but at first it was confined to niche markets, as the safety conditions of the power 
plants were not adequate. Therefore, it took a whole century before the first example 
of “modern” ORC was realized by D’Amelio at University of Naples (Tartiére et al., 
2017, Bronicki, 2017). Specifically, the developed ORC utilized solar energy as heat 
source for single stage turbine running with ethyl chloride as working fluid. Finally, 
it is in the 1960s that ORC finally bloomed thanks to the research work of Tabor and 
Bronicki (founder of Ormat technologies) at the National Physic Laboratory in Israel 
and of Angelino, Macchi and Gaia (this last founder of Turboden Ltd.) at Politecnico 
di Milano. 

ORMAT was founded in 1964 and Turboden in 1970. These two companies are 
still today the biggest players in the ORC market. In more recent years, many new 
companies have been established; a list of the major manufacturers, as well as the total 
number of installed units and capacities is shown in Tab. 2.2. 
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Table 2.2 List of manufacturers, data updated at 31* December 2016 (Tartiére et al., 2017) 


Manufacturer кө ee Manufacturer xin М 
АВВ 2 3.8 gt- Energy 2 0.7 
Tech 
Adoratec 23 16.4 Johnson Control 1 1.8 
БЕ ce 20 3.6 Kaishan 40 272 
tional 
Canetix /CETY 50 6.3 Opcon 3 2.0 
DurrCyplan 6 1.2 Orcan 16 0.3 
Electratherm 55 3.14 ORMAT 1102 1701 
Enerbasque 3 0.13 Rank 5 0.07 
Enertime 2 1.6 TAS 17 143 
Enex 1 9.3 TMEIC 1 1 
Enogia 11 0.26 Triogen 37 5:2 
Exergy 34 300 Turboden 267 363 
General Electric 6 101 UTC Powr 10 2.8 
GMK 18 5.3 Zuccato 21 1.7 


ORC units. The many fields of applications were also considered and resumed as 
shown in Fig 2.6. One of the main and unparalleled advantages of ORC technology is 
that it is adaptable to any heat resource, allowing its application in many fields: from 
heat recovery applications at gas turbine discharge (Invernizzi et al., 2018, Chacartegui 
et al., 2009, Al-Sulaiman et al., 2010), or in internal combustion engines (Vaja et al., 2010) 
or industrial waste heat (Bonilla et al., 1997, Desideri et al., 2016), to energy conversion 
from renewables such as biomass (Schuster et al., 2009), solar (Zhai et al., 2009, Lentz et 
al., 2006), geothermal resources (Heberle et al., 2006, Fiaschi et al., 2017), as well as to 
micro-scale CHP units (Dong et al., 2009, Qiu et al., 2011,Tchanche et al., 2011, Tempesti 
et al., 2012). Furthermore, the ORC technology allows also the harnessing of ocean 
thermal power gradient (OTEC) (Bombarda et al., 2013). 
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Fig. 2.6 Market share per application and per manufacturer (Tartiére et al., 2017) 
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As can be noted from Fig 2.6, Ormat is the leader in ORC technologies and the 
principal application to which this technology is associated is energy conversion in 
geothermal power plants (especially in the USA where the total installed capacity of 
geothermal ORC plants is about 750 MW). Nonetheless, another application, which 
is emerging in the last few years, is related to industrial waste heat recovery. In this 
field, ORMAT is still the principal manufacturer, but the supremacy is not defined, as 
both General Electric and Turboden owns a considerable share of installed power. 
Actually, the main application in the heat recovery field is recuperation from Diesel 
engines or gas turbines exhaust gases, as shown in Fig. 2.7. It is to be remarked that 
there is still plenty of room for expansion in this sector, as witnessed by the low share 
in energy intensive industries such as metal, cement and glass sectors. Finally, the 
evolution of the total installed ORC capacity is being displayed in Fig 2.8. Particu- 
larly, it is interesting to note how the global energy prices, represented by the crude 
oil prices directly influences the ORC market. Before XXI century, ORC technology 
was basically only applied to geothermal application, but as the crude oil prices started 
rising and environmental issues emerged, these new technologies started expanding 
in the other sectors, especially in the biomass heat recovery applications. 


ICE or Gas Turbine 65% 


Waste to energy 8.8% 


Metals 7.5% 


Cement & Lime 6.6% 


~ Glass 4.7% 


Biogas 1.5% 
LNG 1.3% 
Petroleum & coal 1.3% 
Landfill ICE 0.9% 
Chemical industry 0.8% 


Others 2% 


Fig. 2.7 Installed capacity share in heat recovery application (Tartiére et al., 2017) 
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Fig. 2.8 Installed ORC capacity from 1984 up to 2017, different application highlighted (Tarti- 
ere et al., 2017) 


3.2 ORC architectures 


ORC technology presents several degree of optimization in the design process. 
The main ones can be resumed to 5: operating parameter selection (mainly pressure), 
working fluid selection, component selection, cycle architecture and control strategy. 
Differently from steam Rankine cycle, ORC can easily incorporate components such 
as re-heaters or recuperators, as well as to develop supercritical cycle at relatively low 
heat sources temperature. These features allow a possible enhancement in cycle effi- 
ciency, as well as the feasibility of cost-effective solutions. 

The simplest structure of an ORC is often called basic ORC and is composed of 
4 main components: two heat exchangers (an evaporator and a condenser), a pump 
and an expander, as shown in Fig. 2.9a. In order to increase the thermal efficiency of 
the simple basic cycle, various modifications have been proposed, from the simplest 
with a recuperator (Fig 2.9b), to multi pressure configuration (Fig 2.9c), flash (Fig 
2.9d) and ejector layout (Fig. 2.9e). 

The ORC with recuperator enhances thermal efficiency though the utilization of 
a recuperative heat exchanger at turbine exit. The exhaust gases at the turbine pre heat 
the fluid at evaporator inlet, allowing thus a reduction of the heat input to the ORC, 
while maintaining the same level of power output (Braimakis et al., 2017). 

The multi pressure level configuration has been introduced in order to decrease 
the irreversibility in the heat transfer process at the evaporator. The match between 
the higher-pressure level and the heat source is improved and this allows an increase 
in thermal efficiency. It is to be remarked that an accurate optimization of the two— 
pressure level is fundamental in order to obtain the lowest level of irreversibility. The 
drawback of this configuration, which enables the global thermal efficiency of the 
power plant, is the increase in the complexity of the layout, and consequently of power 
plant costs (Manente et al., 2017, Shokati et al., 2015). 

The ORC flash cycle is used when zeotropic mixtures are employed as working 
fluid. The advantage of zeotropic mixtures is that evaporation does not occur at fixed 
temperature, allowing hence a better heat transfer process. Nonetheless, as for the 
multi-pressure solution, even if there is an increase in thermal efficiency, there is also 
an increase in power plant complexity and costs (Lee et al., 2016). 


23 


Micro turbo expander design for small scale ORC 


Ejector type ORCs have been investigated with the aim of increasing the expan- 
sion ratio of the turbine, allowing therefore a higher power production. The simplest 
configuration is obtained introducing an ejector at turbine exhaust and a second stage 
evaporator. The second stage evaporator allows the vaporization of the primary fluid 
of the ejector, which consents the reduction of the pressure of the turbine discharge 
stream (Li et al, 2012). Other ejector ORC configurations have been assessed; intro- 
ducing the ejector in parallel to the turbine, using a bleed of the turbine as the hot 
primary stream, or utilizing the hot stream at turbine exit as the primary flow of the 
ejector in order to develop a combined ORC-refrigeration cycle. The disadvantages 
of the ejector ORC configuration are the presently low ejector efficiencies and the 
increase of complexity of the cycle when compared to simple and recuperative ar- 
rangements (Chen et al., 2018, Palacz et al., 2018, Zhang et al., 2016). 

Not only ORC layout can be listed as possible architecture modification, but also 
advanced thermodynamic cycles need to be taken into account. Indeed, several solu- 
tions, which utilize basically the same components configuration, can take place, such 
as trilateral cycles, super critical cycles or trans—critical cycles. Trilateral cycles, or 
often-called triangular cycle, employ the same components of a basic ORC, with the 
exception of the expander, which is specially designed to work in two-phase condi- 
tions. The main advantage of a triangular cycle would be the possibility of optimal 
matching of evaporator temperature profiles. Nonetheless, actually, there are not effi- 
cient expanders which work in two-phase region (Yari et al., 2015). 

Both trans-critical (TC) and super—critical (SC) cycles bypass the two-phase re- 
gion when heated, allowing a good thermal match between working fluid and heat 
source. The difference between TC and SC cycles is that in the former one, conden- 
sation still takes place in the two-phase region. These two cycles configurations allow 
an increase in power production; however, the thermal efficiency is usually compara- 
ble or even a bit lower than other cycle configurations, due to the increase of heat 
demand (Astolfi et al., 2018). Finally, Tab. 2.3 resumes the possible cycle architecture, 
as well as the challenges related to each configuration. It seems clear that, as it often 
occurs in engineering design, there is not an optimal solution, but each configuration 
has advantages and disadvantages which should be carefully weighted and assessed 
depending on the specific application. 


Table 2.3 Possible ORC architecture (Lecompte et al, 2015) 


Cycle Modifications Challenges 
: 1 ficial ifl li 
Recuperative Extra heat exchanger Only pane ие 
limit of flue gasses 
Flash Added Separator, throttling valve, Performance comparable to basic 
mixer, extra heat exchanger ORC; many extra components 


Multiple pumps and heat exchang- 


Multi pressure 
ers 


Many extra components needed 


Availability of high efficiency 


Trilateral Two phase expander 
two-phase expanders 
Trans/Super Super critical fluids High pressure; working fluid sta- 
critical bility 
Ejector Added evaporator, ejector Ejector efficiency 
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Fig. 2.9 ORC layout. а) basic, b) recuperative, с) multi pressure, d) flash, e) ejector 
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3.3 Working fluid selection 


As previously stated, the thermodynamic cycle of an organic Rankine cycle ex- 
ploits the same concept of a steam Rankine cycle, nonetheless, the utilization of or- 
ganic fluids adds a further degree of freedom for the system design. Indeed, lifting the 
limitation of water or air as working fluids, it is possible to select the most appropriate 
fluid, which guarantees the optimal compromise between thermodynamic efficiency 
and technical system configuration, as well as plant costs. An example of the increased 
flexibility given by the fluid selection is given by the possibility of exploiting super- 
critical cycle configuration at low temperature of the heat source. Furthermore, there 
is an increased level of versatility in the selection of pressure and density levels within 
the system, which are almost independent from the upper cycle temperature. 

Organic fluids have a lower boiling point than water and higher molecular com- 
plexity, which are properties that make ORCS suitable for small-medium size power 
plants (1—5000 kW) and for heat recovery applications. Nevertheless, organic fluids 
need to comply with many constraints, which may be very tight depending of the field 
of operation. Indeed, optimal organic fluids should be (Colonna et al., 2015): 


e Non-toxic, inflammable, non-corrosive, cost-effective; 

e Detain low (or better zero) global warming potential (GWP) and ozone de- 
pletion potential (ODP); 

e Thermally stable and compatible to sealing material; 

e Good lubricant; 

e Proper heat transfer properties; 


The possibility to choose the most suitable working fluid depending on the appli- 
cation guarantees (i) increased component efficiencies, (11) usually no vacuum con- 
denser (iii) higher cycle performance compared to gas cycle and steam Rankine cy- 
cles. These favourable features gave rise to numerous research studies on the assess- 
ment of the optimal working fluid selection (Colonna et al., 2015, Bao et al., 2013, Qiu, 
2012, Quoilin et al., 2013). 

Organic working fluids can be categorized not only by their molecular structure, 
but also by the slope of the saturation vapour curve (right side), which is a pivotal 
feature of these working fluids. There are three types of working fluid: a wet working 
fluid is called one, which possess a negative slope vapour saturation curve È < 0); 
a dry working fluid conversely, is one which is characterised by a positive slope va- 
pour saturation curve È > 0); finally, an isentropic working fluid is the one which 
has an almost infinite slope. Fig. 2.10 displays some of the most common organic 
working fluids saturation curves. Particularly, it is possible to distinguish wet fluid 
(such as Water, Ammonia or R134a), from dry fluids (such as n—pentane, n-hexane 
or MM) and isentropic fluids (such as R245fa or R1233zd(E)). From the analysis of 
Fig. 2.10 it is possible to easily understand the advantage of dry and isentropic fluids. 
Particularly, wet fluids require super heating in order to avoid wet expansion, con- 
versely to isentropic and dry fluids, which can therefore be optimized through the 
utilization of a saturated cycle configuration. 


26 


Lorenzo Talluri 


400 
| Water 
300 
o n-hexane MM 
= n n-pentane 
R1233zd(E) 
100 
0 
-2 0 2 6 8 10 


4 
s [kJ/kgK] 


Fig. 2.10 Saturation curves of some common organic working fluid compared to water 


The slope of the saturation curve is not the only thermodynamic parameter to be 
taken into account when selecting the optimal working fluid for an ORC. Tab. 2.4 
resumes the principal thermodynamic and physical properties, which need to be con- 
sidered, while designing an organic Rankine cycle power plant for a specific applica- 
tion. 

Another important aspect that needs to be taken into account when selecting the 
working fluid for a specific application is that it has to be easily obtainable and have 
reasonable cost. Good availability and low cost fluids are those already used in refrig- 
eration and chemical industries. Keeping in mind that the assessment of the proper 
working fluid for a specific application needs to be considered in the design process 
of an ORC, an interesting classification of optimal working fluids to be used as func- 
tion of the heat resource temperature level has been carried out in 0, and reported in 
Fig. 2.11. 
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Fig. 2.11 Optimal working fluids as function of heat source temperature (Wang et al., 2013) 
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Table 2.4 Principal thermodynamic and physical properties of an organic fluid (Bao et al., 
2013, Qiu, 2012, Quoilin et al., 2013) 


Property 


Effect 


Vaporization latent heat 


Density 


Specific heat 


Critical temperature 


Boiling Temperature 


Molecular Mass 


Molecular complexity 


Viscosity 


Thermal 
conductivity 


Evaporating pressure 
Condensing pressure 


High temperature 
stability 


High safety 
level 
Environmentally 
Green 


At fixed conditions, at a higher latent heat corresponds a 
higher specific work output; conversely, a lower latent heat 
allows a better thermal match between working fluid and 
heat source, as most of the vaporization process occurs at 
variable temperature. 

High vapour density is fundamental if compact, cost effec- 
tive expander has to be designed. Low vapour density means 
higher volume flow rate and therefore higher expander di- 
mensions. 

Low liquid specific heat allows lower specific work required 
by the pump. 

High critical temperature allows higher thermal efficiency. 
The negative side is low pressure at condenser. Furthermore, 
reduced density, related to working at densities much lower 
than the critical one, leads to higher expander dimensions. 
High boiling temperature leads to higher thermal efficiency, 
if fluids of the same “family” are compared. 

High molecular mass allows higher power production at re- 
duced rotational speed, which leads also to a positive effect 
on turbine efficiency; nonetheless, high molecular mass ne- 
cessitate of high heat transfer area. 

Molecular complexity is directly linked to the slope of the 
vapour saturation curve. Low molecular complexity is linked 
to wet fluids, while higher molecular complexity is charac- 
teristic of dry fluids; high molecular complexity allows 
higher efficiency at reduced turbine size, moreover it fa- 
vours the efficiency of the regenerative cycle configuration. 
Low viscosity in both liquid and vapour phases allows high 
heat transfer coefficients and low friction losses in pipes 


High thermal conductivity allows better heat transfer 


Higher evaporating pressure leads to higher costs and com- 
plexity of the layout of the power plant. 
Condensing pressure should be higher than atmospheric in 
order to avoid air infiltration into the system. 
Maximum heat source temperature limited by the chemical 
stability of the working fluid, which tends to deteriorate at 
high temperatures 


Low toxicity and flammability required 


Low ODP and GWP 
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3.4 Expander assessment 


The utilization of an ORC is a sound solution when the system is composed by 
efficient, reliable and low-cost components. A critical component for ORCS is the 
expander, as it often does not combine low cost and reliability. As the ORC power 
and hot temperature resources vary strongly depending on the application, the variety 
of expanders, which can be considered, is also wide. Turbines (axial or radial) as well 
as volumetric expanders (scroll, screw, vane or piston) are employed depending on 
hot source and power range. 

Axial turbines are often used for plants with power production between 500 kW 
and few MWs (Talluri et al., 2017), while radial turbines are better suited for the lower 
power ranges (50—500 kW), due to their low degree of reaction and therefore their 
capability of dealing with large enthalpy drops at low peripheral speeds, allowing the 
adoption of a single stage design (Fiaschi et al., 2012, Fiaschi et al., 2014, Fiaschi et al., 
2016). Finally, for very small and micro power range applications (500 W to about 50 
kW), volumetric expanders, like scrolls, screws, vane or piston, are usually utilized, 
although their efficiency is limited by leakages, friction and heat transfer losses (Garg 
et al., 2016, Lemort et al., 2009, Ziviani et al., 2014). Several studies have been performed 
on the assessment of expander selection, depending both on power and volumetric 
expansion ratio (Branchini et al., 2013) or application (Quoilin et al., 2012). As can be 
noted both from Figs. 2.12 and 2.13, at higher power ranges the predominant technol- 
ogy is the axial turbine, while for the small-micro power range volumetric expanders 
dominate. For intermediate power ranges (20—200 kW) the specific application and 
boundary conditions are the determining factor for expander selection. 
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Fig. 2.12 Expander technologies comparison as function of power and volumetric expansion 
ratio (VER) (Branchini et al., 2013) 
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Fig. 2.13 Expander technologies comparison as function of power and application (Quoilin et 
al., 2012) 


ORC axial turbines have essentially the same design as steam turbines. Nonethe- 
less, due to the strong difference in thermo physical properties, some peculiar features 
need to be taken into account: 


1. High molar mass implies low sound velocity, which means that when design- 
ing the stator, a particular attention has to be taken in order to limit the exit 
nozzle Mach number to reduce shock losses; 

2. Organic fluids realize smaller enthalpy drops then steam, this affects in a dif- 
ferent way all kind of turbine losses; 

3. High density and small specific volume allow a compact design of the ex- 
pander; 

4. When utilizing organic fluids that are flammable or explosive, special care 
needs to be taken so that no leakages to air are present. 


Considering the above-mentioned features, a one-stage axial turbine is often uti- 
lized for high mass flow rate systems, while for lower mass flow rates radial turbines 
are preferred. The main advantages of radial turbines for low mass flow rate applica- 
tions can be resumed as follows: 


1. High efficiency levels even at off-design conditions (obtained by means of 
variable-geometry inlet guide vanes); 

2. The downscaling of radial turbines is favoured when compared to axial tur- 
bines due to the lower sensitivity to geometric inaccuracies of the blade pro- 
file; 

3. Radial turbines are simpler to manufacture than axial turbines (Bao et al., 
2013); 

4. Higher peripheral speeds than axial turbine are possible, which guarantee a 
higher enthalpy drop per stage; 


Axial and radial turbines are actually not appropriate for micro and small power 
generation units, mainly due to the very high rotational speed, which dramatically 
increase when the turbine power output is lowered. Therefore, for low power genera- 
tion positive displacement expanders are predominant. These technologies will be ex- 
amined in Section 2.3.5, as they can be considered as “direct competitors" of the Tesla 
turbine. 
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3.5 Micro expanders 


As previously stated for small—micro power generation, volumetric expanders are 
actually the only alternative. The principle of operation of positive displacement ex- 
pander is the decrease in pressure through an increase in volume. Commonly, the ex- 
pander is composed by a stator and one or more rotors directly linked to the expander 
shaft. While the shaft rotates, the geometry of the rotor defines a series of closed 
working chambers, where the fluid is enclosed. These chambers increase in volume, 
reducing thus the pressure of the fluid, which transfers energy through shaft momen- 
tum exchange (resulting from the pressure developed on moving surfaces). The major 
difference between volumetric expanders lies fundamentally in the mechanism that 
determines the variation of the volume of the working chambers, as well as the timing 
of the thermodynamic processes (Lemort et al., 2017). 

The principal positive displacement expanders are: reciprocating piston, screw, 
scroll and vane. These expanders hold different characteristics, which distinguish one 
from another. The first feature is the type of motion of the rotor. As shown in Fig. 
2.14, three main motion categories exist: reciprocating, orbital and rotary (Lemort et 
al., 2017). Each expander belongs to one of these classes. Apart from the Roots ex- 
pander, all other volumetric expanders possess a built-in volume ratio where the 
working fluid expands. Another difference between the expanders is represented by 
the presence (reciprocating piston) or not (scroll, screw, vane, roots) of valves. The 
absence of valves enhances the reliability of the expander. 


Volumetric 


Reciprocating Orbital Rotary 


Axial piston Crank train Scroll — Trochoid Vane Roots Screw 


EM ©@@ BF 


Fig. 2.14 Volumetric expander categories as function of motion type (Lemort et al., 2017) 


Scroll expanders are composed by two off-axis spiral-shaped profiles one of 
them is fixed while the other one is turning, delineating a spiral trajectory. The fluid 
enters from an opening, which is placed at the centre of the fixed profile, and after 
expanding through the expansion chambers developed by the rolling of the profile, it 
exits radially at the discharge chamber. As scroll compressors are very common in 
refrigeration and air conditioning applications, this expander is often obtained from a 
compressor just by reversing the direction of rotation and the inlet and outlet of the 
fluid. Fig. 2.15 displays the fluid expansion process inside the assessed technology. 
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Fig. 2.15 Working fluid evolutions inside a scroll expander (Lemort et al., 2017) 


As typical for volumetric expanders, scroll technologies operate on a fixed volu- 
metric ratio, typically between 1.5 and 5 0. Moreover, the optimal sizing for an ex- 
pander is indeed different from that of a compressor, so that using an existing com- 
pressor just reversing the direction of shaft rotation results in general in a poor effi- 
ciency. Two kinds of losses arise if the operating system volume ratio does not match 
the expander design volume ratio; which are under-expansion and over—expansion 0. 
These losses can hinder the efficiency of the expander, limiting the expansion ratio. 
Other assessed losses involve friction losses, supply pressure drop, internal leakages 
and heat transfer 0 0. Scroll expanders can be further categorized in two: compliant 
and cinematically rigid. Compliant scroll tolerates a further movement in a provided 
direction, which let them be adapt for transient operation or for two—phase flow con- 
ditions. A very attractive feature of scroll expanders is that they are very reliable as 
the total number of moving components is quite low. Furthermore, according to 0, 
scroll expanders can reach quite high efficiency (>70%), at optimal pressure ratio and 
rotational speed. 

Screw expanders can be categorized in two main types, single screw and twin— 
screw expanders, as shown in Fig. 2.16. 


Fig. 2.16 Single and twin screw expander (Zhang et al., 2014, Zywica et al., 2016) 


Twin-screw expanders are more common in organic Rankine cycle applications, 
but the research in single screw is increasing in recent years, due to the favourable 
feature of having a better balance than twin-screw, which allows achieving an aug- 
mented bearing lifespan. The fluid enters from one side of the screw and exit on the 
opposite side, passing through a sequence of supply, expansion and discharge cham- 
bers, which are determined by the turning of the screw(s). The rotational speed that 
can be achieved with screw expanders is quite high, which implies the necessity of 
the utilization of gearboxes and control speed equipment. Screw expanders can be 
lubricated or oil-free. Oil-free expanders have the advantage of being mechanically 
simple machines and therefore allow simple manufacturing process; conversely, they 
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hold a main disadvantage, that is, higher internal leakage losses. Nonetheless, screw 
expanders present high efficiency values (up to 80%) with a power range between 1 
and 200 kW (Zhang et al., 2014). A further favourable feature of these expanders is 
their good ability to tolerate two—phase flows. As scroll expanders, also screw ex- 
panders can be easily obtained from a reversed compressor, just changing the direc- 
tion of rotation. The maximum built in volume ratio of screw expander is about 8. 

The Roots expander is similar in principle to a screw expander with a 2D geom- 
etry, but it is not as common. Research studies on this kind of expander are quite few. 
The typical built-in volume ratio is close to one and the power range is typically be- 
tween 1 and 30 kW. The advantage of this expander is that it can handle easily two— 
phase flow conditions. 

Reciprocating piston expanders for organic Rankine cycles are derived from the 
well-known design concept developed in various areas, such as automotive, energy 
or petrochemical industries. The most used design of volumetric expander relies on 
the crank mechanism, but other mechanism can be utilized, such as free pistons 
(Preetham et al., 2016) or slash-plate (Haiqing et al., 2006). The main advantage of piston 
expander is that they are suitable for applications where high temperature and high— 
pressure ratio are required. The built in volume ratio of a piston expander can be as 
high as 15. The efficiency value is around 70% and they are suitable for low displace- 
ment and low power applications. 

The piston expander is the only one between volumetric machines for micro gen- 
eration in organic Rankine Cycle that requires the utilization of valves. This feature 
increases the complexity of the machine, lowering its reliability. Furthermore, there 
are also some other drawbacks, such as lubrication issues and high manufacturing 
costs. 

The Rotary vane expander is characterized by a rotary motion of a vaned cylin- 
drical rotor within a housing, which acts as stator. The geometry of the expander is 
appreciably simple, which guarantees low manufacturing costs. The main advantages 
of the vane expander are the flat efficiency curve throughout an extended range of 
operating conditions, the low operating speed, the possibility of accepting two—phase 
flows and the high expansion ratios achievable (up to 10). Conversely, as most volu- 
metric machine, the vane expander is limited in temperature. Indeed, high inlet ex- 
pander temperatures are not possible, as they would cause an increase in friction losses 
in the vanes due to thermal expansion of the rotor, which finally could end up to be 
locked. 

When going through the selection of an expander for micro organic Rankine cy- 
cles, many parameters need to be assessed, such as efficiency, pressure ratio, ability 
to tolerate two phase conditions, rotational speed, reliability and cost. Which is the 
fundamental feature for expander selection would be determined by the specific ap- 
plication, but some guidelines can be outlined. Particularly, the comparison between 
various types of micro expanders for ORC applications is resumed in Tab. 2.5. As 
previously discussed, in the very small power range, radial turbines are not suitable, 
and actually volumetric type machines are the only alternative. Among volumetric 
machines, scroll and rotary vane expanders аге more suitable for very small-scale 
applications, whereas screw and reciprocating piston expanders belong to a higher 
power output range. Therefore, in this context, the Tesla turbine may represent a direct 
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competitor to scroll and rotary vane expanders, as, if properly designed, it holds the 
same characteristics of moderate rotational speed (if relatively high rotor diameter is 
utilized), low manufacturing cost and suitability to very different fluids and applica- 
tions. Furthermore, conversely to most of volumetric expanders, it does not require 
lubrication, which may be very important in several applications. Section 2.4 will be 
centred on the Tesla turbine, from the principle of operation to the actual state of the 
art. 


Table 2.5 Comparison of micro expanders for ORC applications (Bao et al., 2013, Lemort et 
al., 2013, Zywica et al., 2016) 


Power Rotational 
Type range need [rpm] Cost Characteristics 
Scroll 1-10 <10,000 Ly + High efficiency, low cost 
expander — Lubrication requirement 
+Flat efficiency curve at 
Screw Me- off-design 
— < 
expander Brem Poon’ dium — Difficult to manufacture, 
lubrication 
А { +Mature technology, high pressure 
po E. 20100 «200 е ratio 
P р — Heavy weight, complex 
È + i 
Rotary vane ex 1—5 <10,000 Low Low cost and low noise 


pander —Small power range, lubrication 
+Light weight, mature technology 
50-500 5,000-80,000 High —High cost, low efficiency in 
off-design 
+Low cost, low noise, moderate 
efficiency, reliable 
— Few prototype tested 
(very low TRL) 


Radial inflow 
turbine 


Tesla turbine 0.5-10 <10,000 Low 


4. The Tesla turbine 


4.1 Principle of operation 


A competitive technology to the actual available micro—expanders is the Tesla 
turbine, which is a viscous bladeless turbine. This concept was first developed by 
Nikola Tesla at the beginning of the 20" century, but it went through a long period of 
indifference due to the run towards large size centralized power plants. Only recently 
it found a renewed appeal, as its features make it suitable for utilization in small and 
micro size systems, like ORC applications, where low cost components become very 
attractive for the exploitation of residual pressure drop. The first description of the 
turbine (also called friction or disk turbine) was given in the patent submitted by Tesla 
(Tesla, 1913) (Fig. 2.17). This type of radial expander is characterized by the absence 


34 


Lorenzo Talluri 


of rotor blades, which are replaced by multiple parallel flat disks; a little gap separates 
the rotor disks from the related stator parts, which consist of one or more tangential 
nozzles. The working fluid accelerates, expands through the nozzles and enters, al- 
most tangentially, in the gaps between the disks, where it depicts a spiral centripetal 
path. The working fluid moves from the inlet to the outlet radius due to the difference 
in pressure determined by friction and by the exchange of momentum, and exits from 
openings made on the disks at the inner radius. Tesla turbines seem to have several 
advantages when compared to conventional expanders for low power generation, as 
their relatively simple structure allows a straight manufacturing process, as well as 
low cost, reliability, modularity, and versatility. On the other hand, Tesla turbine pre- 
sents values of efficiency lower than those of conventional turbines. The most critical 
aspects are the design of the nozzle and of the jet velocity profile. 


Fig. 2.17 Figs. from Tesla patent, 1913 (Tesla, 1913) 
4.2 Literature review 


In this Section the available literature research studies will be briefly assessed. A 
chronological order will be followed. 

1950 

After Tesla pioneer work, it was only in the 1950s that the Tesla turbine was fur- 
ther investigated. Particularly, Leaman AB, (1950) designed and built a prototype of a 
Tesla turbine utilising air as working fluid. The rotor diameter was of about 0.13 m; 
the maximum assessed power output was of 87 W and the maximum documented 
efficiency was 8.6%. 

1952 


Armstrong JH, (1952) carried out an accurate study on the design and realization 
of an experimental test rig to investigate the power and efficiency of the disk turbine. 
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He conducted a test campaign focusing on the identification of the critical issues of 
the machine, using steam as working fluid; furthermore, different nozzle configura- 
tions were analysed. A valuable result was the understanding of one of the causes of 
inefficiency. Indeed, it was found that the nozzle flow strongly affects the perfor- 
mance of the turbine. 


1961 


Beans EW, (1961) performed both a numerical and experimental investigation on 
the performance of the friction disk turbine. The performance model developed was 
reliably predicting the performance of the turbine. The performed campaign of inves- 
tigation allowed assessing a maximum turbine efficiency of over 24% and a power 
near 1.8 kW. 


1965 


Rice W, (1965) realized one of the first analytical/numerical models of the flow 
inside the Tesla turbine based on simplified Navier-Stokes equations, assuming a 
steady, incompressible and inviscid flow; moreover, he designed and tested six dif- 
ferent disk turbines operating with air. Rice W discovered that the best efficiencies 
were achieved with small size turbines operating at low flow rates, in contrast with 
conventional bladed turbines, and suggested the application to small power range, 
exploiting the qualities of low cost, ease of manufacture, low noise level and reliabil- 
ity. 

1966 


Beans EW, (1966) performed a throughout investigation of a 0.15 m scale air Tesla 
turbine, following the work previously developed in 1961. The performance investi- 
gation included the assessment of various channel width (from 0.6 to 12.7 mm), var- 
ious rotational speeds (from 4000 to 18000 rpm) and various turbine inlet pressure 
(from 1.7 to 3.75 bar). The assessed efficiencies ranged between 7 and 25%. The 
maximum obtained power was of 1.8 kW at 3.75 bar of supply pressure and 1.5 mm 
of channel width. 


1967 


Matsch L and Rice W, (1967a) developed a mathematical model of the potential 
flow of an incompressible fluid between two circular disks. The main feature of the 
analysis was the weighting of the partial admission effects on the flow developing 
inside the disks. A follow up of the previous work included the effect of rotation of 
the disks with a flow at low Reynolds number (Matsch L and Rice W, 1967b). 


1968 


Boyd KE and Rice W, (1968) developed a numerical model on the laminar incom- 
pressible flow of a Tesla turbine rotor with full peripheral admission. Particularly, 
various conditions were investigated, changing the three main input parameters: flow 
rate, Reynolds number and mean tangential velocity component at inlet. 

Matsch L and Rice W, (1968) formulated an asymptotic fully developed flow so- 
lution of the flow of an incompressible fluid inside two rotating disks. Particularly, 
the asymptotic solution is function of two main parameters, Reynolds number and 
mass flow rate. 
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1970 


Adams R and Rice W, (1970) performed an experimental investigation of the flow 
characteristics inside a Tesla turbine rotor. The focus on the analysis was the compar- 
ison between experimental measurements of pressure throughout the disks and the 
numerical model developed in (Boyd KE and Rice W, 1968). Experimental results were 
faithfully predicted by the numerical model. 


1971 


Boyak BE and Rice W, (1971) developed an integral solution for the three—dimen- 
sional laminar flow inside a Tesla turbine. The developed model was validated against 
experimental data and the matching with static pressure distribution was found to be 
very good. The model did not take into account compressibility nor heat transfer, but 
it was stated that these improvements could be easily achieved. 


1972 


Lawn MJ, (1972) exploited the numerical models developed in (Boyak BE and Rice 
W, 1971) in order to draw the performance maps of a Tesla disk turbine. The perfor- 
mance maps were realized for various dimensionless tangential velocities at several 
Reynolds numbers. Efficiency as high as 8096 was predicted. 


1974 


Lawn MJ and Rice W, (1974) presented the performance maps of Tesla turbines 
highlighting how efficiency deeply depends on geometry and rotational speed. Partic- 
ularly, nozzle direction was assessed as a fundamental parameter. It was also claimed 
that the Tesla turbine could be designed with a high efficiency configuration for any 
fluid. 

Pater et al., (1974) performed an experimental campaign in order to investigate 
the behaviour of the flow inside two corotating disks. Static pressure alongside the 
disks was measured, as well as flow visualization (through the utilization of a dye) 
was performed. Previous analytical models developed both for pumps and turbines 
were thus verified, showing excellent agreement between analytical and experimental 
results. 


1975 


Bassett CE, (1975) developed an integral solution of the compressible flow of 
Tesla turbines. The developed model was validated against published experimental 
data, demonstrating a good agreement between numerical and experimental results. 
The model took into account several parameters, such as inlet Mach number, Reyn- 
olds number, rotor diameter, mass flow rates. The maximum efficiency was obtained 
for the minimum assessed value of flow rate between the disks. It was also found that 
the optimal rotor radius ratio is between 0.25 and 0.3. 


1976 
Garrison et al., (1976) presented a numerical model for the prediction of laminar 
compressible flow between Tesla turbines. Particularly, the model efficiency predic- 


tion reached values as high as nearly 9096, with low mass flow rates and an optimal 
rotor inlet/outlet diameter ratio between 0.2 and 0.4. 
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Steidel R and Weiss H, (1976) performed an experimental investigation of a Tesla 
turbine for geothermal applications. The tests were performed in wet conditions, with 
a vapour fraction between 6 and 15%. The maximum turbine efficiency observed was 
of a 6.8% at 4000 rpm; the maximum assessed power was of about 2.8 kW at 4000 
rpm. 

1978 


Truman CR etal., (1978) developed a numerical model of a laminar flow between 
two rotating disks of a vapour containing liquid droplets. This type of flow condition 
is typical of geothermal applications. Vapour—liquid interaction was modelled on the 
assumption of a drag force coefficient. 


1979 


Truman CR et al., (1979) modified the model developed in (Truman CR et al., 1978) 
in order to simulate flows containing solid particles. The results showed that for very 
small particle sizes the behaviour of the fluid would not change; but when the particle 
size increases (>4 um) the solid particles would not follow the trajectory of the fluid 
and could also not enter the rotor. 


1990 


Allen JS, (1990) exposed a closed—form model for velocity and pressure inside a 
Tesla turbine rotor, achieved through the assumption of fully developed boundary 
layer. A closed formulation for evaluation of torque was also formulated. The detailed 
solution of the Navier Stokes equation reduction was presented. 


1991 


Rice W, (1991) presented a state of the art analysis of the Tesla turbine. Particu- 
larly, the available models in literature were collected and discussed; an assessment 
on the possible applications was also carried out. High efficiency, as high as 95%, 
were claimed for very small mass flow rate, as well as several field of application, 
from small shaft power application to use of very viscous fluids or geothermal appli- 
cations. 


2001 


Sandilya P et al., (2001) developed a numerical assessment of the fluid dynamic 
and mass transfer in a Tesla turbine rotor. The effect of rotation on the shape of ve- 
locity was highlighted and compared to the stationary case. Furthermore, the flow 
pattern, as well as its effect on mass transfer coefficient, was carefully assessed. 


2002 


Patel N & Schmidt DD., (2002) carried out an experimental campaign on a bound- 
ary layer turbine using biomass combustion gases as working fluid. The 40 hours 
campaign was executed in order to verify the effect of deposition, erosion and corro- 
sion of the Tesla turbine due to the substances present in the flue gases. The maximum 
obtained efficiency was of 1196 with 3.2 kW of power output and a rotational speed 
of 6284 rpm. The same turbine was tested also with steam as working fluid, allowing 
a power production of 9.3 kW with an efficiency of 13.7% at 6500 rpm. 
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2003 


Lezsovits F, (2003) investigated the possibility of utilizing a Tesla turbine in a 
decentralized power generation system based on biomass. The selection of the turbine 
was justified by the requirements of reliability, robustness and cost effectiveness; and 
on the claim that turbine efficiency was not a fundamental parameter for the specific 
application. 


2004 


Ladino AFR, (2004a, 2004b) developed a computational fluid dynamic assessment 
of a Tesla turbine of a 100 mm rotor outer radius, with air as working fluid. Efficiency 
of around 20% was obtained for low mass flow rate when both laminar and turbulent 
model were applied. The visualization of the path lines, as well as the highlighting of 
pressure losses in the nozzle-rotor interaction was reported. 


2006 


Couto et al., (2006) presented a simple and straightforward model, which allowed 
assessing the optimal number of channels of a Tesla turbine in order to be applied to 
a specific application. 


2007 


Batista M, (2007) provided an analytical solution of a steady state, incompressible 
flow of fluid evolving between two rotating disks. The model was obtained through 
the application of asymptotical evolution of unknown components of velocity and 
pressure. 

Bloudicek P and Palousek D, (2007) presented a simplified design procedure for 
the design of a Tesla turbine. The process design, as well as the realization of an ex- 
perimental campaign with water as working fluid was assessed. The effectiveness (ac- 
tual machine performance over theoretic machine performance) was evaluated and 
the maximum reached value was 55% with a 13.63 W power output. The maximum 
obtained power output was of about 58 W. A significant data, giving evidence to the 
most attractive feature of the turbine, was the total cost of the prototype, which was 
of about 124 €. 


2008 


Deam et al., (2008) evaluated the benefits of utilizing Tesla turbines instead of gas 
turbines for small-scale applications. The development of scaling laws for gas turbine 
was refined and the results from the model were compared with experimental data on 
Tesla turbine. The main result was the demonstration that for small—sized turbines, 
viscous turbines are more efficient than conventional bladed turbines, as the losses 
are quite high. The assessed experimental turbine efficiency was 23.5%. 

Lemma et al., (2008) performed a comprehensive experimental and numerical 
study on a 50 mm rotor Tesla turbine. The assessed performance of the turbine was 
over 20%, claiming that the main causes of losses were mainly parasitic losses and 
specifically bearing losses. If these losses could be avoided, or at least reduced, the 
turbine efficiency could reach values close to 40%. 
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Valente A., (2008) applied the Tesla turbine concept as an equipment for pressure 
reduction of Hydrocarbon gases in a near isothermal process. Experimental tests con- 
firmed the wanted characteristic of near isothermal expansion. This allows the ex- 
pander to be used as a substitute of a Joule-Thomson valve in hydrocarbon pipelines. 


2009 


Crowell R, (2009) is the first author to take Tesla turbine as a possible expander 
candidate for small ORC applications. Particularly, the study assessed an integrated 
system for residential applications, exploiting solar hot water collectors as heat source 
with the utilization of a micro ORC in order to produce electricity. The fluid utilized 
in the simulation was Care30 (a refrigerant blend of Isobutane and Propane with sim- 
ilar characteristics to R134a). The developed model of the turbine was quite simple, 
as the aim was to assess the power output at the shaft. A turbine efficiency of 30% 
was assumed, which guaranteed a daily power production of 3.2 kWh. 

Hoya GP and Guha A, (2009) designed and manufactured a flexible test rig for 
Tesla turbines. They carried out several experimental analyses, comparing various 
measurement methods and developed a new, simple and cheap approach (angular ac- 
celeration method) for measuring torque and power output, which overcame the dif- 
ficulties associated with the determination of very low torque at very high rotational 
speed. 

Guha A and Smiley B, (2009) investigated the nozzle, recognizing it as the source 
of the major irreversibility, according to their test results; they demonstrated that a 
careful design of the nozzle could reduce the nozzle losses by 40-50%. They showed 
that utilizing a plenum chamber could result in a total pressure loss of less than 1%; 
furthermore, it allowed also a considerable enhancement in the uniformity of the jet 
and an improved overall efficiency of the Tesla turbine. 

Lampart P et al., (2009) developed a throughout CFD investigation on different 
Tesla turbine dimensions with Solkatherm SES36 as working fluid. The assessed ef- 
ficiency of the turbines ranged from 30 to 50% depending on turbine size and rota- 
tional velocity. For a 0.32 m shaft diameter the assessed efficiency was around 50% 
and the power output obtained was between 1.5 and 5 kW. 


2010 


Carey VP, (2010a) realized a one-dimensional idealized model of momentum 
transfer in the rotor, and used it to predict the turbine efficiency in a 4 kWe solar 
Rankine cycle combined heat and power system, using steam as working fluid. His 
model resulted similar to the one developed by Rice W, (1965), using body-forces to 
represent the wall shear effects, assuming laminar flow: according to it, under optimal 
design conditions, 75% turbine isentropic efficiency was predicted. 

Carey VP, (2010b) exploited the two-dimensional model developed in (Carey, 
2010a) and compared the possible achievable results that could be obtained with CFD 
analysis; concluding that 1D/2D model is a useful tool for preliminary design of Tesla 
turbines. 

Emran et al., (2010) developed a simple model for estimating the proper sizing of 
a Tesla turbine dynamometer and generator. Tesla turbine static torque model was 
assessed, both analytically and experimentally. 
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Puzyrewski R & Tesch K, (2010) developed a 1D analytical model of the Tesla 
rotor calibrating it through 3D CFD analysis. The calibration through 3D CFD al- 
lowed to realize very simple and reliable 1D model, which provided the prediction of 
high efficiency (above 75%) when the gap between the disks is very tight (0.1 mm). 

Romanin V et al., (2010) investigated the possible power and efficiency enhance- 
ment that could be obtained by a commercial Tesla turbine. Particularly, the compar- 
ison between experimental data and numerical simulation where performed. The max- 
imum assessed experimental efficiency was 10.1%, while the numerical optimization 
of the turbine predicted that efficiency higher than 75% could be reached. 


2011 


Batista M., (2011) re-organised the work presented in (Batista, 2007), providing an 
analytical solution of a steady state, incompressible flow of fluid evolving between 
two rotating disks. 

Choon et al., (2011) performed a computational fluid dynamics optimization anal- 
ysis on a Tesla turbine with water as working fluid with the aim of exploiting the 
energy hold within the household water supply. The CFD analysis allowed the up- 
grade of an existing prototype, resulting in an increase in efficiency from 6.8 to 10.7%, 
with a pressure drop of 0.04 m. The main results were the possibility of applying this 
technology to such a niche application, as the expander is suited to produce power at 
very low pressure drops. 

Cirincione N, (2011) designed and realized an ORC waste heat recovery system 
with a Tesla—hybrid turbine, which was claimed to be able to reach isentropic effi- 
ciency levels above 70% with steam; the last not reported tests were carried out with 
R245fa. 

Emran TA, (2011) exploited and further developed the results obtained in (Emran 
TA, 2010), improving the torque model of a Tesla turbine, as well as comparing the 
analytical solution with experimental tests. As a further result, the developed model 
allowed the design of an improved Tesla turbine with optimal dimensions. 

Ho-Yan BP, (2011) applied the model developed in (Rice W, 1965) in order to es- 
tablish a preliminary design approach for Tesla turbine for Pico Hydro applications. 
The obtained performance was over 8096 for a 300 W power output turbine, working 
with a pressure head of 20 m and a water flow rate of 2.5 l/s. 

Krishnan V et al., (2011) realized some micro-turbines (1-2 cm diameter rotors) 
using commercially available technologies, and tested them with different nozzles and 
rotors configurations, achieving almost 40% shaft mechanical efficiency. 

Lampart P & Jedrzejewski L, (2011) further developed the analysis presented 
(Lampart et al., 2009), considering the Tesla turbine as the expander for a micro- 
power plant of 20 kW heat capacity. The maximum efficiency obtained from a com- 
putational fluid dynamic analysis was 51% with a mass flow rate of Solkatherm 
SES36 of 0.13 kg/s. 

Podergajs M, (2011) provided a short report on Tesla turbine historical back- 
ground, flow model and applications. In the report, the author highlighted the range 
of application where the turbine could be an attractive solution. 

Romanin VD and Carey VP, (2011) developed an integral solution scheme of the 
fluid dynamics of incompressible, steady flow inside a Tesla turbine rotor. Enhanced 
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drag effect, due to surface micro structuring was also taken into account. The enhance- 
ment of efficiency due to surface micro grooving could be effectively relevant, as high 
as 9% enhancement when compared to smooth surface. The maximum predicted ef- 
ficiencies of the turbine were over 80%. 


2012 


Borate HP and Misal ND, (2012) carried out a performance analysis of the Tesla 
turbine, with a specific focus on the effect of surface finish and the space between the 
disks. The assessment was carried out considering water as working fluid. Experi- 
mental and theoretical analysis allowed the characterization of the major inefficiency 
(from the nozzle) and that, if direction grooving on disk surface is applied, an increase 
in efficiency between 5 to 6% can be achieved. 

Peshlakay A, (2012) performed a numerical and experimental analysis of a Tesla 
turbine, comparing different nozzles using air, water and steam as working fluids; 
achieving a rotor efficiency of 95% (+ 9.5% uncertainty) and a global turbine effi- 
ciency of 31%. 

Romanin VD et al., (2012) modified the model developed in (Romanin VD and Carey 
VP, 2011) in order to match the results obtained in (Krishnan V et al., 2011). The modi- 
fied analytical model results were compared with CFD analyses and experimental re- 
sults. The correlation of the results was fairly good, as, on average, the experimental 
results efficiency was 45% lower than the expected analytical model results. The dis- 
crepancy was mainly due to the not proper modelling of the nozzle. 

Romanin VD, (2012) presented a comprehensive study, which collected analytical 
perturbation model results, computational fluid dynamic investigation and experi- 
mental analysis. The obtained results were presented in (Romanin V et al., 2010, Krish- 
nan V et al., 2011, Romanin VD and Carey VP, 2011, Romanin VD et al., 2012). 

Sengupta S and Guha A, (2012) presented an analytical model resulting from the 
reduction of Navier-Stokes equations. They assumed a steady, incompressible and 
laminar flow, introducing the viscosity of fluid and a velocity gradient near the walls. 
This mathematical theory represents a simple but effective method of predicting the 
performance and efficiency of a Tesla turbine. 

Van Wageningen T, (2012) presented a concept review of the Tesla turbine to be 
utilized as the engine for a flapping wing mechanism. Tesla turbine seemed not to be 
suitable for this application when compared to two stroke and heat engine. 


2013 


Bao G et al., (2013) defined and validated a numerical model using CFD tools to 
describe the flow boundary layer; then the model was applied to different organic 
working fluids to obtain the related performance curves, concluding that best perfor- 
mance can be achieved with thin gap width, turbulent flow and fluids with high kine- 
matic viscosity. 

Deng Q et al., (2013) proposed and improved the model developed in (Guha A and 
Smiley B, 2009) and compared the analytical formulation with computational fluid dy- 
namic simulations. Various nozzle configurations were taken into account. The high- 
est efficiency (43.6%) of a 0.05 m rotor diameter turbine was reached when the lowest 
number of nozzle (2) was taken into account, as well as the highest rotational velocity 
(24,658 rpm) and the largest pressure ratio (0.8). 
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Guha A and Sengupta S, (20/3) further polished the model presented in (Sengupta 
Sand Guha A, 2012) and moreover they investigated, for the first time, the roles of each 
force affecting the power and pressure fields, physically explaining several fluid dy- 
namics behaviours like flow reversal and the shape of complex relative path lines. 

Guimaraes LNF et al., (2013) considered the Tesla turbine as a possible concept 
to be applied to an Emergency Core Cooling Systems of a standard power plant. They 
performed an experimental campaign with air as working fluid, reaching a 55 W 
power output value. 

Gupta HE and Kodali SP, (2013) presented a state of the art review of the Tesla 
turbine. The assessed performance, as well as the other main results of investigation 
of previous research projects were highlighted, taking into account also the patents 
developed. A special focus has been put on the assessment of possible applications of 
Tesla turbines. 

Khan MUS et al., (2013a) developed a modified configuration of a Tesla turbine, 
which merged both the aspect of a bladeless and a bladed turbine. An experimental 
campaign with water as working fluid was carried out. The average power output of 
the turbine was of 100 W. This new turbine concept was developed in order to exploit 
the potential energy of water flows present in canals or rivers. In a successive work, 
Khan MUS et al., (2013b) explored the possible applications where the Tesla turbine 
could be utilized. Particularly, the presented assessment of applications of a Tesla 
turbine was very wide: from power generation coupled with renewable energy 
sources, to the utilization in irrigation channels, up to the use in hybrid electric vehi- 
cles. 

Krishnan VG et al., (2013) described a method to scale Tesla turbines to millimetre 
dimensions, also providing the assessment of loss mechanisms. 4096 efficiency was 
expected even when millimetre scale turbine was considered. 

Placco GM et al., (2013) considered a modified Tesla turbine as the expander for 
a passive thermal circuit as a part of an Emergency Core Cooling System of nuclear 
power plant. An experimental campaign with air as working fluid was conducted, 
achieving power higher than 100 W with low mass flow rates. 

Sengupta S and Guha A, (2013) investigated the three-dimensional flow field and 
the flow path lines within a Tesla disk turbine, comparing the results obtained from 
the analytical theory and the computational fluid dynamics; moreover, they investi- 
gated the operating parameters that affect the shape of the path lines within the rotor 
and the local balance of the various forces. 

Yang Z et al., (2013) described an innovative method to measure and predict the 
mechanical power output of a Tesla turbine. The method is based on the rotational 
inertia of the turbine, as well as on the friction in the bearings. It is called “dynamic 
dynamometry" and allows the realization of power curves as function of rotational 
speed. 

Zhao D and Khoo J, (2013) designed and tested a 40 mm bladeless turbine for the 
harvesting of energy from air and rainwater applications. 0.5 W power output was 
obtained with air as working fluid at 3300 rpm. 


2014 


Guha A and Sengupta S, (2014a) developed a similitude study on the flow of the 
Tesla turbine. The scaling laws were obtained through the Buckingham Pi theorem, 
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which lead to the definition of 7 fundamental non-dimensional numbers (rotor radius 
ratio, aspect ratio, tangential speed ratio at rotor inlet, flow angle at rotor inlet, dy- 
namic similarity number, power coefficient and pressure drop coefficient). A further 
study of Guha A and Sengupta S, (2014b) demonstrated that the application of the 
Euler turbomachinery equation is consistent only if local velocity mass—averaged val- 
ues are considered. 

Hasan A and Benzamia A, (2014) investigated the effects of temperature and mass 
flow rate variations on the performance of a Tesla turbine of 100 mm outer rotor di- 
ameter, through the means of CFD analysis. The performance of the turbine was found 
to be dependent on the temperature, as the boundary layer thickness changes due to 
temperature effects. 

Pandey et al., (2014) developed a computational analysis of a 1 kW Tesla turbine 
for pico hydropower applications. A 127 mm rotor outer radius was considered for 
the calculations, obtaining a total of 777 W for a 9 disks configuration of the turbine, 
reaching an efficiency of 77.7%. 

Schosser et al., (2014) carried out an assessment on the design and optimization 
of a test facility for 3D tomographic PIV measurements of the flow field inside a Tesla 
turbine rotor. A throughout CFD optimization analysis was carried out on the test rig 
components and specifically on the Tesla turbine rotor. 

Shimeles S, (2014) performed a design an optimization assessment of a Tesla disk 
turbine, applying an incompressible steady state computational fluid dynamic analy- 
sis. Various rotor configurations were simulated, changing both geometric and fluid 
dynamic conditions. The maximum predicted rotor efficiency was of about 27%, uti- 
lizing air as working fluid. The maximum total power output obtained was of over 6.5 
kW at 11500 rpm. 

Siddiqui MS et al., (2014) presented a computational fluid dynamic assessment of 
an air Tesla turbine. The turbine dimensions were selected in order to compare the 
results with experimental data provided in (Rice, 1965). Different flow configurations 
were analysed both applying laminar and К—є turbulence schemes. The maximum ef- 
ficiency obtained was over 20% with the lowest flow parameter assessed. Efficiencies 
of numerical and experimental results were comparable. 

Singh A, (2014) derived a closed form expression of the flow field and pressure 
distribution inside a Tesla turbine rotor by simplifying Navier-Stokes equations. A 
careful assessment of the flow behaviour as function of the Reynolds number was 
carried out. 

Thawichsri K and Nilnont W, (2014) performed an experimental investigation on 
two Tesla turbines, differing in rotor outlet diameter (the first one of 120 mm, the 
second one of 75 mm). The working fluid expanding through the turbine was Isopen- 
tane. Three different hot source temperatures were analysed (70, 80 and 90 °C). For a 
0.05 kg/s mass flow rate and a rotational speed of 3000 rpm, the total effectiveness 
was of 36%, with higher power output given by the bigger turbine. 

Zhao D et al., (2014) performed an experimental campaign, as well as a computa- 
tional design assessment, of three Tesla turbines for the harnessing of energy from 
rainwater. Firstly, an experimental campaign with air as working fluid was conducted; 
finally, a rainwater application was assessed. The overall energy conversion assessed 
was approximately between 2 and 3%. 
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2015 


Baginski P and Jedrzejewski L, (2015) focused their studies on the assessment of 
the dynamic analysis of the Tesla turbine. Both strength and modal analysis were car- 
ried out through the means of two commercial software (Abaqus and Ansys). The 
results obtained for the assessed geometry (the same one used in (Lampart P and 
Jedrzejewski L, 2011), which consisted of a rotor composed by 11 disks, with 100 mm 
outer diameter and 2 mm thickness), showed that no resonance effect would present 
at start-up and that the first mode shape presented at 387 Hz. 

Guimaraes LNF et al., (2015) deepened the work presented in (Guimaraes et al., 
2013, Placco et al., 2013), developing a Tesla turbine for space application, claiming 
that this technology could be utilized both in Brayton and Rankine cycle, with the 
optimal characteristics which make it suitable for space application, such as low pres- 
sure drop and almost no maintenance required, as well as the possibility of using it 
with any kind of fluid. 

Holland K, (2015) designed, built and tested a 92 mm rotor diameter Tesla turbine. 
The tests were carried out with air as working fluid and the maximum obtained effi- 
ciency was of 8.5% at 3 bar inlet total pressure and 10,000 rpm. 

Kólling A et al., (2015) applied the Tesla turbine concept to small-distributed 
power generation from biomass resources. Particularly, in the study, the design, con- 
struction and tests of a 175 mm rotor diameter Tesla turbine were presented. The max- 
imum power obtained was of 1.4 kW when a saturated steam at 9 bar and 0.11 kg/s 
mass flow rate was utilized. The assessed isentropic efficiency was of about 30%. 

Krishnan V, (2015) exploited the work presented in (Krishnan V et al., 2011, Romanin 
VD and Carey VP, 2011, Krishnan V et al., 2013 ) providing a Matlab tool for generating 
the design of micro to small Tesla turbines. Analytical, computational and experi- 
mental results were obtained, demonstrating the soundness of concept of Tesla turbine 
for sub- Watt and Watt range applications. 

Neckel AL and Godinho M, (2015) realized and tested ten convergent-divergent 
nozzles to improve the injection efficiency of the working fluid; different geometries 
were first investigated with one and two-dimensional approaches and then an exper- 
imental test campaign with air as working fluid was carried out. 

Nedelcu D et al, (2015) performed a numerical simulation both on stress evalua- 
tion and the flow calculation of the Tesla turbine rotor. The aim was to develop and 
manufacture a turbine working with air for teaching application. 

Raje a et al., (2015) presented a review of the principle of operation, design con- 
sideration, performance parameters and possible applications of the Tesla turbine. 

Ruiz M, (2015) applied the Tesla turbine concept to a completely different appli- 
cation, which does not involve power production, but instead utilized as a heat sink. 
Numerical as well experimental results were obtained both with single phase and two— 
phase conditions (water as working fluid). The milestone to be exploited for Tesla 
turbine power applications is the development of an accurate heat transfer model in 
two-phase flow conditions. 

Ruiz M and Carey VP, (2015) presented the experimental results obtained when 
utilizing a Tesla turbine as a heat sink, employing the work developed in (Ruiz, 2015). 
Heat transfer, as well as pressure drop of a micro tesla turbine with water as working 
fluid was analysed. 
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Schosser C and Pfitzner M, (2015) developed a throughout CFD analysis, focusing 
on the velocity profile of the flow inside a Tesla turbine with air as working fluids. 
Particularly, they stated that laminar CFD results were better approximated by a fourth 
order polynomial function, compared to parabolic profile. 

Thawichsri K and Nilnont W, (2015) performed an experimental comparison be- 
tween a centrifugal turbine and a Tesla turbine working with Isopentane at low tem- 
perature (70—90?C) heat resources. The experimental assessment demonstrated, from 
one point that the centrifugal turbine was on average 3096 more performing than the 
Tesla turbine, but also that the Tesla turbine is a very cheap and easily to manufacture 
turbine. The maximum reached efficiency of the Tesla turbine was of about 1295, the 
specific power output of 35 kJ/kg and the cycle efficiency of 8.696, when the heat 
source was 90?C. 


2016 


Bankar N et al., (2016) designed a Tesla turbine for micro power application. Par- 
ticularly, they called the new design approach as “hybrid Tesla turbine", due to the 
introduction of some changes, such as the grooving of the disks and utilization of 
polycarbonate disks of thin thickness. 

Herrmann-—Priesnitz В et al., (2016) developed a fluid dynamics model on the 
structure of the boundary layer of a Tesla disk turbine. Inviscid core region was cou- 
pled to boundary layer formulation in order to develop a sound analytical model, 
which was compared to numerical simulations run in OpenFoam environment. 

Jose R et al., (2016a, 2016b) realized a comprehensive project on the design, the- 
oretical and experimental analysis of a Tesla turbine utilizing water as working fluid. 
The effect of surface roughness, as well as disk spacing and number of disks was 
analysed. The experimental results highlighted that a spiral grooving on disks could 
enhance the turbine efficiency of about 5-6%. 

Joshi KN et al., (2016) coupled the Tesla turbine design configuration with a Pel- 
ton turbine, obtaining therefore a “hybrid Tesla wheel turbine". A computational fluid 
dynamic assessment on Tesla and hybrid Tesla Pelton turbine was carried out. The 
Hybrid turbine allowed a 5% increase in efficiency compared to the bladeless turbine. 
A rapid prototyping technique was also showed as a possible solution for Tesla-Pelton 
turbine fabrication. 

Qi W et al., (2016) performed a numerical investigation on the influence of disk 
spacing on the performance of Tesla turbines. The analysis conducted considered 7 
different spacing (from 0.1 to 1 mm) highlighting the changes in performance (power 
and efficiency), as well as on the mass flow rate of the turbine. Several rotational 
speeds were also considered and optimal values of disk spacing, as well dimensionless 
tangential velocities were obtained when air as working fluid was considered. 

Schosser C et al., (2016) performed a throughout investigation of the flow field 
inside a Tesla turbine with air as working fluid, through the means of a stereoscopic 
3D-PTV measurement technique. The main result was the assessment of the laminar 
velocity profile inside the channels, which slightly differs from parabolic distribution 
and is better described by a 4" order polynomial. 

Schosser C, (2016a) exploited the study presented in (Schosser C and Pfitzner M, 
2015, Schosser C 2016b) and carried out a throughout experimental and numerical in- 
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vestigation of a Tesla turbine working with air as working fluid. Analytical and nu- 
merical results were also compared and an efficiency prediction of over 50% was 
obtained. 

Sengupta S and Guha A, (2016) investigated the performance of a Tesla turbine 
when nanofluids were utilized as working fluid; an increase of power output of 30% 
appears to be possible when the volume fraction of nanoparticles is increased from 0 
to 0.05. Therefore, the maximum predicted efficiency was of about 56% when water 
based nanofluid with ferro—particles were utilized. 

Zahid I et al., (2016) analysed the performance of an air driven Tesla turbine. The 
claimed reached efficiency was of about 58% for a turbine with an external rotor di- 
ameter of 100 mm, run with an inlet pressure of 272 kPa. 


2017 


Alrabie MS et al., (2017) applied an analytical solution in order to assess the per- 
formance of a small-scale hydro Tesla turbine. The design of experiment was carried 
out and the definition of the optimal prototype geometry was obtained, assessing var- 
ious geometric parameters (inlet, outlet rotor diameters, disks spacing, nozzle number 
and location), as well as the most significant thermodynamic parameters (density, vis- 
cosity). Finally, an 80 mm rotor outer diameter turbine was selected as the proper 
design in order to obtain 0.47 W power output at 1197 rpm. Similarity scaling law 
equations were applied to assess the performance of a 3 mm rotor outer diameter tur- 
bine working with Ethylene Glycol. 

Damodhar R et al., (2017) designed and fabricated a cm-scale Tesla turbine uti- 
lizing water as working fluid. Design specifications were given; as well experimental 
tests were conducted. The maximum power obtained was of 0.12 W at 140 rpm, with 
a corresponding efficiency of 20.97%. 

Guha A and Sengupta S, (2017) performed an optimization assessment ofthe Tesla 
turbine based on the non-dimensional scaling laws presented in (Guha and Sengupta, 
2014). The role of each non-dimensional parameter was assessed and the physical 
explanation of its effects on the performance of the turbine is given. As an illustrative 
example of the methodology applied, a turbine efficiency of 68.8% was designed 
when a tangential inlet velocity ratio (the ratio between tangential velocity and rota- 
tional speed at rotor inlet) of 1.45 was considered. 

Li et al., (2017) explored the performance characteristics of the Tesla turbine 
through experimental and numerical assessment. An 11 mm outer rotor diameter, with 
1 mm disk spacing and a total of 13 disks, Tesla turbine was fabricated and tested with 
water as working fluid. The experimental assessment was carried out at various, rota- 
tional speeds, inlet pressures and mass flow rates; a maximum efficiency of about 896 
was obtained. The numerical simulation assessed efficiencies higher than 4096 for 
rotational speeds between 700 and 1000 rpm. 

Lisker R et al., (2017) developed an analytical model, which takes into account 
thin film condensation in a Tesla turbine rotor. The influence ofthin film condensation 
on velocity fields and mass flow rate was highlighted when condensing steam was 
utilized as working fluid. 

Mandal A and Saha S, (2017) investigated the performance of a cm-scale Tesla 
turbine for micro-air vehicles application. Numerical simulations on 20 disks, 0.1 m 
outer rotor diameter air Tesla turbine were performed. Different rotational velocities 
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were assessed, obtaining maximum rotor efficiency of 89.56% for an inlet velocity of 
10 m/s and a maximum power output of 21.3 W for an inlet velocity of 20 m/s at 1500 
rpm. 

Polisetti S et al., (2017) fabricated and conducted a performance study on an air 
driven Tesla turbine. A 95 mm outer rotor diameter, 1.27 mm thick disks Tesla turbine 
was fabricated and various configurations were tested, varying number of disks, dis- 
tance between the disks, number of nozzles and pressure at rotor inlet. The results 
were presented in the form of rotational speeds. The maximum reached rotational 
speed was 25,324 rpm for a 4 disks configuration at 6 bar inlet pressure. 

Schosser C et al., (2017) presented a comparison between analytical and numerical 
models on the flow field developing inside a Tesla turbine single channel. Several 
analyses were conducted, showing the influence of non-dimensional parameters, such 
as dimensionless friction factor, inlet conditions and geometry parameters. Inflow ef- 
fect was analysed, as the discrepancy between theoretical and numerical analyses was 
found to be relevant. 

Shah V and Dhokai S, (2017) carried out experimental tests of an air driven Tesla 
turbine. The turbine was composed by 8 disks of 95 mm rotor outer diameter, with 
disk spacing of 1 mm. The maximum power obtained was of 137.83 W at 4500 rpm, 
with 4 bar of total inlet pressure. 

Shah V et al., (2017) conducted a literature review on the principle of operation 
and on the parameters affecting the performance of the Tesla turbine. 

Song J and Gu CW, (2017) developed a 1D model of a Tesla turbine working with 
organic fluids. Furthermore, a thermodynamic analysis of an organic Rankine cycle 
utilizing R245fa as working fluid, exploiting the Tesla turbine as the expander was 
carried out. The assessed isentropic efficiency of the turbine reached a maximum 
value of 45% when a 15° oriented nozzle was utilized and an evaporation temperature 
of 370 K was taken into account. The maximum power output of 1.24 kW was ob- 
tained when the heat source temperature was of 355 K. The maximum thermal cycle 
efficiency obtained was of about 8%. 

Song J et al., (2017) improved and used the one-dimensional Tesla turbine model 
presented in (Song J and Gu CW, 2017) to predict the efficiency of a small scale ORC 
power plant adopting various working fluids and operating conditions: at design point, 
the ORC with R245ca released 1.25 kW power output at 4% thermodynamic effi- 
ciency. 

Thiyagarajan V et al., (2017) applied the Tesla turbine concept to solar power 
refrigeration. A prototype of 50 mm outer rotor radius was fabricated and some pre- 
liminary tests carried out. Based on the data obtained, the theoretical time to convert 
1 kilogram of water in ice was calculated (6.5 hr). 

Umashankar M et al., (2017) assessed the possible application of the Tesla turbine 
in cogeneration of heat and power systems (CHP). Cogeneration concept and turbine 
analytical approach were assessed. 

Variava JM and Bhavsar AS, (2017) performed an experimental and numerical 
analysis on Tesla turbine utilizing water as working fluid. The experimental investi- 
gation allowed assessing turbine efficiency as high as 6096 with a 0.56 kg/s mass flow 
rate and a power output of 15.18 W. 
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2018 


Sengupta S and Guha A, (2018) investigated the effect on non—uniform admission 
inside a Tesla turbine rotor, taking into account discrete inflows effect, finite disks 
thickness and radial clearance. An extensive computational fluid dynamic assessment 
was catried out assessing both the variation in power and efficiency as function of the 
number of nozzles, clearance, rotational speed and disks thickness. Small disk thick- 
ness, straight disk edge, optimum radial clearance and high number of nozzle were 
recommended as an optimal design solution. 

Sengupta S and Guha A, (2018) improved the one-dimensional model developed 
in (Song J and Gu CW, 2017, Song et al., 2017) and compared the predicted performances 
with the experimental results obtained by Rice, (1965) with air as working fluid. 

Traum MJ et al., (2018) extended the work developed by Carey VP, (2010) in order 
to assess the sensitivity of shaft power to a selected design variable. A differentiable 
closed form analytical model was applied, starting from the one developed in (Carey 
VP, 2010) and applying some substantial changes, such as: identifying model Reynolds 
number limitation; defining closed-form torque and power output expressions; intro- 
ducing simple nozzle geometric height parameter and showing the importance of 
checking maximum turbine rotational speed. 


SERG (Sustainable Energy Research Group) contribution to Tesla turbine 


state of the art 


The SERG group of the department of Industrial Engineering of University of 
Florence contributed to the state of the art on Tesla turbines with the following re- 
search. 

2017 

Manfrida et al., (2017) revisited the original concept presented by Tesla, improv- 
ing the stator layout and applying a modular design to the turbine, increasing its flex- 
ibility. The evaluation of the turbine performance with two different organic working 
fluids (R245fa and n-Hexane) was also carried out, analysing the influence of the 
most significant turbomachinery non-dimensional parameters (ns, ds, ф, y). A prelim- 
inary design of a prototype was also presented, and the expected efficiency with the 
assessed working fluid highlighted. Maximum efficiency (about 40%) was obtained 
for n-Hexane as working fluid, for low mass flow rates. 

2018 

Manfrida G and Talluri L, (2018) revised the mathematical model of the Tesla 
turbine, upgrading the model presented in (Carey VP, 2010, Sengupta and Guha, 2012), 
taking into account compressible real fluid properties. A performance analysis on the 
Tesla turbine rotor was presented, allowing the understanding of the principal param- 
eters influencing turbine efficiency and reaching a design of an air Tesla turbine rotor 
with over 90% efficiency. 

Manfrida G et al., (2018) further polished the work conducted in (Manfrida et al., 
2017), presenting an upgraded concept of the Tesla turbine for Organic Rankine cycle 
applications. Concept design and performance assessment of a manufactured proto- 
type were performed. The fluid assessment comprehended 5 different organic work- 
ing fluids: R245fa, R134a, SES36, n-Pentane, n-Hexane. High turbine efficiency 
(over 5596) can be reached for n-Hexane as working fluid. Higher power production 
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was reached with refrigerant fluids, and SES36 resulted as the best compromise be- 
tween power production and efficiency. 

Talluri et al., (2018) performed a design and optimization of Tesla turbine for 
ORC applications. An innovative rotor model was presented, as well as the assessment 
of each component losses was carried out. Not only the rotor, but also stator, still 
chamber and diffuser model were designed. As an outcome of the methodology pro- 
cess, three different turbine configurations were assessed; the maximum total to static 
efficiency assessed was of 64% when a 0.5 m rotor diameter was considered and n— 
Hexane was utilized as working fluid. The power production per single channel of 
such configuration was of 58 W, which allowed a very compact axial design. 
Statistics 

To conclude the literature review, a resume of the available literature is presented 
in Tab. 2.6. The previously described works were resumed highlighting the type of 
analysis that was carried out and which kind of fluid was considered. 

Fig. 2.18 displays the number of publication on Tesla turbines (updated in August 
2018). It has to be remarked how the trend is drastically increasing in the last few 
years; this is directly linked to the high attention that micro power generation gained 
on the energy market. Fig. 2.19 collects the most active authors assessing the Tesla 
turbine. Rice W (Arizona State University, Tempe) is the author that most investigated 
the Tesla turbine, both from analytical and experimental point of view. Of particular 
relevance are also the work developed by Guha A and Sengupta S (Indian Institute of 
Technology, Kharagpur) and Carey VP (University of California, Berkeley), who de- 
veloped analytical, computational and experimental research, properly explaining the 
flow behaviour inside a Tesla disk turbine. Finally, Fig. 2.20 resumes the working 
fluids that have mainly been considered when developing a study on Tesla turbine. 
Air is the principal fluid investigated, but of particular importance is the emerging 
interest in recent years of the application of Tesla turbines in organic Rankine cycles. 
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Fig. 2.18 Number of publications (Journal paper, conference paper, Ph.D. and M.Sc. available 
thesis) from 1950 to nowadays (updated in August 2018) 
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Fig. 2.19 Most active authors on Tesla turbine research 
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Fig. 2.20 Working fluid taken into account in literature studies 
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Chapter 3 
Methodology and Models 


As discussed in Section 2.4.1, the first concept of the Tesla turbine was developed 
by Tesla in 1913 (Tesla, 1913). Conversely, to traditional turbomachines, which exploit 
the pressure difference that is produced when a fluid flows around a row of blades; 
the Tesla turbine generates power through the frictional interaction between the evolv- 
ing fluid and the bladeless rotor. This peculiar rotor, characterized by the absence of 
blades, is the main feature of the Tesla turbine. Indeed, due to the configuration of the 
rotor, as well as to the principle of operation, the Tesla turbine is also called bladeless 
turbine, viscous turbine, boundary layer turbine or friction turbine. 

Differently from conventional turbines, the rotor is composed by a series of par- 
allel flat disks with a very small gap between them. The admission of the flow in the 
rotor occurs through one or more nozzles, which allows the fluid to enter from the 
external radius of the disks and to exit from the openings made on the disks at the 
inner radius. Inside the rotor, the fluid depicts a spiral centrifugal path, due to the 
interaction between the viscous forces and exchange of momentum. Typical configu- 
rations of the Tesla turbine consider only the rotor as part of the turbine, and often, 
nozzles are just converging pipes positioned tangentially to the rotor outer radius, as 
displayed in Fig. 3.1. 
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Fig. 3. 1 Typical configuration of Tesla turbine a) (Armstrong, 1952), b) (Rice, 1965), c) (Stei- 
del and Weiss, 1976), d) (Carey, 2010) 
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Conversely to the typical configurations, the Tesla turbine described and analysed 
in this work consists of several components: an external toroidal plenum chamber, a 
stator with fixed nozzles and a bladeless rotor composed by parallel thin disks fixed 
to the rotating shaft; Fig. 3.2 displays the schematic of the ORC prototype, nonetheless 
the same nomenclature is used for the air prototype. 


Stator 


Fig. 3. 2 Schematic of Tesla turbine 
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1. Objectives and structure 


A 2D thermo-fluid dynamic model was developed in Engineering Equation 
Solver (Klein and Nellis, 2012). Each procedure is thus explained. 


1.1 Nozzle Design 


Blade design is a procedure, which enables the definition of the stator geometry. 
The stator is designed like a circular crown, of a certain thickness “Hs” in the 
axial direction, on which a series of channels (“ог nozzles") are obtained, defined by 
the blades profile. 
The required inputs for the procedure are: 
e Outlet diameter of the stator (r7); 
e Inlet and outlet angles (ao, a); 
e =Total length of the camber (mot), which can be determined imposing the in- 
let/outlet diameter stator ratio at 1.25 as suggested in (Glassman, 1976, Ventura 
et al., 2012). 


Another feature of the code is that the camber angle distribution can be adjusted 
changing the exponent n of a predefined power law. 


e train 
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Fig. 3. 3 Blade design code interface 


Once all the inputs are fixed, the code draws blades on the circular crown, defining 
the main parameters such as width of the throat section, chord, pressure and suction 
side exit angles). The maximum feasible number of passages is then calculated (Zmax) 
and partialized (Z;). The partialization consists of the occlusion of several volumes 
present between two blades, due to the low mass flow rate required by the turbine 
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rotor. The final geometry of the stator is then obtained, composed by a small number 
of real channels, located inside a thick volume. 

Specifically, the program can include the coordinates of NACA or any other blade 
profile, allowing the construction of pressure and suction sidewalls. The user can 
adapt the geometric inputs and the drawn profile will be visible both in the configu- 
ration of a single blade and as section of stator, as shown in Fig. 3.3. 

Blade design main equations 
In the lookup table, the coordinate of the selected NACA profile with chord value 
equal to 1 are stored. Absolute camber line values are thus determined: 
mj = (т. E Mr(i-1)) * Mtot (3.1) 
Where: 
m; — absolute discretised value of camber line 
m; — relative discretised value of camber line 
mio: — total length of camber line 

The camber line is then curved through an exponential law, which sets the local 
camber angle. 

а = Ag + (ag — 04): mp (3.2) 

In order to find the coordinate of each point of the camber line, trigonometry anal- 
ysis is applied. From the law of cosines, it is possible to determine the radius and from 
the law of sines the wrap angle (Fig 3.4). 

ri = (11-1)? + (mj)? — 2: mj ri, * cos(a;)) (3.3) 


Ө; = 0;_, + arcsin (=. sin(0;)) (3:4) 


Fig. 3. 4 Trigonometric representation of stator blade 
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It is then possible to find the Cartesian coordinate of the camber point with a sim- 
ple coordinate transformation. 
X; = rj;sin(0j) (3.5) 
Y; = rj: cos(6;) (3.6) 
The chord and stagger angle (the angle between the chord and radial direction) 
values are thus found: 


(Xo-X1) 3.7 

B — arctan (So) (3.7) 
— Or 3.8 
Chord = =) (3.8) 


The coordinate of suction and pressure side can be easily determined knowing the 
blade thickness (which is a function of the chord): 


: (3.9) 
Xipsyss = Xi + (Ce. costa) 
aE (3.10) 


Yi psyss = Yi + (Se sint) 


Once the coordinates of suction and pressure side of each blade are known, it is 
possible to determine the distance between two blades at each point and therefore to 
determine the throat length. 


1.2 Stator model 


As the most original component of a Tesla turbine is the bladeless rotor, simple 
nozzles are commonly used instead of properly designed stator vanes. The Tesla tur- 
bine stator purpose is to generate the necessary tangential flow stream at rotor inlet 
and to convert the pressure energy ofthe flow in the plenum chamber to kinetic energy 
at nozzles output. The reduction of cross sectional area for a subsonic flow produces 
a favourable pressure gradient and an acceleration of the fluid (Anderson, 2001), avoid- 
ing wall separation; as a result, the efficiency of nozzles is usually very high, often 
exceeding 96% (Shepherd, 1956, Shames, 2002). Anyway, for small size nozzles, where 
the throat width is lower than 3 mm (as is the case for Tesla turbine), the boundary 
layer might occupy a significant portion of the cross sectional area (Butenko et al., 
1976), generating increased viscous losses. In these cases, the flow is laminar 
(Re« 10?) and the total pressure losses decrease with increasing the Reynolds number 
(Butenko et al., 1976). 

As discussed in several papers (Armstrong, 1952, Rice, 1965, Hoya and Guha, 2009, 
Guha and Smiley, 2009, Neckel and Godinho, 2015), the stator is commonly acknowledged 
to be one of the critical components of the Tesla turbine. It is indeed recognized 
among the main reasons for the low efficiency of the machine obtained in several 
experimental tests, as it is the source of high total pressure losses. Specifically, Guha 
and Smiley, (2009) carried out an experimental and computational study of the Tesla 
turbine nozzles and concluded that the stator is responsible for total pressure losses 
up to 35% of the total inlet conditions. 

In the present study, the stator design follows the approach derived from radial 
expanders vaned stators (Fiaschi et al., 2012, Fiaschi et al., 2014, Fiaschi et al., 2016, 
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Rohlik, 1975, Whitfield and Baines, 1990), also accounting for partial admission due to 
the reduced flow rates. In (Fiaschi et al., 2014), the design guidelines for a radial ORC 
turboexpander are defined using a zero dimensional model: an input dataset (thermo— 
fluid dynamic variables, dimensional and non—dimensional parameters) is provided 
by the designer, while the outputs of calculations are the basic geometry, velocity 
triangles and nozzle efficiency. 

In order to carry out a parametric analysis, a set of geometric and thermodynamic 
parameters, needs to be defined. The formers are: stator blade angles, number of noz- 
zles, length and height of the channel and nozzle geometry profile; which allow the 
definition of the full geometry (in particular the throat width and the chord length). 
The main thermodynamic conditions to be fixed are the inlet total pressure and tem- 
perature and the mass flow rate. 

Two iterative processes were implemented, the former one on mass flow rate, the 
latter one on two stator loss coefficients. The static pressure at throat section is taken 
as the iterative parameter for the matching of the mass flow rate, which will ultimately 
assume the imposed input mass flow rate. 

The first loss coefficient (Cw) was taken from (Fiaschi et al., 2016, Whitfield and 
Baines, 1990, Dixon, 2005) and it is only dependent on Ф, which is the ratio between 
the real and isentropic velocity at stator exit. The range of this parameter is between 
0.9 and 0.97. The second loss coefficient (Cr) (Ventura et al., 2012) is on the other hand 
dependent on geometry values, such as chord, pitch, stator height and stator exit angle, 
as well as on the Reynolds number of the flow. The two coefficients are compared 
and until they are not equal, an iterative process on flow exit velocity runs. 

Stator main equations 

The thermo-fluid dynamic model for the calculation of the fluid behaviour into 
the Tesla turbine stator assumes real fluid Equations of State (EoS). For this reason, 
all the thermodynamic properties were evaluated as functions of couples of local var- 
iables (typically, p and 7) using the Engineering Equation Solver (EES) EoS library 
data. Fig. 3.5 displays the enthalpy—entropy diagram of the Tesla expander stator and 
its nomenclature. 

The thermodynamic of point 00 is defined from input values poo and Too. Point 0 
is defined through iteration on density (first guess value is роо) and application of 
mass flow rate definition in order to find vo. 

m = Vg ' Po ` Ag (3.11) 

The root finding interpolation method was applied in order to find the correct 
expansion ratio in the stator. P; is taken as the variable for the interpolation method; 
convergence is reached when the imposed mass flow rate and the calculated one are 
equal, as shown in the diagram of Fig. 3.6. Eqns. (3.12) and (3.13) are main equations 
of the applied scheme. In order to set the right starting point, the upper bound for P; 
was assumed as equal to Po — 100 Pa; the lower bound was set assuming an incom- 


) 


: ) : (Р, (down) — Р, (ир)) (3.12) 
ир 


pressible flow expansion. 


(т — 


Р, (up) = P (up) +-———“ 
(Mabwn —m 
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Fig. 3. 5 Stator Enthalpy-Entropy diagram 


An iterative process was also implemented on two loss coefficients Су and Cr in 
order to assess the right velocity at stator outlet and therefore the efficiency of the 
nozzle. Eqns. (3.14) and (3.15) display the calculation of the loss coefficients. 


hy T his 
n= 3 | (3.14) 
=V. n 
ME 
0.05 tan(a,) | cos(04) 
ба = Reoz C" pitch TRU y (3.15) 
chord 
Where: 


H,- Blade height at stator outlet ; a,— Absolute angle at stator outlet; 
Re — Reynolds number (based on blade height and absolute velocity at stator out- 
let) 
After the two iterative process point 1 and 01 are defined and, finally, stator effi- 
ciency is calculated as: 
_ (oi ha) 


Е : (3.16) 
LI gi 
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Geometric: Н,, Zstat, 
ThroathWidth, a,, 
Chord, D,, Dy 


Geometry Definition 


Thermodynamic and 
Kinematic of Point 
1:18: 01:01 


Fig. 3. 6 Flow Diagram of Stator model 
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1.3 Stator-rotor coupling model 


Assumptions 

When connecting the three-separated developed flow models (stator, rotor and 
gap pressure losses), the following assumptions were applied: 

e Conservation of total enthalpy between the stator outlet and the rotor inlet, 
because the transformation may be considered adiabatic and without any 
work transfer. 

e Invariable static enthalpy, and consequently constant absolute velocity уу, 
while the flow direction changes during the gap crossing. It means that the 
pressure losses may be treated like an isenthalpic throttling process into a 
valve. 

Model characteristic 

The stator and rotor models only take into account of the distributed pressure 
losses inside the components, but they do not consider those concentrated into the 
stator-rotor gap. The passage of the flow from the throat of the nozzle to the gap first 
and to the rotor channel later, involve an abrupt cross-section enlargement followed 
by a contraction and generates large flow pressure losses. 

Both loss coefficient (abrupt enlargement and abrupt contraction) were obtained 
from standard incompressible flow loss treatment theory (Borda Carnot coefficient 
and polynomial fitting of empirical data (Idel 'chik, 1960)). In order to take into account 
compressibility effects, an iterative procedure was implemented in EES to compute 
average values of density between inlet and outlet sections. 

Stator—rotor losses main equations 

For an incompressible fluid, the concentrated pressure losses are calculated by the 
definition of a loss coefficient, depending on the system geometry and the flow con- 
ditions, which reduces the kinetic energy of the fluid (/del’chik, 1960): 


1 
Apo 76:5 pv (3.17) 
For the case of Tesla turbine, it can be written as: 


1 1 
Ap = Apen + APcon = 3 SenP Vi + 3 SconPW2 (3.18) 


Where Ape, is the pressure loss occurring immediately after the throat section 
(abrupt enlargement) and Apcon is the pressure loss related to the flow entering the 
rotor micro-channels (relative flow contraction), while Cen and Gcon are their respective 
loss coefficients. 

The loss coefficient for abrupt enlargement (Gen) was modelled as an incompress- 
ible Borda-Carnot coefficient (Idel’chik, 1960), according to Eq. (3.19): 


А, 
Cada (3.19) 
out 
Where A, is the throat cross section, here having rectangular shape: 
Ain = Ly: Hs (3.20) 


Aout is the cross section on the disks (transversal) covered by the flow jet: 
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tana, sina, 


]/cos a, — б : tana, — Hs (3.21) 


Aout = fr tan aay 


Eq. (3.19) can be used for a turbulent flow with a uniform velocity profile; when 
these assumptions aren’t checked, some numerical and graphical correlations availa- 
ble in (Idel’chik, 1960) can be used. 

The total pressure loss for abrupt expansion can be calculated using the velocity 
immediately upstream the enlargement. A parametric analysis showed that the pres- 
sure loss is mainly influenced by the velocity v; (second order law), followed by the 
density and, finally, by geometric parameters (throat width TW and gap extension G). 

The loss coefficient for abrupt contraction (боп) was obtained through a polyno- 
mial fitting of empirical data (Idel'chik, 1960): 


A At Гес Жоо 
бооп = 0.126 (=) + 1.0296 (==) — 1.279 (=) 
in X in in (3.22) 
— 0.12097") + 0.5 
Ain 

Where Ain is the total cylindrical surface, including disks and channel surfaces: 
Ain = 2107; : H; (3.23) 

While Aou only takes into account the cylindrical surface of the channel: 
Aout = 212: b: Neh (3.24) 


The velocity used in this case is the radial component of the relative velocity w,2 
(normal to passage section), immediately after the contraction. The main geometric 
parameters influencing the pressure loss are the height ofthe nozzle throat, the thick- 
ness and the number of rotor channels. Generally, the pressure loss for abrupt enlarge- 
ment is far higher than the one for abrupt contraction at rotor inlet, because the veloc- 
ity v; influencing the former is higher (typically by a factor 10). 

Eq. (3.18) could be used if the flow Mach number was less than 0.3, thus the flow 
could be treated as incompressible, but it is not applicable for the flow at nozzle out- 
put/rotor inlet sections. For this reason, two iterative calculations were implemented 
into the EES code in order to achieve average values of density between the input and 
the output sections (both for abrupt enlargement and abrupt contraction), which allows 
the recovery of compressibility effects. 


1.4 Rotor model 


The first developed model for the rotor flow was derived from (Rice, 1965, Carey, 
2010), applying some remarkable improvements. Specifically, real and compressible 
fluid behaviour was considered (rather than ideal and incompressible). The variable 
density and the other thermodynamic functions were taken as fluid properties, de- 
pending on the local variables (for example pressure and temperature). As for the sta- 
tor model, the fluid properties were locally evaluated using EES EoS library data. The 
assumptions of steady, laminar and two-dimensional flow were kept, as well as the 
viscous forces treated as body forces acting on the flow at each position (Carey, 2010). 
It allowed simplifying and thus numerically solving the fundamental Navier-Stokes 
equations, expressed in cylindrical coordinates. 
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Assumptions 
When developing the first rotor model, the following assumptions were made: 

a) Steady, laminar flow. 

b) The viscous force is treated as a body force acting on the flow at each 
(7—0) position. 

c) Two-dimensional flow: 

e р, = 0; 
e т, = constant across the channel ; 
e vg = constant across the channel. 

d) Radial symmetric flow field, uniform at the inlet (r — ro). The flow 
field is thus the same for any 0, therefore the derivative 0/00 = 0 for 
all flow variables. 

e) (0p/00) negligible compared to wall friction forces. 

General Flow Equations 
Taking into account the above assumptions, the fundamental Navier-Stokes equa- 


tions in cylindrical coordinates are reduced to: 


Continuity 
1 д(тру) _ 
va ru (3.25) 


Momentum, r-direction 
д 2 д 

vote = 2 (2) f (3.26) 
Momentum, Ó-direction 


v, Y ре = fy (3.27) 


дг г 


Momentum, z-direction 
1 (ӧр\ — 
-=(2) =0 (3.28) 


The integration of the reduced continuity Eq. (3.25) results in rpV, = costant. 
Furthermore, knowing the mass flow rate inside each channel, it follows that locally: 


Mc 


Vr = 


Е 2пгрр (3.29) 
Formulation of the viscous shear stress 

Considering a fluid element between the two disks defining the flow channel, a 
control volume V, can be defined with base surface А, and height b. The fluid wetted 
area is Ay, = 2A,. Therefore, the hydraulic diameter Оһ is equal to 2b. Consequently, 
Ve 2Ve 4Ve 


Ae TUB mu Aw b (3.30) 
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For laminar flow, the wall shear effect can be expressed as a function of a friction 
factor ¢ and of the relative velocity of the flow. Eq. (3.31) displays the expression of 
the wall shear stress, decomposing the relative velocity in its two components. 


t= ew? = 2 [(vo — wr)? + v,?] (3.31) 
Considering U = (Uo / rg): r and ( — 24 / Re as usual for laminar flow 
between parallel plates: 
_ 24u | 24p — 24u 
S PWDn o PWDh o pDny(ve- œr)?+vr? (3.32) 
So that: 


ту = m Ve — wr)? + у,2 (3.33) 


The force resulting from wall friction force is given by the product of the wall 
shear with the wetted area: 


F = ZEL We — en? + у? (3.34) 


The wall friction force has a tangential and a radial component, which influence 
the torque and the radial pressure gradient, respectively. 
Solution of the rotor flow 

Fig. 3.7 shows the local velocity triangle of the fluid element inside the rotor. 


Wo В 


тҮҮ: 


Fig. 3. 7 Local velocity triangle 


The radial component of the friction force is given by: 
Е, = F cos(p) (3.35) 


Where f is the angle between relative velocity and the radial direction. The value 
of cos(6) can thus be calculated as: 
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Wr | Vr 
соѕ(В) = vira (3.36) 
Substituting Eq. (3.36) in Eq. (3.35), a compact expression of the radial force 
component is obtained: 


Е, = VE v, (3.37) 


Dividing Eq. (3.37) by the mass of the fluid element between two disks, the body 
force term in the radial direction can be expressed as: 


= =; (3.38) 


Proceeding іп ће same way for the tangential direction, the wall friction force is 
given by: 


Fo = —F sin(B) (3.39) 
i _ We _ (vg- wr) 
sin) = 2 = ae (3.40) 


Similarly, substituting Eq. (3.40) in Eq. (3.39), a compact expression of the tan- 
gential force is obtained: 


Fy = — (ve - or) (3.41) 


The body force in tangential direction is thus given by: 


3 
fg = =r — wr) (3.42) 


In order to determine the local pressure, Eq. (3.38) is substituted in Eq. (3.26): 


OV; vg? 1 (2) 12u 
к же ш Sucre 4 
Уг or r р \дг + pb2 Ы: (3 3) 


Using Eq. (3.29), the local derivative (дъ, )/дт can be expressed as: 
-=--v (3.44) 


Finally, substituting Eq. (3.44) in Eq. (3.43), the pressure gradient in radial direc- 
tion is given by: 


Op — 12u( п, p( me Y р 2 

(2) 5 b? (=<) i r С) T r Ve (9:39) 
Likewise, in order to compute the tangential velocity, Eq. (3.42) can be substituted 

in Eq. (3.37): 


a r 12 
VPE = E — or) (3.46) 
Obtaining finally: 
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Ove 24utr(ve— or) ve 
or = ^ bn — X ps (3.47) 
Which determines the profile of vg(7). 

Eq. (3.47) was implemented in EES environment and numerically solved by ap- 
plying a step forward method (centered finite difference): the rotor channel was dis- 
cretized in radial direction with a predefined number of equal steps. Increasing the 
discretization steps allowed a higher accuracy prediction of the local thermodynamic 
variables; nonetheless, the increase of accuracy was obtained at the cost of an in- 
creased computational time. Therefore, two—hundred discretization steps were se- 
lected, as it was found as a reasonable compromise between accuracy of results and 
calculation time. 

This equation set allows the calculation of the local values of pressure and veloc- 
ity, both in absolute and relative coordinates. Finally, the rothalpy conservation Eq. 
(3.48) was applied to calculate the local value of static enthalpy: 

w^ u 


hah = tS (3.48) 


A further upgrade to the model, compared to (Rice, 1965, Carey, 2010) was the as- 
sumption of variable viscosity in Eqns. (3.45, 3.47): it was locally evaluated as a flow 
property, as a function of temperature and pressure, using the EES fluid library data. 
The results obtained showed a limited reduction of the viscosity values (lower than 
2%), because of the modest variation in fluid temperature. However, this upgrade al- 
lowed more accurate results at negligible additional calculation time. 

The rotor model was completed by the calculation of the performance indicators 
like power output, total to static efficiency and the non-dimensional fluid dynamic 
parameters, such as load and flow coefficients. 

Upgraded Rotor model 

An upgraded model for the rotor flow was derived from (Guha and Sengupta, 2013). 
The main difference compared to the previous model, was the assumption of viscous 
flow in place of the equivalent body forces along the radial and tangential directions. 
Adopting this new approach, the Navier-Stokes equations in cylindrical coordinates 
were simplified again; in this case, the body forces were assumed negligible, while 
the viscous terms were still present. Thus, the Navier-Stokes equations were reduced 
to: 

Continuity equation: 


lO0(prw,) | à 

ras (3.49) 
r-Momentum equation: 
Ow, 3 wi 1idp д?м„ 

Wr dr -Q il Um ap ws (3.50) 
0—Momentum equation: 
дуу w,w 0?w. 

wiz = “+ 20w, = va (3.51) 


70 


Lorenzo Talluri 


z-Momentum equation: 


an 0 (3.52) 


The present model introduces an axial velocity profile, so that the relative veloc- 
ities in r and 0 directions may be expressed as: 


We(r,z) = We,6(R)G(z) (3.53) 
wr(r,z) = Wr,5(R)H(2) (3.54) 
Where: 
pios ae. aaj 
2 r2 
Wo(r,z). _ Wr(r,z) 
GT WD Tw) 


G(z) and H(z) are the variations of tangential and radial velocities respectively 
along z direction within the boundary layers. 

Following the procedure outlined in (Guha and Sengupta, 2013), it was initially as- 
sumed that the velocity profile of the fully developed flow was laminar, thus para- 
bolic. Accordingly, G(z) and H(z) could be expressed as: 


Z Z 
G (2) = HG) = 6c (1-7) (3.55) 
and: 
w,(t,2) = W,- ez E z) (3.56) 
Wa(r,Z) = Wo: ez = z) (3.57) 


Integrating the differential form of the 0-momentum and r-momentum equations 
between z=0 and z=b/2, and applying the boundary conditions reported іп (Guha and 
Sengupta, 2013), which assumed maximum velocity value at mid channel and zero ve- 
locity at the walls, it was possible to calculate the gradient of relative tangential ve- 
locity and static pressure in radial direction. 


дм 5 10v 1 
1dp Ow, 6 we 6 12 
Сао ar ‘TN r+ 20% + px er (3.59) 


These equations were also implemented in EES environment and numerically 
solved by applying a step forward method (second-order centered finite difference). 

In order to generalize the mathematical model of the flow, a coefficient for the 
parabolic velocity profile was defined, still under the assumption of laminar flow con- 
dition. Accordingly, G(z) and H(z) can be expressed as: 
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G (2) = H(2) =a? (1-2) = az - a? (3.60) 


Where the coefficient “а” is set to 6 in (Guha and Sengupta, 2013). 

Following the above calculation steps, the reduced 0 and r momentum equations 
were achieved and implemented into the developed EES calculation code, in the same 
way as the previous case: 


Wa _ n 60v n 1 
dro a Gina? P Da (3.61) 
1dp ðw, a? a wê а? 2a 


vq c ge ТОШ 20wg = + T 30 р (3.62) 

In order to further generalize the mathematical model of the flow, a general law 
for the definition of the velocity profile for turbulent flows was introduced. The re- 
quired constraints were the zero and maximum velocity at the walls and at mid chan- 
nel height, respectively. 

The turbulent flow velocity profile law is suggested in (Munson et al., 2013); ac- 
cordingly, the applied power law distribution was adopted for the definition of G(z) 
and H(z): 


1 
z\" (n+1)-(2n+1) 
ee 
5 
Following the previous calculation process, with the only difference in the inte- 
gration of velocity distribution functions, it was possible to calculate the gradient of 
relative tangential velocity and static pressure in radial direction and to integrate them 
in the EES calculation code with a step forward method: 


G(z)=H(z) = | 1 (3.63) 


дм _ (п + 2): n? 
дг — (n + 1)?(2n + 1) 
у 4(2п + 4)n? RE (3.64) 
w, n?b?(n4 1)2n* 1) r va 
1dp dw, n (0+1)? * (20 + 1)2 n 
—— = зуу «= hd E 20 
p dr Ут эг n+2 4n^ OECD OR 
(n+1)*(2n+1) wg n 
"—— + —: 3.65 
2n? r n+2 69) 
(n+ D'-Qn* D', ^ (n+1):(2n+1) 
4n* T pb? 2n? 


For fully developed turbulent flow, n—7 is the most used in literature. For this 
reason, the law is also called the one—seventh power law velocity profile. 

Furthermore, for laminar flow profiles, reference (Munson et al., 2013) suggests the 
following expression for the functions G(z) and H(z): 
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G(z) = H(z) = 2:| 1- (3.66) 


NI O| м 


Which, after the required passages of integration and derivation, corresponds to 
the generalized parabolic formula when the coefficient “a” is equal to 8. 

The last upgrade of the here proposed turbulent rotor model was the estimation of 
the right power law exponent: as reported in (Munson et al., 2013), for turbulent flows 


(generally Reynolds >10*) it can be determined as a function of the Reynolds number. 
(w:2b) 


Rotor Reynolds number is calculated as: Re = at each discretization step. In 


this way, a logarithmic law for n was implemented: 
n = 0.7823 In(Re) — 2.0013 (3.67) 


When the Reynolds number is lower than 2000, the laminar expression derived 
from (Guha and Sengupta, 2013) was adopted. 

When the Reynolds number is between 2000 and 10°, corresponding to a transi- 
tional flow regime, a parabolic distribution with lower value of coefficient “a” can be 
used rather than a power low equation, with the exponent n obtained from an approx- 
imate exponential law. 

In Eq. 3.60, the coefficient “а” was introduced in the parabolic velocity profile in 
order to generalize the model. This coefficient is suggested to be set to 6 in (Sengupta 
and Guha, 2012, Guha and Sengupta, 2013). A more accurate match with CFD results 
was achieved for the fully developed region with “a” equal to 8, as suggested by (Mun- 
son et al., 2013). 

Fig. 3.8 shows the effects of the “а” coefficient on the velocity profiles: high val- 
ues of “а” determine pronounced parabolic distributions, typical of fully developed 
laminar flows, while low values are related to transitional and not fully developed 
laminar flows, characteristic of the entry region. In order to calculate the length of the 
entry region, Eq. 3.68 was applied (Carey, 2010). 

e 3.68 
© 50 ( u ee 

It was found with a comparison with CFD analysis that in the entry region the 
correct value of the “а” coefficient is around 4, as there, the flow is not fully devel- 
oped. This has been confirmed in all analyses conducted, for all fluids (R404a, R134a, 
R245fa, R1233zd(E)) at various rotational speeds (1500, 3000, 4500) and total inlet 
temperatures (60, 80, 120, 140°C) at a total pressure inlet of 1.16. Nonetheless, in 
order to be certain that the “а” coefficient hold the proposed value, an extensive ex- 
perimental campaign, involving visual measurements of the flow field inside a Tesla 
turbine would be required; however, in literature there are no available researches on 
this topic. 
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Fig. 3. 8 Velocity distribution inside the channel for different values of the “а” coeffiient 
The infinitesimal length of the fluid trajectory can be assessed through (3.69). 
ds = yd@2 + dr? (3.69) 
Where d0 = wọ dt is tangential component and dr = w, dt is the radial component, 


as shown in Fig. 3.9. 

When Y^ ds; is lower than the entrance region length, the coefficient for not fully 
developed flow is used; when УЛ ds; is higher than the entrance region length, the 
laminar flow coefficient is adopted. 


Fig. 3. 9 Infinitesimal trajectory 
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Disk profiling 
Another feature of the developed in-house code is the possibility of simulate disk 


profiling. Disk profiling refers to the potentiality of sharpen the outer section of the 
disk, in order to reduce the abrupt contraction losses. 

An exponential law (3.70) has been implemented in order to simulate this feature. 
Another possibility in order to simulate disk profiling is to define a priori the profile 
of the disk and store it in lookup tables, which will be used as input in the procedure. 

ex 
bi = bmin + (bmax — Баһ) * (=e) (3.70) 

The percentage of the disk, which will be sharpened, can be imposed in the code. 
After the sharpening of the edge, the disk maintains a constant profile. 

Rotor Streamlines 

In order to determine the trajectory of the flow inside the rotor, the velocity com- 
ponents in differential form must be considered: 


r do = vgdt (3.71) 
dr = v,dt (3.72) 
The combination of Eqns. (3.71) and (3.72) implies that: 
аө Vg 
(E = (3.73) 
dr Жї 
For the relative trajectory, following the same procedure, it сап be found that: 
d 
(=) = (3.74) 
dr WI 


The determination of the absolute and relative streamlines is of paramount im- 
portance, not only because it allows the visualization of the flow trajectory, but also 
because it is a way to assess the “filling” of the rotor channel (absolute streamlines in 
Fig. 3.10a) and to verify the eventual flow reversal (relative path lines in Fig. 3.10b). 
Specifically, the reverse flow condition occurs when the absolute tangential velocity 
is lower than the rotor peripheral velocity, resulting into a negative relative tangential 
velocity, as it is noticed on the red relative streamline in Fig. 3.10b. 


C---4500rpm \ 


Ж "d 


270° 


Fig. 3. 10 Absolute and relative rotor streamlines with R404a working fluid: (a) Absolute at 
1500 rpm, 4 nozzles; (b) Relative streamlines at various rpm 
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Non adiabatic model 

Initially, Eq. (3.48) (rothalpy conservation) was applied in order to calculate static 
enthalpy at each discretization step. Indeed, in steady state operation, Tesla turbine 
can be considered as an adiabatic expander; however, for transient, and especially at 
start up condition, the turbine rotor cannot be considered as adiabatic. Therefore, Eq. 
(3.48) has been modified to Eq. (3.75). 


-1 42 3.75 
Бе (3.75) 

The heat transfer model in order to determine the heat exchanged between the 
disks and the fluid assumes a constant temperature of the disks. This assumption is 
reasonable, if one takes into account that the velocity at which the heat conduction 
takes place is extremely faster compared to the convection “velocity”. 

In order to determine the convection heat transfer coefficient, the correlations pre- 
sent in (Cengel, 2006) for the internal laminar flows of rectangular shaped pipes were 
applied. Furthermore, the model, as for the viscous case, takes into account the ther- 
mal entry region. In this region, the Nusselt number, conversely to the fully developed 
laminar flow, is not constant, but it is a function of Reynolds and Prandtl numbers. 
The complete set of equation is therefore resumed in Eqns. (3.76) — (3.80). 


Nu = 2.98 for fully developed laminar flow (3.76) 

(0.065+кегрг2') (3.77) 
Nu = 2.98 + —————— for entry region 

10.04: (керг)? 
h; = = (3.78) 
h 

dA = n(r2., – г?) (3.79) 
dQ = hj * dA * (T; — Taisk) (3.80) 


For steady state flow, an iterative procedure has been implemented, which impose 
the temperature of the wall as the area weighted mean temperature of the fluid through 


the rotor. 
T _ ХТА (3.81) 
mean T YN dA 


1.5 Diffuser 


At the rotor outlet, the fluid trajectory is curved and the velocity mainly assumes 
an axial direction, while a considerable amount of the kinetic energy is lost. The pres- 
ence of a diffuser can partially recover this energy through a gradual section enlarge- 
ment, which is able to decrease the fluid velocity, thus converting a fraction of the 
kinetic energy into pressure. The typical solution is the installation of a conical dif- 
fuser at the turbine outlet. The design concerns the choice of the outlet diameter, the 
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length and the diffusion angle; the latter is particularly important to reduce the pres- 
sure losses related to the wall separation of the fluid. 

The model provides the calculation of the axial velocity at diffuser inlet through 
the continuity equation, while the tangential and radial components are conserved 
from the rotor to the diffuser output. Here, the continuity equation ensures the reduc- 
tion of axial velocity due to the increased cross-section, while the radial and tangential 
components were calculated by assuming the conservation of their angle with the axial 
velocity. These velocity components are responsible for a 3D swirled flow. 

Inside the diffuser, the total enthalpy is conserved, while the total pressure is re- 
duced due to the wall friction inside the diffuser: 


1 
Pos = Poa — aip PaVi (3.82) 


This equation is similar to the one used for the abrupt enlargement, but the loss 
coefficient бант was obtained through a polynomial fitting of experimental data avail- 
able in (Idel’chik, 1960) as a function of the diffuser in/out area ratio and the diffusion 
angle. 

А, 
баш = (1— A ni (3.83) 
out 

A further step was the assumption of swirled flow: in (McDonald, and Fox, 1971), 
a large number of diffusers were analysed and tested to evaluate the influence of vor- 
ticity on the performance, concluding that swirled flow could increase the diffuser 
efficiency. Therefore, an incremental coefficient, which depends on the inlet/outlet 
area ratio of the diffuser, on the swirl ratio (i.e. the ratio between tangential and axial 
velocity), and on the diffuser angle, can be applied to the efficiency, defined as: 


mre Ps — P4 © 
diff ^ “7 swirl 
5 pave (3.84) 


The pressure recovery and the outlet static pressure are increased due to the 
swirled flow. 

Therefore, when the velocity, total pressure and total enthalpy are known, it is 
possible to completely define the fluid state at inlet and outlet of the diffuser. 


1.6 Performance indicators 


In order to assess the performance potential of a Tesla turbine, the right design 
parameters need to be defined. Non—dimensional parameters were selected, following 
common practice in turbomachinery (Whitfield and Baines, 1990, Dixon, 2005). There- 
fore, the flow and load coefficient can be expressed as: 


alge (3.85) 


_ work = Vo2U2 r Vo3U3 


The specific speed and the non-dimensional specific diameter are given by: 
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on 
N, = rpm * p2 078 
v2 : (3.87) 
(hz — з) + 2. 
> 0.25 
(hz — з) + + 
8 
3.88 
D, = D; 0.5 ( ) 


The total to static efficiency of the turbine is defined as: 


work №205 — Мз Оз 
Пав з) eee 
Moreover, critical design parameters for the Tesla turbine were identified in the 
geometrical ratios (D3/D2) and (b/D2); for output conditions the exit kinetic energy 
and the absolute flow angle, which should be as low as possible, were identified as 
critical performance indicators. The exit kinetic energy is presented in non-dimen- 
sional form as: 


V3 
m Ekin,3 _ 2 (3.90) 
3 oe 2 К 
„л 


And the exit fluid angle can be calculated as: 


_1 [Уз 
a = tano! (=) 
"a (3.91) 


2. Mechanical Design 


2.1 Static analysis 


Plenum chamber 

The thickness of the plenum chamber is calculated assuming it as a cylinder and 
applying therefore the theory of thin pressure vessels, assuming that the thickness is 
significantly less than the diameter. Eq. (3.92) allows obtaining the minimum thick- 
ness that the plenum chamber has to have; provided that the inner radius, the operating 
pressure and the permissible load are fixed. 


= Pamm ‘Tint 
thin = —— — — (3.92) 


Oamm 
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Rotor disks 

In order to verify the mechanical resistance of the rotor disks, rotating disks theory 
was employed, considering that centrifugal, thermal and pressure forces act on them. 
These forces act both radially and tangentially; however, since the maximum stresses 
develop within the internal radius, the former are zero, while the latter are maximum. 
In addition to the geometrical and thermo-fluid dynamic parameters of the machine, 
itis necessary to know the main characteristics of the material, in this case aluminium; 
these are resumed in Tab. 3.1. 


Table 3. 1 Aluminium properties (at 20?C) 


Density Poisson Young module Thermal expansion coef- Yield tension 
[kg/m?] ratio [GPa] ficient [°C] [MPa] 
2700 0.33 70 0.0000234 180 


The equations provided by rotating disk theory for the calculation of the tangential 
forces are reported. For the stresses due to inertial effects, or centrifugal force, Eq. 
(3.93) is applied (Bhandari, 2014). 


34v rir 1+3v 
c= pw? (r2 + r2 - +2 — —— r?) (3.93) 


For stresses related to the thermal gradient along the disk Eq. (3.94) is employed 
(Bhandari, 2014). 


BET Та 5 rirs 
05-3 2 (2 Hg Ferg + ——— 2r (r2 +гз) (3.94) 
3-I5 rs r 


For stresses linked to boundary conditions, i.e. pressure inside and outside the 
disk, Eq. (3.95) is used (Bhandari, 2014). 
p2r? 


2. .2 
T2? — T3 


r$ 

o=- (5+1 (3.95) 

These formulas do not take into account a safety coefficient, so values of temper- 
ature and pressure higher than the possible achievable conditions during tests have 
been considered, both for the air (10 bar and 150°C) and the ORC (25 bar and 200°C) 
prototypes. 
Power shaft 

The power transmission shaft has to be dimensioned in order to withstand the 
torsional forces; the shaft is therefore schematised as a solid cylinder subjected to 
torque. The torque defined by the ratio between power and rotation speed, as shown 


in Eq. (3.96). 
Power 
T= T (3.96) 
The shear stress acting on the shaft is given by Eq. (3.97) (Bhandari, 2014). 
16: T 
USES (3.97) 
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This value (further multiplied by a safety coefficient) has to be lower than the 
maximum allowable tangential stress, which, according to the Tresca criterion, is ob- 
tained from the yield strength as reported in Eq. (3.98) (Bhandari, 2014). 


1 
Tsn = 2 Osn (3.98) 


Shaft screws 

A further mechanical sizing concerns the screws connecting the shaft to the outer 
plate, which in turn is connected by threaded cylinders to the disks; these screws are 
passing through the disk and are captive on the shaft, therefore they are loaded by 
compression as most of the screws used for these purposes. However, unexpected 
shaft tightening could lead to the severing of the screws; to avoid this risk, their sizing 
is done following shear theory and particularly the theory of bolted unions. 

The tangential stress associated with the torque is calculated through Eq. (3.99). 

Fp 

T7 Dip ip Ap (3.99) 

Where F; is the exchanged force, n, the number of bolts, nj; the number of cutting 
planes and A; the resistant area of the screw. 

This tension (again multiplied by a safety coefficient) has to be lower than the 
tangential stress leading to yield, calculated through Tresca criterion. 
Plenum chamber screws 

The external case was sized to withstand the pressure inside it, but the number 
and size of the screws needed to keep it close still have to be defined. The screws can 
be dimensioned for static failure or for separation, and this last condition being more 
stringent than the first one is preferred for the selection of the screws. 

The dimensioning of a threaded connection is obtained through Eq. (3.100). 


АРТЕ (3.100) 


The left term indicates the overall section of the bolts; while in the right term there 
are the force (Е) acting on the screws, the admissible load (op) and the coefficients K 
and C. The force is calculated by multiplying the internal pressure by the section of 
the box, while the admissible load is known from the type by the class of screws used 
(usually 8.8) and by the safety coefficient (1.25). The standards (Bhandari, 2014) sug- 
gest to assume a coefficient K=0.75, while the coefficient C is calculated through Eq. 
(3.101). 

Ep 

EZ 8 E, +E, (3.101) 


Where E»—E,. represents Young's modulus of steel, so C=1/9. 


C 


2.2 Dynamic analysis 
As thermodynamic conditions, as well as the maximum rotational speed of the air 


Tesla turbine during tests condition were not planned to reach critical conditions, dy- 
namic analysis has been carried out only for the ORC Tesla prototype. 
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Simplified Rotor Geometry 

In order to reduce the computational cost of the simulations, a simplified rotor 
geometry has been considered. Particularly, the simplified geometry has been built 
with the same dimension of the disks (outer diameter of 0.216 m, inner diameter of 
0.055 m) and shafts of the real prototype, but with a reduced number of disks (from 
60 to 11), as shown in Fig. 3.11. The real geometry of the prototype will be throughout 
discussed in Section 3.3.2. 


a) b) 


Fig. 3. 11 Simplified Tesla rotor geometry 


Calculation model in SolidWorks Software — Ideal Bearings Hypothesis 

The calculation of the simplified rotor geometry was performed using SolidWorks 
Simulation frequency analysis tool. The main mesh and solver parameters are shown 
in Tab. 3.2 and Fig. 3.12. Several meshes were created, with the objective of investi- 
gating the effects of elements size on the obtained results. In order to attain a proper 
compromise between accuracy and computational time, the 228,738 nodes mesh was 
finally selected. 

A material selection process for the prototype was carried out. Stresses on rotating 
disks were evaluated through the application of Eqns. (3.93) — (3.95). An aluminium 
alloy (7075) was finally selected. This was mainly due to the lower weight compared 
to iron cast alloys. 


Table 3. 2 Settings used for Software calculations 


Direct Sparse Solver + Soft Spring 
stabilization 
Solid Mesh — Parabolic tetrahedral 
solid elements 
Number of elements 115663 
Number of Nodes 228738 


Solver Settings 


Mesh parameters 
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Mass Properties 


Model Mass [kg] 0.77 
Moment of inertia Ixx [kg*mm?] 2134 
Moment of inertia Туу [kg*mm?] 2134 
Moment of inertia Izz [kg*mm?] 4035 


Fig. 3. 12 Mesh distribution of simulated simplified rotor geometry 


Fig. 3.13 displays the resulting total amplitude for the first mode. Particularly, 
Fig. 3.13 а—с is obtained for 3 different velocities of the expander (3000, 4500 and 
6000 rpm). As it can be noted, if the first mode frequency is reached, the total defor- 
mation of one of the central disposed disk (different for each rotational speed) would 
be too high and therefore it would break. The positive trait is that the resonance fre- 
quency rises with rotational speed; therefore, for the assumed test conditions the res- 
onance frequency is not reached. (The issues arise for rotor velocities values between 
12000 and 15000 rpm.) 


a) b) 
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$ 


Fig. 3. 13 Total Amplitude for ideal bearing stiffness; a) 3000 rpm; b) 4500 rpm; c) 6000 
rpm; d) twisted result of first mode at 6000 rpm 


Ball bearing stiffness model 
The radial deep groove ball bearing stiffness was determined through the appli- 


cation of the model presented in (Dyrobes, 2018), which was developed in (Gargiulo, 
1980). 


К. = 0.0325E06 - /D- F- Z2 · (cosa)? (3.102) 


Where: 

e К. is the Radial Stiffness in [N/m]; 
e Fis the external radial force in [N]; 
e Dis the ball diameter in [m]; 

e 715 the number of rolling elements; 
e cis the contact angle in [rad]. 


On the other hand, the double row angular ball bearing stiffness was determined 
on the basis of the model developed by Guay P. and Frikham A. (Gargiulo, 1980). 

The model is based on the calculation of the ball stiffness, which is built on Hertz 
theory (Hertz, 1882). Therefore, the developed model starts from the preliminary cal- 
culation on curvature to finally obtain the ball stiffness. 

Particularly, the curvature parameters required are shown in Tab. 3.3. 

Applying Hamrock and Anderson model (Hamrock, 1973, Hamrock, 1983) for a sim- 
plified solution of the ellipse elongation, it is possible to calculate all the parameters 
required for the ball stiffness calculations. Eqns (3.103) to (3.105) are used for calcu- 
lating ellipse elongation, elliptic integral of the first and second kind, respectively. 


x = 1.18 - p°598 — 0.19 (3.103) 
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F(k) = + q: (1 + In(p) ) (3.104) 
q 

$(к) 214- 

(к) ie (3.105) 


Table 3. 3 Curvature parameters 


Symbol/Equation 


Parameter 
Ball Diameter D 
Raceway groove curvature radius r 
r 
Raceway conformity f= D 
Contact angle a 
Pitch Diameter din 
: а D: cosa 
Dimensionless parameter y {дел 
m 
Equivalent curvature radius 


Inner raceway Outer raceway 


D 
Rye = (1+) 5 


D 
Ray 
D f.-D 
ИЕ ala Ri = ——— 
(2—1) 21) 


Finally, the stiffness of ball/raceway contact, as well as the ball stiffness, can be 


found respectively from Eq. (3.106) and Eq. (3.107). 


т 25г (3.106) 
кекек 
3 
“2 
1 1 3.107 
K, =| 5+5 (3.107) 


Where e and i suffixes are used for external and internal raceways. 
Axial and radial stiffness can thus be computed applying Eq. (3.108) and Eq. 


(3.109), with the assumption of constant contact angle. The assumption of constant 
contact angle is justified for contact angles higher than 25° (Gargiulo, 1980). 


_3 LR кеу? (3.108) 
ka = 5° 2° Ky: (sina)? (5) 
3 5 1 
k; = ——c':Z:K,: (cos а)2 · (е · tana)2 (3.109) 
к= n ` (cos a)2 + ( ) 
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For paired bearing the axial stiffness is worth twice the single bearing axial stiff- 
ness, therefore Eq. (3.108) should be multiplied by 2. 

Tab. 3.4 resumes the main characteristics, while Tabs. 3.5 and 3.6 display the 
values of radial and axial stiffness for the Tesla prototype bearings. In Tab. 3.4, the 
preload of the bearing has been calculated as suggested in (SKF, 2018a); that is taking 
the minimum load on ball bearing as 0.01C, where C is the basic dynamic load rating. 
In Tab. 3.6, very low preload value, of 50 N, has been considered in order to evaluate 
the effect of preload on modes frequencies. 


Table 3. 4 Bearings main parameters (SKF, 2018b, SKF, 2018c) 


Deep groove ball bearing SKF 61817-2RZ 

Pitch diameter dm 97.5 

Ball diameter D 7.14 
Number of rolling elements Z 24 

Double row ee contact ball bear- SKF 3207 A-2ZTN9/MT33 

Pitch diameter dm [mm] 53.5 

Ball diameter D [mm] LETI 
Contact angle [?] 30 
Number of rolling elements Z 18 


Table 3. 5 Bearings stiffness for ORC Tesla prototype; preload = 0.01C 


Deep groove ball bearing SKF 61817-2RZ 
Axial Stiffness К. [N/m] - 
Radial Stiffness k, [N/m] 2.36E+08 

Double row icr contact ball bear- SKF 3207 A-2ZTN9/MT33 
Axial Stiffness К. [N/m] 3.62E+08 
Radial Stiffness k, [N/m] 5.44E+08 


Table 3. 6 Bearings stiffness for ORC Tesla prototype; preload = 50 N 


Deep groove ball bearing SKF 61817-2RZ 
Axial Stiffness k, [N/m] — 
Radial Stiffness k, [N/m] 6.95E+07 
Double row angular contact ball bear- 
m SKF 3207 A-2ZTN9/MT33 
Axial Stiffness k, [N/m] 8.41E+07 
Radial Stiffness k, [N/m] 1.26E+08 
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Calculation model in SolidWorks Software —Bearings with Real Stiffness 

Same settings as the ones used for ideal bearing stiffness calculations were used 
also for the case with real bearing stiffness. Figs. 3.14 and 3.15 display the resulting 
total amplitude for the first mode, when real bearing stiffness is applied (Fig. 3.14 
with Р = 0.01 C and Fig. 3.15 with Р = 50 N). Particularly Figs. 3.14 and 3.15 are the 
analogue of Fig. 3.13. As can be noted from the comparison of Figs. 3.13, 3.14 and 
3.15, the effect of real bearing stiffness model slightly influences the resonance fre- 
quency as well as the shape of deformations. The only significant change is that more 
disks display high deformation, even if with a lower amplitude. The preload value 
slightly affects the modes frequency values. 


a) b) | 
d) | 


Fig. 3. 14 Total Amplitude for real bearing stiffness; a) 3000 rpm; b) 4500 rpm; c) 6000 rpm; 
d) twisted result of first mode at 6000 rpm 
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Fig. 3. 15 Total Amplitude for real bearing with low stiffness; a) 3000 rpm; b) 4500 rpm; c) 
6000 rpm 


Fig. 3.16 displays the comparison between ideal and real bearing stiffness modes 
frequencies. For low rotational speeds the frequencies obtained are very close, for 
higher rotor velocities, the influence of bearing stiffness is more marked, as it slightly 
lower the modes frequency values. 
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Fig. 3. 16 Ideal vs. real bearing stiffness modes frequencies 


The modal analysis of the rotor of the ORC Tesla turbine prototype, through the 
meand of a commercial solver (Solidworks simulation) was performed; first assuming 
perfect bearing stiffness, after, a simplified model to assess the stiffness of the selected 
bearings has been developed. The modal analysis has been carried out for different 
rotational speeds, both for ideal and real stiffness values. Two different stiffness val- 
ues have been obtained, one with a preload value recommended by the bearings man- 
ufacturer and another with a very low preload. The influence of preload in the modal 
analysis is negligible. Frequency of the modes calculated with real bearings stiffness 
are very close to the ones calculated with ideal bearing stiffness. To conclude, for test 
case conditions (velocities lower than 9000 rpm), the design prototype seems not to 
be experiencing resonance. 


3. Prototypes design: from thermodynamic considerations to realization 


In order to reach the final prototypes geometry, various factors were taken into 
account. First thermodynamic optimization was carried out, but without forgetting to 
perform a mechanical analysis of the optimal thermodynamic configuration, choosing 
therefore the right compromise and also keeping in mind manufacturing process and 
test benches limitations. 
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3.1 Air Tesla turbine 


Design constraints and limitations 

The air Tesla turbine prototype parts were built through different processes. Par- 
ticularly, the stator was made with a 3D printing technique in ABS material. There- 
fore, one of the boundary conditions of the design was the printable area of the 3D 
printer. 
3D printer and ABS material 

3D printing is a technique based on the overlapping of layers of plastic material 
heated to relatively high temperature. At the proper temperature, the material becomes 
malleable, and therefore it is deposited on a plate, layer after layer, in order to produce 
the desired piece. The contact with air provides the cooling and therefore the stiffening 
of the material, returning it to solid state. 

The printer used to make the stator of the air Tesla prototype is the Anet A8 model 
manufactured by Prusajr (Fig. 3.17). Tab. 3.7 resumes its main specifications. 


Table 3. 7 3D Prusajr Anet A8 characteristics 


5: Max: noz: Max plate Nozzle Printing Possible 
Printing Area zle Tem- : + RACE 
tempera- diameter velocity printing 
[mm] perature б А 
[°С] ture [°C] [mm] [mm/s] materials 
PLA, 
220x220x240 250 100 0.4 100 ABS, TPU, 
PVA, PP 


The main boundary condition given by the 3D printer was the maximum printing 
area. It was found while using the 3D printer that, in order to have a good quality of 
the printed stator, the maximum external stator diameter had to be of about 150 mm. 

The material used for printing the stator is ABS, an acronym for acrylonitrile— 
butadiene-styrene, a thermoplastic polymer frequently used in industry and, thanks to 
its thermal and mechanical properties, particularly suitable for three—dimensional 
moulding. ABS in fact has a melting point of 100°C, which allows the resistance of 
the printer to work at a temperature between 215°C and 250°C, and the material to 
cool rather quickly once deposited. ABS has good impact resistance, high toughness 
and rigidity. Thermal resistance is good a low temperature, however, for power pro- 
duction it determines a limitation, as it fluid inlet temperature affects the machine 
performance. Manufacturers recommend using ABS at temperatures below 85°C, 
which prevents the air from being heated to higher temperatures. The main properties 
of ABS are resumed in Tab. 3.8. 


Table 3. 8 Main ABS properties (at 20 °C) 


Density Poisson ratio Young module Yield strength Bs 
3 
[kg/m°] [GPa] [MPa] [MPa] 
1060 0.35 2500 45 70 
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Friction co- Melting tem- Temperature a expan- Relative per- 
efficient on rature [°C] range of opera- sion coeffi- mittivi 
steel регаш tion cient[°C-] ty 
0.5 85/130 -40 / 80 0.00005 3:5 


Fig. 3. 17 3D Prusajr Anet A8 


The other parts of the air Tesla turbine prototype were built through the different 
processes. The outlet casing is made by Plexiglas material and it has been made 
through a milling process. The rotor disks are made of aluminium through laser cut- 
ting and the shaft are made of steel (through turning). The reason for choosing metal 
materials for the rotor lies in the greater stresses and deformations to which is sub- 
jected; in particular, the disks, due to their limited thickness and the mechanical and 
thermal loads (even if limited) to which they are subjected, could run the risk of un- 
dergoing high deformations (verified through the application of the model discussed 
in Section 3.2.1), altering the shape of the channels. Nonetheless, the choice of these 
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materials allowed not having further restrictions given by the materials or the manu- 
facturing process. 
Test bench 

The air Tesla turbine prototype was designed to a specific test bench, which com- 
prised a centrifugal compressor (MICO 90/8 model produced by SAVIO SRL), whose 
speed can be regulated with a variable-frequency inverter. Fig. 3.18 shows the per- 
formance curves of the centrifugal compressor. As can be noted, the reachable pres- 
sure is not high (max 1.85 bar), therefore while designing the turbine, this had to be 
taken into account. 


Condizioni al'aspramone `- Sueton conditions 
Conditions à 'aspination : - Bedingungen zu der Zugivafi у 90/8 
Poar= 760 [аан] Temp=0f'C] Dens. = 1292 fg] COMPRESSIONE 
Compression - Compression + Verdichtung 
KPa 
Г [— 1. 
—— = “SL 
ia ~ M 


Fig. 3. 18 Centrifugal compressor performance curves 


The turbine furthermore needs to be connected to an electric motor (brushless 
type) that controls the rotational speed through a servo drive. By means of a continu- 
ous exchange of electricity between the electric motor and the grid, the brushless is 
able to operate as generator or as brake, according to the power produced by the tur- 
bine. The maximum allowable rotational speed by the generator is 3000 rpm. The 
torque and, consequently, the power produced by the turbine is measured by a torque 
meter (Lorenz, nominal torque 10 Nm) connected between the turbine and the brush- 
less motor through flexible couplings. 

The limitation on rotational speed is an important feature to take into account, as 
the Tesla turbine prototype, and generally Tesla turbines, as it will be discussed deeply 
in the results Section, perform well at relatively high rotational speeds. 

Design objectives 

Taking into account the boundary conditions provided by materials, manufactur- 
ing processes and test bench, some design objectives where selected in order to have 
significant tests results. Particularly, a 1096 thermodynamic efficiency and a 100 W 
power output have been taken as reachable objectives. Having therefore stated the 
limitations and objectives, the proper dimensioning of the turbine can be carried out. 
Plenum chamber 

The turbine entry is made by 4 different “Rilsan pipes" in order to have a homo- 
geneous distribution inside the plenum chamber. The toroidal plenum chamber allows 
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having low velocities at stator inlet, which, as demonstrated in (Guha and Smiley, 2009, 
Neckel and Godinho, 2015), are beneficial to the turbine efficiency. 

The mass balance equation states that for a fixed mass flow rate and thermody- 
namic conditions, in order to have low velocities, a wide area is required. The first 
thing to dimension is therefore the required section area of the plenum chamber, which 
will allow the flow to slow down, passing from 4 circular duct of 3/4" to a toroidal 
chamber of greater section. As the lower boundary of the plenum chamber is given by 
the external diameter of the stator, the first geometry parameter that has to be defined 
is the height of the plenum chamber (which will determine the internal radius of the 
carter). Assuming the maximum achievable mass flow rate of 60 g/s at 1.8 bar and 
60°C, the velocity at entry of the plenum chamber is 28 m/s. In order to have velocities 
in the plenum chamber lower than 1.5 m/s the height of the plenum chamber should 
be around 35 mm. Therefore, assuming maximum external stator diameter achievable 
with the 3D printer, an internal carter diameter of 0.22 m is required. 

The thickness of the external carter is sized by the theory of thin pressure vessels; 
assuming that its value is significantly lower than the diameter. Eq. (3.92) allows ob- 
taining the minimum thickness that the component must have. From Eq. (3.92), the 
minimum required thickness for the test conditions would be of 3 mm. Nonetheless, 
it was decided to design the turbine in order to withstand 10 bar of internal pressure, 
in case future tests in other test benches could be run. Therefore, the minimum re- 
quired thickness obtained from Eq. (3.92) is of 13 mm and taking into account the 
diameter of the screws, the final minimum radial thickness of the plenum chamber 
was set at 20 mm. 


Table 3. 9 Final external carter dimensions 


Carter external dimen- Carter internal Carter radial mini- Height of plenum 
sions [m] radius [m] mum thickness [m] chamber [m] 
0.26x0.26x0.1 0.11 0.02 0.035 


Fig. 3.19 displays the Solidworks drawing and the realized carter of the prototype 
in PMMA material, while Tab. 3.9 resumes the carter dimensions. 


Fig. 3. 19 External carter of the turbine, SolidWorks drawing and real prototype 
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Stator 

The design of the stator, due to the external diameter limitations has been quite 
straightforward. Indeed, applying the inlet/outlet diameter ratio recommended by 
(Glassman, 1976, Ventura et al., 2012), ithas been possible to select the internal diameter, 
which allows highest efficiencies. 

Knowing internal and external radius of the stator, the last and most important 
things to determine are number of nozzles, the height and width dimensions, as well 
as the exit angle. The nozzle exit angle has been taken as the maximum that the man- 
ufacturing process could achieve, or 85°. Indeed, in order to obtain maximum effi- 
ciencies, a completely tangential fluid to the rotor is required. The 3D printing limited 
the throat width at 0.8 mm, as the uncertainty of the printing process is of 0.2 mm. 

Therefore, the last dimensions to determine have been stator height and number 
of nozzles. Tesla turbines efficiency is high when high Mach number and low mass 
flow rates are reached. Particularly, increasing height of the nozzle, as well as the 
number of nozzles, would bring to a penalty in efficiency. Nonetheless, low number 
of channels would not guarantee and homogeneous flow; therefore 4 nozzles per stator 
have been selected. Nozzle number optimization is being further assessed in Section 
4.1.1.2. 

The stator height has been defined so that each stator would provide the flow to 
ten channels; therefore, the height of the stator is obtained from the sum of 9 disks 
thickness plus 10 rotor channels width. Tab. 3.10 resumes the main dimensions of the 
nozzle and Fig. 3.20 displays the Solidworks drawing and the realized stators of the 
prototype in ABS material. 


Table 3. 10 Final stator dimensions 


Stator external di- Stator internal di- Number of noz- 
Chord length [m] 
ameter [m] ameter [m] zles 
0.152 0.126 0.018 4 
Stator inlet angle [°] Stator outlet angle Nozzle throat Nozzle throat 
[°] width [m] height [m] 
0 85 0.0008 0.012 


VA 


Fig. 3. 20 Turbine stator, SolidWorks drawing and real prototype 
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Rotor 

Larger is rotor diameter, higher is the power produced par channel, therefore, the 
maximum possible rotor diameter, taking into account a 1 mm gap between stator and 
rotor was chosen. 

Once the diameter is chosen, in order to optimize the rotor, two main parameters 
have to be selected, the rotor outlet/inlet diameter ratio and the width of the channels; 
furthermore, the total number of disks has to be defined as well. 

The most significant parameter of a Tesla turbine is the width of the rotor chan- 
nels, which deeply influences the performance of the machine. Particularly, tight gaps 
are required for high efficiency and performance. Therefore, an analysis on the per- 
formance of the machine varying rotor gap at various inlet conditions was carried out. 

Fig. 3.21 display the variations of total to static efficiency (a—c—e) and power pro- 
duction par channel (b-d-£), for three different mass flow rates and four different total 
inlet pressures. At low pressure and mass flow rate optimal channel width is a bit 
larger (in the range of 0.4—0.5 mm), but as pressure and mass flow rate increases, a 
tighter gap is preferred (in the order of 0.3 mm). Higher pressure allows higher total 
to static efficiency at fixed geometry and mass flow rate. This is directly linked to the 
velocity decrease at nozzle stator exit. Indeed, higher pressure means higher density 
values, and accordingly to mass balance equation, in order to maintain the same mass 
flow rate at fixed section, velocity has to be decreased. 

Theoretically, Tesla turbines work better at high velocities, but it has to be kept 
in mind that the maximum reachable rotational speed is limited to 3000 rpm, there- 
fore, optimal tangential velocity ratio is not optimized. 

Power produced per channel is directly linked by throat Mach number value. Par- 
ticularly, as the throat section is defined as the sum of the disks thickness plus the 
channels width, increasing the gap between the channels directly augment the throat 
section and therefore lower throat Mach number values are reached. Consequently, in 
order to have a proper power production, tight channels are preferred. 
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Fig. 3. 21 Efficiency (a-c-e) and power par channel production (c—d—f) for a Tesla turbine 
working with air at various pressure and mass flow rates, as function of channels width 


The other fundamental parameter of Tesla turbine rotor is the ratio between rotor 
inlet and outlet radius. Fig. 3.22 displays the efficiency (a-c-e) and the power pro- 
duced per channel (b—d-f), at various mass flow rate and total inlet pressures, as func- 
tion of the rotor outlet/inlet radius ratio. Optimal values of efficiency are reached for 
a rotor radius ratio of about 0.4, for all mass flow rates and pressures. On the other 
hand, the power produced per channel is increasing as the rotor radius ratio decreases, 
as it implied that the exchange area is increased. Therefore, for optimal efficiency a 
rotor radius ratio of 0.4 has to be chosen, while for maximum power a rotor radius 
ratio of 0.1 is preferred. Taking into account these considerations, the selected design 
ratio chosen was of 0.25, for two main reasons: (i) to have high efficiency with rela- 
tively high power production; (ii) to not have bearings at turbine outlet with high di- 
ameter, as low diameter bearings produce lower mechanical losses. 
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Fig. 3. 22 Efficiency (а—с—е) and power per channel production (c—d—f) for a Tesla turbine 
working with air at various pressure and mass flow rates, as function rotor outlet/inlet radius 
ratio 


Finally, the total number of channels had to be chosen. At design condition (60 
g/s at 1.8 bar inlet pressure), the power produced per channel resulted of 5.5 W. There- 
fore, in order to obtain the desired value of 100 W, 20 channels would be enough. 
Nonetheless, taking into account mechanical losses, and keeping in mind the possi- 
bility of further tests at higher pressure, double the number of the minimum channels 
was selected (40). Tab. 3.11 resumes the main dimensions of the rotor and Fig. 3.23 
displays the Solidworks drawing and the realized rotor of the prototype in aluminium 
and steel materials. 


Table 3. 11 Final rotor dimensions 


Rotor inlet diameter [m] Rotor outlet diameter [m] Channel width [m] 
0.125 0.031 0.0003 
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Fig. 3. 23 Turbine rotor, SolidWorks drawing and real prototype 


Full assembly 

Once the three main parts of the machine have been dimensioned, all other com- 
ponents needed to be assessed, from the rotor shafts, to the bearings, and the screws 
of the carter. Each component has been verified according to the formula displayed in 
Section 3.2.1 and the final dimensions of each component have been resumed in Tab. 
3.12. 

Finally, each drawn component has been assembled in a final assembly, as shown 
in Fig. 3.24, in order to verify if some inconsistencies were present. Once the final 
drawing was completed, all components were ordered and assembled in Linea Labor- 
atory of University of Florence. 


Table 3. 12 Final prototype dimensions 


А оа Сар ва: Power shaft Washers dimen- 
Bearings ber of chan- tor/Rotor . i 
diameter [m] sions 
nels [m] 
SKF 1205 
EKTN9 
SKF 61811 — 40 0.001 0.02 0.003x0.006x0.0003 
2RZ 
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Fluid side ; : 
аа. Rötòrscrews Disks-shaft Carter Sheet gasket thick- 
ian Screws Screws ness 
0.055 M3 M5 MIO 0.00015 


Fig. 3. 24 SolidWorks drawing and real prototype 
3.2 ORC Tesla turbine 


As for the air case, in order to proceed to the design of the ORC prototype, the 
boundary conditions have to be firstly analysed. The limitations on the turbine design 
in case of the ORC prototype were given mostly by the test bench and by manufac- 
turing issues when dealing with nozzle outlet width. 

Test bench 

The ORC prototype was designed in order to work in a test bench available in 
University of Florence (Linea Laboratory) working with R404a as working fluid. The 
scheme of the test bench is showed in Fig. 3.25. 
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TEST BENCH CONCEPT Replaced by Mixing between the cold two-phase 
fluid (4) and the superheated vapor from the 
compressor discharge (5) 
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Fig. 3. 25 Scheme of University of Florence Test bench (Linea Laboratory) 
The test bench is composed by the following components: 


e Avolumetric compressor of the H5 series, model H3400CC manufactured by 
Dorin, which has the task of compressing the vapour of the working fluid to 
the required pressure; as known, once the rotational speed is set, the volumet- 
ric compressors are able to compress a constant volumetric flow of gas over 
a wide range of compression ratios. The model used, operating at a frequency 
of 50 Hz, is able to process a flow rate of 102.35 m?/h absorbing a maximum 
electrical power of 35.2 kW; these limits have to be kept in mind when as- 
sessing the mass flow that will evolve in the turbine and the maximum inlet 
pressure to the expander. 

e The expander to be tested, placed in parallel to a bypass valve for hot gases. 
The torque meter used to measure the torque is the same as that used in the 
air test bench, therefore the torque must be less than 10 Nm to avoid damage 
to the instrument; the brushless motor will also have the same limit as seen 
in the analysis of the air machine (3000 rpm). However, in this test a speed 
reducer will be placed in series that will allow the machine to operate up to 
4500 rpm, allowing to obtain a greater production of power and a better effi- 
ciency. Finally, the axial thrust will influence the choice of bearing type and 
size. 

e  Anair-cooled condenser, i.e. a gas-air heat exchanger capable of cooling the 
vapour of the working and bringing it to saturated liquid conditions. The con- 
denser is designed to exchange a maximum amount of heat of 25 kW, while 
the use of air in ambient conditions (approximately 15?C), requires that the 
saturation temperature is at least 5-10°C higher, to which is associated a min- 
imum value of saturation pressure: the diagram 7-s in Fig. 3.26 shows that, 
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at a saturation temperature of 23°C corresponds a pressure of 12 bar, which 
will constitute the minimum value of Poo. 

e A throttling valve, through which the cooled vapour is expanded and brought 
back to the same pressure conditions as the fluid expanded in the turbine. 

e А storage tank, where the flows of the two branches mix and exchange heat. 
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Fig. 3. 26 R404 Temperature—Entropy diagram, highlighting saturation pressure (12 bar) at 
23°C 

The test bench (Fig. 3.27) compresses the entire flow of fluid to a minimum pres- 
sure of 12 bar, after which the flow is separated into two currents: the former expands 
into a turbine, making the required pressure drop, the latter is cooled and condensed 
under saturated liquid conditions, then expanded into a throttling valve to the turbine 
discharge pressure and finally mixed. The hottest fluid leaving the turbine will return 
the entire flow rate to conditions of overheating, ensuring a temperature of about 10°C 
higher than saturation conditions, in order to ensure the absence of small droplets at 
the entrance of the compressor. 

Therefore, once the upper and lower pressure levels are fixed, the efficiencies of 
the two machines and the flow rate of fluid required in the turbine are determined; 
using a code developed in EES environment, it is possible to obtain the total flow rate 
processed by the compressor, the power absorbed and the heat exchanged in the con- 
denser. It is thus possible to evaluate whether the experimental conditions respect the 
limitations imposed by the test bench; the program has also been implemented in such 
a way that, once the diameters of the ducts are known, it is possible to calculate the 
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relative speeds and pressure losses through an appropriate internal EES function 
which determines pressure drops. 


o> 
è 


Receivers Compressor 


Fig. 3. 27 University of Florence Test bench (Linea Laboratory) 


Developing a parametric analysis varying the three fluid dynamic parameters (up- 
per and lower pressure and mass flow rate), it is possible to obtain the operating range 
of the turbine: for each inlet pressure the machine discharge pressure and the evolving 
flow rate are varied, the results are excluded for which the heat exchanged in the con- 
denser exceeds 25 kW and the work carried out by the compressor is greater than 35 
kW. Fig. 3.28, referring to an inlet pressure of 12 bar, shows the trend of the flow rate 
processed by the compressor as a function of the discharge pressure (or equivalent, 
the compressor inlet pressure) and the flow rate used in the turbine; the green curve, 
on the other hand, refers to the maximum flow rate that the compressor is able to 
dispose of, calculated as a product between the volumetric flow rate (V =102.34 m?/h) 
and the density of the inlet refrigerant, a function of pressure. This line defines the 
operating range in which the test bench is able to operate, thus identifying a minimum 
limit to the turbine exhaust pressure; this limit is more stringent the greater the flow 
rate processed by the expander. 
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Fig. 3. 28 Max. and actual compressor flow rate, as function of turbine mass flow rate and 
lower pressure of the test bench for an upper pressure (Poo) of 12 bar 


When higher inlet pressures are considered, such as the case presented in Fig. 
3.29, by operating the turbine with a flow rate of 0.1 kg/s, the fluid could expand up 
to atmospheric pressure, while a flow rate of 1 kg/s could expand only up to 8 bar, 
limiting the pressure drop to 4 bar. At higher inlet pressures, the flow rate of fluid 
elaborated by the compressor does not change for the same value of (mp); however, 
the limits linked to the maximum heat exchanged in the condenser are more stringent, 
so much so that in the graph of Fig. 3.29, referring to an upper pressure (Poo) of 24 
bar, there is no intersection between the operating lines. 
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Fig. 3. 29 Max. and actual compressor flow rate, as function of turbine mass flow rate and 
lower pressure of the test bench for an upper pressure (Poo) of 24 bar 
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Design objectives 

As for the air expander design, the boundary conditions have been taken into ac- 
count, considering materials requirements, manufacturing processes and test bench 
operation range; therefore, some design objectives have been selected in order to ob- 
tain significant tests results. Particularly, a 40% thermodynamic efficiency and a 500 
W power output have been taken as reachable objectives. Having therefore stated the 
limitations and objectives, the proper dimensioning of the turbine can be carried out. 
Components design 

The design of each component was carried out with the same procedure applied 
to the air Tesla turbine prototype. Tesla turbines, as will be further assessed in the 
result Section, produce high power when inlet rotor diameter is greater, while keeping 
high efficiency level if proper thermodynamic conditions are provided. The aim of 
the design was to manufacture an efficient turbine with relatively limited dimensions. 
The maximum possible outer dimension that could not be exceeded by the casing was 
of 0.4 m, imposed by test bench characteristics. The dimensioning of the turbine, 
keeping in mind the maximum allowed outer casing diameter started from the rotor. 
Rotor 

The design of the rotor started from the evaluation of the external diameter. Var- 
ious total inlet pressures, mass flow rates and rotational speeds were assessed, in order 
to find an optimal configuration, which could comply with the test bench and the set 
objectives. As shown in Figs. 3.30 and 3.31, efficiency and power per channel are 
higher at low pressure and high rotor diameter. Furthermore, at a fixed mass flow rate, 
a higher rotational speed allows to achieve higher efficiency at lower rotor inlet di- 
ameter. The power developed par channel is fairly low for the mass flow rate consid- 
ered, while the efficiency values are in the range of the other competing technologies 
for micro power generation. The external rotor diameter has thus been selected trying 
to have a relatively high power production at limited dimension and high efficiency. 
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Fig. 3. 30 Efficiency and power per channel of a Tesla turbine working with R404a, as function 
of total inlet pressure and external rotor diameter, for a rotational speed of 3000 rpm 


The design point conditions have been taken to be a rotational speed of 4500 rpm 
and a mass flow rate of 0.25 kg/s. Considering therefore these conditions, the balanced 
selection of the radius was a value of 0.108 m, which could guarantee an efficiency 
over 40% for a very wide test range, with a power production close to 500 W at 12 
bar of total inlet pressure, when a total of 60 rotor channels is considered. 
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Fig. 3. 31 Efficiency and power per channel ofa Tesla turbine working with R404a, as function 
of total inlet pressure and external rotor diameter, for a rotational speed of 4500 rpm 


After selecting the outer rotor diameter, the other fundamental parameters, as for 
the air Tesla prototype are the width of the rotor channels and the rotor outlet/inlet 
radius ratio. Fig. 3.32 displays the efficiency (a) and the power produced per channel 
(b), at various total inlet pressures, as function of the rotor channel width. Optimal 
values of efficiency and power per channel are reached for the smallest assessed width 
(0.08 mm); nonetheless, after a market assessment, spacers of 0.08 mm were hard to 
find, therefore a rotor channel width of 0.1 mm was selected, as it still guaranteed 
high efficiency and high power par channel. Fig. 3.32 also displays the efficiency (c) 
and the power produced per channel (d), at various total inlet pressures, as function 
of the rotor outlet/inlet radius ratio. Optimal values of efficiency are reached for a 
rotor radius ratio between 0.25 and 0.4, tending to the lower value when total inlet 
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pressure increases. On the other hand, the power produced per channel is steadily 
increasing as the rotor radius ratio and the pressure decrease; which is directly linked 
to the increase in exchange area. Therefore, for optimal efficiency a rotor radius ratio 
of 0.25 has to be chosen, while for maximum power a rotor radius ratio of 0.1 is pre- 
ferred. Taking into account these considerations, the selected design ratio chosen was 
of 0.25, for the same main reasons discussed for the air Tesla turbine prototype: (i) to 
have high efficiency with relatively high power production; (ii) to not have bearings 
at turbine outlet with high diameter, as low diameter bearings produce lower mechan- 
ical losses. 
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Fig. 3. 32 Efficiency and power per channel of a Tesla turbine working with R404a, as function 
of total inlet pressure and (a-b) rotor channel width, (с-а) outlet/inlet rotor radius ratio, for a 
rotational speed of 4500 rpm 


Finally, Tab. 3.13 resumes the main dimensions of the realized rotor, while Fig. 
3.33 displays the Solidworks drawing and the realized rotor of the prototype in alu- 
minium and steel materials. 


Table 3. 13 Final rotor dimensions 


Rotor inlet diameter [m] Rotor outlet diameter [m] Channel width [m] 
0.216 0.055 0.0001 
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Fig. 3. 33 ORC Turbine rotor, SolidWorks drawing and real prototype 


Stator 

As for the design of the air Tesla prototype, the stator geometry is straightfor- 
wardly defined once the geometric and fluid—dynamic boundary conditions are im- 
posed. The design hypotheses are summarised briefly below: 

e  Inlet/outlet diameter ratio = 1.25, as recommended in (Glassman, 1976, Ventura 


et al., 2012); 

e Nozzle exit angle = 85° (maximum achievable due to manufacturing pro- 
cess); 

e Throat section = 1 mm (maximum achievable due to manufacturing process). 


The actual number of channels is chosen maximizing the velocity ratio o, (Section 
3.1.2): in the case of the air prototype, the øn coefficient was constantly decreasing as 
the number of channels increased, reducing the stage efficiency. In this case øn is 
initially increasing until the number of nozzles is equal to 5 and then it decreases in a 
monotonous way, as shown in Fig 3.34. Therefore, in order to optimize stator effi- 
ciency, the five channels configuration would be the optimal one; nonetheless, the 
minimum flow rate is reached using two nozzles. A compromise between stator and 
rotor efficiency, as well as considerations related to uniformity, symmetry and distri- 
bution of the jet, it was decided to adopt a stator with four channels, with nozzles 
placed at an angular distance of 90? and a coverage degree (ratio between the segment 
of input to the rotor multiplied by four and the entire circumference) of 10%. 

Tab. 3.14 resumes the main dimensions of the stator and Fig. 3.35 displays the 
Solidworks drawing and the realized stators of the prototype in aluminium material. 
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Fig. 3. 34 Stator efficiency as function of total number of nozzles 


Table 3. 14 Final stator dimensions 


Stator external di- Stator internal di- Number of noz- 
Chord length [m] 
ameter [m] ameter [m] zles 
0.272 0.217 0.059 4 
Stator inlet angle [°] Stator outlet angle Nozzle throat Nozzle throat 
[°] width [m] height [m] 


0 85 0.001 0.001 


Fig. 3. 35 ORC Turbine stator, SolidWorks drawing and real prototype 
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Plenum chamber 

The turbine entry is made by 2 1-inch pipes in order to obtain a uniform distribu- 
tion inside the plenum chamber. The required section area of the plenum chamber, 
which allows the flow to slow down, passing from 2 circular duct of 1" to a toroidal 
chamber of greater section has to be defined. As for the air Tesla prototype, the lower 
boundary of the plenum chamber is given by the external diameter of the stator, there- 
fore, the height of the plenum chamber (which will determine the internal radius of 
the carter) has to be determined. Assuming the maximum achievable mass flow rate 
of 1 kg/s at 12 bar and 80°C, the velocity at entry of the plenum chamber is 22 m/s. 
In order to have velocities in the plenum chamber lower than 1.5 m/s the height of the 
plenum chamber should be around 19 mm, which is lower than 1". Therefore, the 
selected internal carter diameter has been of 0.33 m, in order to have a small space on 
both sides of the 1" entrances. 

The thickness of the external carter is sized by the theory of thin pressure vessels; 
assuming that its value is significantly lower than the diameter. The minimum thick- 
ness is obtained from Eq. (3.92) taking into account maximum achievable test condi- 
tions (25 bar), and it was found to be of 10 mm. Considering also the diameter of the 
screws and a safety coefficient, the final minimum radial thickness of the plenum 
chamber was set to 20 mm. Fig. 3.36 displays the Solidworks drawing and the realized 
carter of the prototype in aluminium material, while Tab. 3.15 resumes the carter di- 
mensions. 


Table 3. 15 Final external carter dimensions 


Carter external Carter internal Carter radial mini- Height of plenum 
dimensions [m] diameter [m] mum thickness [m] chamber [m] 
Ф0.37 x 0.07 0.33 0.02 0.03 


Fig. 3. 36 ORC External carter, SolidWorks drawing and real prototype 
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Full assembly 

As for the air Tesla prototype, once rotor, stator and plenum chamber dimensions 
were defined, all other components needed to be assessed, such the rotor shafts, the 
bearings and the screws of the carter. 

A particular feature, which differentiates the air and ORC expanders, is the mag- 
netic coupling of the latter one. Indeed, in order to have a completely sealed turbine 
configuration, the power shaft has been designed in order to convey power through a 
magnetic coupling. 

Each designed component has successively been verified according to the formula 
displayed in Section 3.2.1 and the final dimensions of each part have been resumed in 
Tab. 3.16. 

Finally, each drawn component has been put together in a final assembly, as 
shown in Fig 3.37, in order to verify if some inconsistencies were present. Once the 
final drawing was completed, all components were ordered and assembled in Linea 
Laboratory of University of Florence. 


Table 3. 16 Final ORC prototype dimensions 


Total num- Gap Sta- Power shaft Washers dimen- 
Bearings ber of chan- tor/Rotor diameter 3 
sions 
nels [m] [m] 
SKF 61817-2RZ 
SKF 3207 А— 60 0.001 0.035 0.004x0.008x0.0001 
2ZTN9/MT33 
Fluid side shaft Rotor Disks-shaft Carter З А 
a O-ring section [m] 
diameter screws screws screws 
0.085 M4 М5 M8 0.00262 
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Fig. 3. 37 SolidWorks drawing and real ORC Tesla prototype 
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4. 3D CFD Model 


The fluid behaviour inside the Tesla turbine rotor of both prototypes was analysed 
through the means of computational fluid dynamics; while only for the ORC config- 
uration stator—rotor coupled simulations were carried out. 

The study of the real gas flow through the rotor of the Tesla turbine was performed 
with three—dimensional CFD simulations, based on compressible Reynolds averaged 
Navier-Stokes equations by the ANSYS Fluent software. The spatial discretization 
method utilized is finite volume with cell-centered approach and the mesh generated 
is hexahedral with local refinement. 


4.1 Air Tesla turbine 


The CFD analyses were performed in order to determine the flow field within the 
disks of the air Tesla turbine prototype. For this purpose, three—dimensional, double 
precision, pressure based, steady and implicit simulations were set. Velocity formu- 
lation was considered in the absolute frame and both laminar and transitional k-kl— 
omega models were analysed. The k-kl-omega was selected as turbulence model in 
order to evaluate the transitional behaviour ofthe flow, given that in this scheme tran- 
sition is not fixed but triggered by velocity fluctuations in the boundary layer. 

The pressure-velocity coupling scheme was set as SIMPLE, with second order 
upwind scheme for momentum and PRESTO! scheme for pressure equation. The ge- 
ometry of the model was created with the software ICEM on a disk with outer diam- 
eter 125 mm, inlet diameter 32 mm and gap between disks of 0.3 mm. A real fluid 
model was considered for air, with Peng-Robinson scheme enabled. 

The parameters set for the developed simulations are summarized in Tab. 3.17 
and the computational domain is shown in Fig. 3.38 where the inlet, the outlet, and 
the periodic boundary conditions are specified. The inlet mass flow was imposed as 
boundary conditions on the inlet, setting also the inlet flow angle. Pressure field was 
imposed on the outlet. The rotation feature of the rotor was applied by considering 
frame motion with fixed frequency of rotation. 


Table 3. 17 Parameters set for ANSYS simulations of air Tesla prototype. 


Fluid Air (real) 
Channel mass flow rate 0.001119 [kg/s] 
Outlet Pressure 131234 [Pa] 
Rotational velocity 18000 [rpm] 
Inlet flow angle 87 [?] 
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Fig. 3. 38 Three—dimensional computational domain of the Air Tesla turbine rotor (very coarse 
mesh) 


A mesh sensitivity analysis was carried out for CFD simulations, using both lam- 
inar and k-kl-omega turbulence models. Several meshes were created, with the ob- 
jective of investigating the effects of elements size and local refinement. In each case, 
the y* was selected to be lower than 1 for being able to evaluate transition effects 
when the suitable turbulence model is applied; the mesh was refined in particular at 
inlet and near walls, where the laminar separation was expected to occur. The imposed 
cell growth ratio was always set to be lower than 1.1. 

Tangential velocity was selected as parameter for mesh independence assessment, 
as this parameter proved to be the most affected by the grid characteristics, especially 
for transitional simulations (Fig. 3.39, where only significant meshes are displayed 
for the sake of clarity). In Fig. 3.39c, it is noticeable that a mesh with at least 800000 
nodes ensures grid independence, while coarser meshes determine errors in absolute 
tangential velocity computation, both at inlet and at outlet. 
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Fig. 3. 39 Results of the mesh sensitivity analysis performed on simulations with transi- 
tional turbulence model 


4.2 ORC Tesla turbine 


Rotor analysis 

As for the air Tesla turbine, the CFD analyses were performed in order to assess 
the flow field behaviour within the disks of the ORC Tesla turbine prototype The 
governing reduced Navier-Stokes equations are discretized with a cell-centered 
method with PRESTO! scheme for pressure, third order MUSCL scheme for density, 
momentum and energy and second order upwind for turbulent kinetic energy, specific 
dissipation rate and momentum thickness. To compute the variable gradients, both the 
Green-Gauss cell-based and the least-squares cell-based methods were applied. The 
turbulence closure was done both with the laminar model and the Langtry-Menter 
transitional shear stress transport model (SST) with second order discretization. 
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In this case k-kl-omega model could not be applied due to the sharp pressure 
gradients at rotor inlet. The equations were implicitly solved by applying the SIMPLE 
approach as pressure correction scheme. 

The gas density was calculated by Helmholtz Free Energy (FEQ) EOS, which is 
valid for different fluids (Tillner-Roth and Baehr, 1994, Lemmont and Jacobsen, 2000, Lem- 
mon and Hacobsen, 2005, Lemmon and Span, 2006). The properties of mixtures like R404a 
were calculated on the basis of those of each single component by using the ideal 
mixing rules; the properties of pure fluids R134a, R245fa and R1233zd(E) were ob- 
tained from literature (TilIner-Roth and Baehr, 1994, Lemmon and Span, 2006, Mondéjar et 
al., 2015). The main properties of pure fluids are summarized in Tab. 3.18. 


Table 3. 18 Main properties of the pure fluids 
MW 


Fluid [kg/kmol] Т. [K] P. [Pa] Pc [kg/m?] 
R134a 102.03 374.21 4.06-10* 511.90 
R143a 84.04 345.86 3.76:106 431.00 
R125 120.02 339.17 3.62-10° 573.58 
R245fa 134.05 427.16 3.65-10° 516.09 
R1233zd(E) 130.5 438.8 3.57-10° 478.9 


The computational domain is shown in Fig. 3.40 where the inlet, the outlet, and 
the periodic boundary conditions are specified. The inlet velocity and the outlet pres- 
sure fields were imposed as boundary conditions on the inlet and the outlet, respec- 
tively. The periodic boundary conditions were introduced for the sidewalls of half a 
quarter of a cylindrical domain in order to speed up computations. The rotation feature 
of the rotor was applied by considering frame motion with fixed frequency of rotation. 
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Fig. 3. 40 Three—dimensional computational domain of the Tesla turbine rotor 


Several diverse meshes, differing on the number of nodes (from 37000 to 1.87 
million), were examined in order to assess the mesh independency of computational 
results. Since tangential velocity is the parameter which is the most dependent on 
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mesh size, particularly at domain inlet, Fig. 3.41 resumes the results of the mesh sen- 
sitivity analysis in terms of inlet tangential velocity and computational time ratio ver- 
sus computational cell number. In order to attain a proper compromise between accu- 
racy and computational time, the 273000 nodes mesh was finally chosen. 
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Fig. 3. 41 Mesh sensitivity of inlet tangential velocity at rotor inlet and computational cost of 
accuracy 


Stator — Rotor analysis 

In order to assess the interaction between stator and rotor, as well as to analyse 
the sources of inefficiency that are triggered in the gap between them, a separated 
approach analysis does not provide accurate results; therefore, a coupled stator—rotor 
simulation was set. 

The numerical analyses of the Tesla turbine with coupling of stator and rotor were 
performed with the same settings specified for the single rotor analyses. The turbu- 
lence closure was still obtained with the Langtry—Menter transitional shear stress 
transport model (SST) with second order discretization. The computational domain is 
displayed in Fig. 3.42. The rotor still presents the same outlet and periodic boundary 
conditions, while the inlet has been specified in the plenum chamber. Indeed, the ge- 
ometry is composed by the stator, the rotor, the plenum chamber and the gap between 
stator and rotor. The inlet and the outlet pressure fields were imposed as boundary 
conditions. The periodic boundary conditions were introduced for the sidewalls of a 
quarter of a cylindrical domain in order to speed up computations. The rotation feature 
of the rotor was applied by considering frame motion with fixed frequency of rotation. 
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Fig. 3. 42 Three-dimensional computational domain of the stator—rotor configurations 


Several diverse meshes, differing on the number of nodes (from 108,100 to 4.45 
million), were examined in order to assess the mesh independency of computational 
results. With the aim to achieve the optimum grid size as a compromise between ac- 
curacy and computational time, the 3.23 million nodes mesh was finally selected, by 
fulfilling the minimum change of rotor inlet tangential velocity of 1*10 m/s. 


5. Test benches setup 


5.1 Air Tesla turbine 


In order to test experimentally the performance of the air Tesla turbine, it was 
arranged on a test bench (Fig. 3.432) equipped with several measurement sensors. The 
turbine was connected to an electric motor (brushless type) that controlled the rota- 
tional speed during the tests through a servo drive. By means ofa continuous exchange 
of electricity between the electric motor and the grid, the brushless motor was able to 
operate as generator or as brake, according to the power produced by the turbine. The 
torque and, consequently, the power produced by the turbine was measured by a 
torque meter (Lorenz, nominal torque 10 Nm, 0.2% FS) connected between the tur- 
bine and the brushless motor through flexible couplings. 

The Tesla turbine was tested with air as working fluid elaborated by a centrifugal 
compressor, whose speed could be regulated with a variable—frequency inverter. Fig. 
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3.43b shows the schematic of the air test rig composed by the compressor, the Tesla 
turbine and the measuring instruments. 

The thermodynamic conditions upstream the Tesla turbine were measured by in- 
stalling a T-type thermocouple and a pressure transducer (Honeywell, 100 psi, 0.1% 
FS) in the space between the casing and the stators. At the outlet, another thermocou- 
ple was installed to measure the temperature downstream the turbine, whereas the 
pressure was imposed as the ambient pressure due to the direct discharge in the envi- 
ronment. The mass flow rate was measured by using a calibrated flange directly con- 
nected to the pipeline between the compressor and the Tesla turbine. The temperature 
before the flange was measured by a T-type thermocouple (copper and constantan 
wires, uncertainty of + 0.5 K) and the pressure difference across the flange was ac- 
quired by a Setra pressure transducer (—5 to 5 PSID, 0.196 FS). 

The experimental measurement chain, the data acquisition and control system are 
showed in Fig. 3.44. All signals were acquired by a National Instruments Field Point 
and processed by a software specifically developed in LabView® environment, which 
was also used to control the servo drive and set the Tesla rotational speed. 

In addition to obtain the thermodynamic boundary conditions, the software com- 
puted the power produced by the expander by the measured torque and the imposed 
rotational speed. The mass flow rate was obtained by elaborating the measured pres- 
sure difference across the calibrated flange and the temperature upstream by using the 
equations implemented in the software after a previous calibration of the flange. 
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Fig. 3. 43 Test bench setup: a) Mechanical connection; b) Test rig schematic 
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Fig. 3. 44 Measurement chain and control system of the test rig 
5.2 ORC Tesla turbine 


Université de Liége test bench 


The schematic of test bench of Université de Liége is displayed in Fig. 3.45. It is 
basically a recuperative organic Rankine cycle, utilizing R1233zd(E) as working 
fluid. It was constructed using standard mass produced components from the HVAC 
industry, as well as from some specific industrial prototypes, such as the evaporator. 
The test bench consists of an evaporator, a condenser, a brazed plate recuperator 
(which however was by-passed during these tests), a pump and a liquid receiver. The 
pump is a Wanner Hydra—Cell piston pump, which is controlled via a frequency in- 
verter. The maximum volume flow rate and outlet discharge pressure are 30.6 l/min 
(at 1450 rpm) and 103 bar (at 750 rpm), respectively. The heat source is an industrial 
heater, which allowed a maximum heat input of about 150 kW. 

The thermodynamic conditions upstream and downstream the Tesla turbine were 
measured by installing T-type thermocouples (copper and constantan wires, very sta- 
ble measuring range between —200 °C and 200 °C, with a maximum error of + 0.5 K) 
and piezoresistive pressure transducers (Keller, 30 bar, 0.2596 FS). The mass flow rate 
was measured by using a precision Coriolis Krohne Optimass 1400Cr flow meter. 

The turbine was connected to an electric motor (Perske D 6800 Mannheim 1) that 
controlled the rotational speed during the tests through a four-quadrant frequency in- 
verter (ABB ACS501-01). The torque and, consequently, the power produced by the 
turbine was measured by a torque meter (Messtechnik DRBK, nominal torque 50 Nm, 
0.5% FS) connected between the turbine and the motor through flexible couplings, as 
displayed in Fig. 3.46. All signals were acquired by a National Instruments platform 
and processed by a software specifically developed in LabView® environment. Tab. 
3.19 resume the measurements ranges and accuracy of the sensors. 
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Fig. 3. 45 Schematic of Université de Liége test bench (recuperative ORC) 


Table 3. 19 Sensors ranges and absolute accuracies 
Quantity Unit Range Maximum error (full scale) 
73- 0.5 
K 
Temperature [K] 473 
Pressure [bar] 0-30 0.3 
Mass flow rate [kg/s] 0-1 2.0 E-3 
Torque [Nm] 0-50 0.25 
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Fig. 3. 46 Mechanical connection of ORC Tesla turbine to electric motor 
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Chapter 4 
Analysis of Results 


The results obtained in this project are resumed in this chapter. The chapter is 
organized in three main Sections, (1) 2D EES model, (ii) computational fluid dynamics 
and (iii) experimental results. 

All types of analysis are fundamental for understanding the principle of operation, 
as well as the research and industrial prospects of the Tesla turbine. The numerical 
model developed in EES allows evaluating the performance, varying lots of degree of 
freedom of the turbine, such as the geometry or the thermodynamic conditions. Com- 
putational fluid dynamics permits to understand flow characteristics, which are diffi- 
cult to simulate in a 2D code; while experimental results are the true litmus paper in 
order to confirm the validity and the perspectives of an emerging technology. 


1. 2D model results 


1.1 Air Tesla turbine 


The results obtained with air as working fluid have been computed in order to 
assess the performance potential of the rotor, as the parametric study displaying each 
component influence on the turbine efficiency and power will be deeply developed in 
the next Section, when organic working fluids will be considered. Nonetheless, a par- 
ametric analysis on the performance of air Tesla turbines as function of the main de- 
sign variables, which are non—dimensionalized following common practice in tur- 
bomachinery (Whitfield and Baines, 1990, Dixon, 2005), has been carried out. 

Rotor output/input diameter ratio (D3/D2) and rotor channel width over inlet rotor 
diameter (b/D2) were recognized as the most significant design parameters for the 
Tesla turbine. Exit kinetic energy and absolute exit flow angle, which should be as 
low as possible, were identified as two critical parameters to assess as they deeply 
influence the performance of the turbine. The trend of rotor efficiency 7 as function 
of (D3/D2) and (b/D3) is shown in Fig. 4.1. 

Rotor efficiency 7 is notably affected by different values of outlet/inlet diameter 
ratio. With decreasing D3/D2, the larger kinetic energy at discharge, due to the higher 
axial component of velocity, appears to be somewhat compensated by the larger rotor 
surface available for momentum exchange between the fluid and the disks. Nonethe- 
less, an optimizing configuration can be selected. Indeed, values of D3/D2 close to 1 
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do not allow the complete conversion of work, while values of D3/D2 close to 0.1 
increase kinetic energy losses at discharge. The right trade-off is therefore necessary. 

The channel width over rotor inlet diameter (b/D2) parameter strongly affects the 
efficiency of the rotor. This is due to the influence of the Reynolds number in the 
laminar flow regime, which can be best explained re-arranging its definition and re- 
membering that b = D,/2 (Eq. 4.1). 


Re = PU Me pun. ше (4.1) 
u 2mrbp ы Tr È 
0.8 


0.4 


9, um > 
0.001 0.0015 0.002 0.0025 0.003 0,0035 0.004 0.0045 


b/D; 


Fig. 4. 1 Rotor Efficiency n vs D3/D» and b/D; 


The load coefficient behaviour ¥ as function of (D3/D2) and (6/D2) is shown in 
Fig. 4.2. The inlet rotor diameter (5/D») strongly affects the load coefficient % fur- 
thermore, the load coefficient is also sharply influenced by (D3/D2). The momentum 
exchange is favoured as the wet area is increased (values of Ds/D» close to 0.1); nev- 
ertheless, the exit diameter should not exceed a certain limit, the penalty being an 
increase of the residual tangential velocity, leading to higher discharge losses. From 
the points of view of rotor efficiency 7 and load coefficient У, values of 0.35 < 
(D3/D2) < 0.45 and 0.005« (b/D?) < 0.015 appear therefore recommendable. 

Another essential parameter, which is definitely influenced by the gap between 
disks, is the absolute exit angle (Fig. 4.3). Lower absolute exit angles are desired in 
order to have an efficient recovery of discharge kinetic energy. A reduction of the gap 
between the disks is certainly beneficial to this end. The decrease of оз for smaller 
values of b is due to the reduction of the tangential component, as well as to the in- 
crease in the radial component of absolute velocity. The reduction of tangential ve- 
locity, as can be noted from Eq. (3.47), is due to the increase of viscous momentum 
transfer for small values of b. On the other hand, the radial velocity increases because 
of the continuity (Eq. 3.25). 


124 


Lorenzo Talluri 


0.3 
0.001 0.0015 0.002 0.0025 0.003 0.0035 0.004 0.0045 


b/D; 


Fig. 4. 2 Rotor Load Coefficient V vs D3/D; and b/D2 


0.001 0.002 0.003 0.004 0.005 
b/D, 


Fig. 4. 3 Exit fluid angle оз and efficiency of the turbine vs b/D2, for D3/D2 = 0.44 


The trend of the absolute velocity and of its components with variable gap is 
showed in Fig. 4.4. The absolute exit velocity displays a minimum, which is deter- 
mined by opposite trends of the two components of velocity (radial and tangential) as 
a function of b/D2. The minimum of the exit kinetic energy Ex;nz corresponds to the 
maximization of the rotor efficiency. 
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Fig. 4. 4 Exit kinetic energy vs non dimensional gap 


Fig. 4.5 displays the absolute exit angle and the efficiency behaviours against the 
flow coefficient Ф. An increase in the flow coefficient leads to a decrease of the ab- 
solute exit angle. This is due to an increase of the radial component of the fluid, which 
therefore turns the fluid in the axial direction. If values of exit flow angles below 50? 
are sought, then a flow coefficient in the range Ф = 0.2 should be selected; under these 


conditions, the rotor efficiency is still high — in the range of 0.94. 
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Fig. 4. 5 оз and n versus Ф 


The analysis of the rotor utilizing air as working fluid allowed confirming the 
possibility of achieving high efficiencies when the turbine is properly designed. Next 
Section will be focused on the influence of not only the rotor but also of the other 
components on the performance of the turbine, taking n-Hexane as working fluid. 
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1.2 ORC Tesla turbine 
1.2.1 Component analysis 


In order to assess the performance potential of the Tesla turbine for organic Ran- 
kine cycle applications, several parametric analyses were carried out: the performance 
parameters were evaluated as functions of the main geometric variables and operating 
conditions. The n-Hexane was adopted as the working fluid for the investigation of 
each component influence on the performance ofthe turbine, due to its favourable low 
expansion ratio and the well suitable thermodynamic critical conditions (e.g. low crit- 
ical pressure 3.034 MPa and high critical temperature 234.67 ?C (Aljundi, 2011)). Fur- 
thermore, hydrocarbons are among the best compromises between environmental con- 
straints (having zero ODP and СҰР and fairly low toxicity (Aljundi, 2011)), and good 
thermodynamic cycle features for use in ORCs. The only drawback might be the flam- 
mability, but the rapidly increasing safety standards, coupled to typically small or mi- 
cro sizes for which Tesla expander technology is generally conceived, make hydro- 
carbons particularly attractive for these applications. Initially, single variable optimi- 
zation was carried out, in order to determine the most critical parameters for the Tesla 
turbine; successively combined parameters optimization was performed, to under- 
stand the mutual influence of the most meaningful parameters both on efficiency and 
power. 

Individual variables optimization 

In this Section, the influence of each single parameter on the performance of the 
turbine is analysed, while keeping constant all the other geometric and thermody- 
namic parameters. 

Rotor inlet diameter 

The rotor inlet diameter is one of the most significant parameters, as it plays a 
primary role on both the power production and the size of the machine. When fixed 
thermodynamic conditions and velocity are assumed at the stator output, a higher rotor 
diameter is associated to a larger throat section. Therefore, according to continuity 
equation, a higher mass flow rate is obtained. The increased of mass flow rate is also 
responsible for an increase of velocity at rotor output (as radial velocity will strongly 
increase as is progress to the exit). It leads to higher kinetic energy losses, which be- 
come unbearable when the outlet Mach number value overcomes the inlet one. High 
kinetic energy losses are not acceptable; therefore, the value of D2 is limited. The work 
output per unit mass of the expander is defined by the Euler equation Eq. (4.2). 

Work = vi? ' U2 — Vi * U3 (4.2) 

Since both the first and second terms increase with larger rotor diameter (increase 
in peripheral velocity), the specific work output presents an optimization value, while 
the overall power output, mainly influenced by the mass flow rate, has a monoton- 
ically increasing trend. The rotor and the total efficiencies are affected by increasing 
the rotor isentropic enthalpy drop and mass flow rate, thus their values decrease in- 
creasing rotor diameter, as shown in Fig. 4.6. 
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Fig. 4. 6 Turbine efficiency a) and turbine power and losses b), versus rotor inlet diameters 


Camber line length 

The length ofthe camber line in the stator channel is another parameter that affects 
turbine performance. An increase of the camber line, while keeping the other geomet- 
rical parameters fixed, is responsible for the reduction ofthe throat section width and, 
consequently, of the mass flow rate. According to Eq. (4.2), the flow velocity at the 
rotor outlet decreases and generates a higher work output. Therefore, an increase in 
the length of the camber line has the same effect of a reduction of the rotor inlet di- 
ameter (Fig. 4.7). 
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Fig. 4. 7 Turbine efficiencies a) and power and losses b), versus length of nozzle camber line 


Height and width of the nozzle throat section 

The height and the width of the throat section directly influence the geometry and 
the variables at the stator output, especially the mass flow rate, which shows a linear 
trend. The maximum value of mass flow rate is determined by the sonic condition at 
the throat section: decreasing the height of the rectangular area, the continuity equa- 
tion ensures a reduced flow rate and, consequently, a lower expansion rate inside the 
rotor. The velocity components are both proportional to the flow rate, Eqns. (3.25) — 
(3.27), so the Mach number increases rapidly, thus increasing the overall efficiency 
(Fig. 4.8). 

The reduction of throat width is responsible for a more than linear decrement of 
flow rate and velocity at rotor outlet that results in a reduction of power output and 
kinetic energy losses at rotor exit, whereas the expander efficiency shows a linear 
increment. 
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Fig. 4. 8 Turbine efficiencies a) and power and losses b), versus height of the throat section 


Rotor Channel height 

The channel height, which is present in Eq. (3.27), is directly proportional to tan- 
gential velocity and inversely proportional to radial velocity. It leads to a maximiza- 
tion of the rotor outlet Mach number and, accordingly, to a minimization of outlet 
rotor pressure. This value of channel height corresponds to the one optimizing both 
rotor and turbine efficiency, as lower velocities are associated to lower kinetic energy 
losses (Fig. 4.9). 
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Fig. 4. 9 Turbine efficiencies a) and power and losses b), versus disk channel height 


Outlet rotor diameter 

The variation of D; presents a significant minimum value of tangential velocity at 
rotor outlet, as can be deducted from Eq. 3.27. Similarly, to the previous analyses, the 
optimizing value of outlet/inlet rotor diameter ratio Ds/D? is characterized by the min- 
imum energy loss and, therefore, the maximum turbine efficiency (Fig. 4.10). 
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Fig. 4. 10 Turbine efficiencies a) and power and losses b), versus in/out rotor diameters ratio 


Thermodynamic conditions 

The effects of the thermodynamic conditions on the turbine performance were 
also considered; particularly, total conditions at the stator inlet, for a determined su- 
perheating level were evaluated. High pressures correspond to high flow rates (con- 
versely to high temperature, which partially reduce it due to a decrease in the fluid 
density) and subsequently to high expansion rates. As reported in (Sengupta and Guha, 
2012), there are specific thermodynamic conditions, which generate negative relative 
velocities at rotor inlet, then a flow reversal (these are linked to the actual velocity at 
throat section). Therefore, for a particular value of inlet total pressure Poo, the relative 
velocity becomes zero. In this condition, the rotor inlet radial velocity component V2 
has the minimum value and maximizes turbine efficiency (Figs. 4.11a and 4.11c). 

The static pressure at the stator output determines the pressure range into the con- 
vergent channel and the mass flow rate. The effect of increasing P2 (which correspond 
to lowering the mass flow rate) is the same as a reduction of Poo, therefore an effi- 
ciency optimizing working point for the machine can be found (Figs. 4.11b and 
4.11d). 
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Fig. 4. 11 Turbine efficiencies versus a) Total inlet pressure, c) Stator outlet static pressure and 
power and losses at various b) Total inlet pressure; d) Stator outlet static pressure 


Rotor peripheral speed 

The variation of the rotor peripheral speed influences the relative tangential ve- 
locity wx, thus a value of peripheral speed above which a reversal flow condition is 
generated exists. At disk exit, the pressure decreases when the rotational speed in- 
creases, due to momentum equilibrium in radial direction (similar to the behaviour of 
pumps where Н = rpm?). Rotational speeds in the range between 4000 and 6000 
rpm seemed to guarantee an optimised value of the expander efficiency, as displayed 
in Fig. (4.12). 
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Fig. 4. 12 Turbine efficiencies a) power output and losses b), versus rotational speed 


Combined variables 

Single parameter optimization is necessary in order to understand the influence of 
each element on turbine performance, but it is not enough in order to achieve a com- 
plete optimization procedure due to the mutual influence of each variable. As reported 
in Fig. 4.13, the combined assessment of both inner and outer disk diameter displays 
an improvement of the expander efficiency (7) at reduced values of stator outlet di- 
ameter D; and rotor outlet diameter D3. This trend is directly linked to the reduction 
of the mass flow rate; power output conversely exhibits an opposite behaviour. 
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Fig. 4. 13 Turbine efficiency a), power output b), between stator—rotor losses c) and kinetic 
energy at rotor outlet d) versus stator output diameter at variable in/out rotor diameter ratio 


A further assessment was developed taking into account the influence of the di- 
ameters and adding an additional parameter, namely the camber line length, which is 
directly linked to the external radius value. In agreement with previous results, the 
efficiency is maximised when the throat section is at the minimum value; it can also 
be pointed out that long camber lines allow better results when coupled with low 
D3/D; (Fig. 4.14). These conditions allow the achievement of relatively high overall 
turbine efficiency (over 60%). 
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Fig. 4. 14 Turbine efficiency a), power output b), stator—rotor losses c), and kinetic energy at 
rotor output d) versus stator exit diameter, for different rotor in/out diameter ratio and camber 


line length 


A comprehensive study of the stator should also include the camber line length, 
the throat length (which are function of the total maximum number of feasible chan- 
nels Zax) and the actual number of channels (Zs). The throat section is directly linked 
to the total number of channel, decreasing quadratically as Zax increase. The curves 
in Fig. 4.15 show that the optimal conditions are obtained for higher turbine dimen- 
sions (ro and D7), due to a reduction of the throat section (because the number of chan- 
nels is maintained constant) and therefore a reduction of mass flow rate. Furthermore, 
Fig. 4.15 displays that for a fixed external 0.35 m diameter, the maximum efficiency 
can be kept constant by reducing both the external radius (larger turbine size), and the 
length of the camber line (therefore keeping constant the throat section). As discussed, 
this is due to the counterbalancing effects, which maintain the mass flow rate constant. 
In this way, a Tesla turbine may achieve a more compact shape and the throat length 
can be reduced. A higher number of stator channels (Zs in Fig. 4.15) is responsible for 
steeper curves and lower global efficiency at fixed values of ro and Zmax. 


a7 0,7 


a) Z=2 b) Z= 10 
né - м so 06 
" Li м 
as к^ um Ў 4, а и 0,8 " " 7 
ар F4 dl т s.s. * [Р 
Z o4 А 4 = 04 à 4 " P а 
- e» fø = ^ à hi Н P 
s Р ^" o» E - ^ 
as .. , 03 fo; . 
é WA 4 Zeax720 —*-D,70.15 ee è n . 
9 ? Zea —а—ру=0.25 . 33 3; Zes720 ~e- Dy=0.15 
Or dé 4 Zea)  —- 29.35 03 2 $ Zax —e-Dı=0.25 
Tax  —— D90.36 
ол ол 
aos 0.1 0.15 a2 9.25 os 0.35 оа 0.05 0.1 0,15 9.2 0.25 as 0,38 оа 
ro [m] ro [m] 


133 


Micro turbo expander design for small scale ORC 


400 400 
с) „^^? d) 07 10 
Za) —9-D,29.45 12—20 —9-D,70.15 
300) 7:40) * Dyes 300) ә +*-D,0.25 
= Zea c 0,=0.35 69 se Dj,0.35 
- 2% È 
= ч ч м — 19 м Ы м 
5 20 2 200 Ú > DI 
LJ Г 
2 1 aa . " PELO 44 1 a 1 
E Е i à . 1 * ^ 
100 LÀ A_N М È 100 à à e à 
^ a м pi 1 E 
ч oe à be à 
we . el n “ La ч AN % à A 
о he Ses n n х Pi 
о 01 615 62 ом оз 036 ол оз 01 оз 02 O28 оз 935 Of 
ro [m] ro [m] 
x 10% x 10 
094 094 
e) А è è t "T 
092 ? ? ^ 097 H 2 3 ^ 
u ju РА ti PG i y [v " 
. f 2 - - 
CUI LI sá | es oa К D my 
— ” -— 
= dl è da = 
& ass . a a ^P em њм 32 & onn 
i i 
e oss, Z,=2 ei ass Z= 10 
4 Zes10 —® 03-915 4 Zea) е Dy @tS 
ons, 7==“4%0 —-D,"025 ояд 7==”40 —ә-Юу"@25 
Zari = [у-@34 Jan  —- D,20.34 
os os: 
9.05 Qi os 92 0.25 °з oss a4 ens [3] os a: ors as °з a4 
ro [m] ro [m] 
xm y 19* 
os as 
2) 472 "20 -9 Di=@.15 h) 7,=10 "20 Dats 
a? fe аа Pi Zea40 —è- Dyes от s e 740 è Dy=0.25 
è Лә“ Dy 34 + a lan"60 2 "M035 
" - 7 - [TI 
- ~ " 
- M5) 4? . 
= — as 
© оз A Es 14 è Y - 
= Ы ^ ^ T u 4? . 
fo E ES ^ 3 = i» 4 
[| 1 м Hi è 5 R 
uw M ч a^ gu .* $ LL 1 
éh * m VM м ы è . 
di 43 WE “м” = a2 oe .* "1. - 
Ste ‚* * * "P A 
° 01 
oo Of as OF OF өз өз оа 6605 dì өз OF ом өз өз оз 
ro [m] го [m] 


Fig. 4. 15 Turbine efficiency, power, losses between stator and rotor, kinetic energy at rotor 
outlet at different stator inlet and exit diameters, with 2 nozzle channels (a, c, e, g) and 10 
nozzle channels (b, d, f, h) 


The analysis of the nozzle throat height can be performed while taking also into 
account the rotor channels height b and the thickness of the disks s. This study was 
performed using a constant value of throat width (TW=1 mm) and 10 channels in the 
rotor. The minimum values of b and s were chosen taking into account potential eco- 
nomical manufacturing processes and possible structural issues. The curves present a 
turbine efficiency optimizing value as function of b for both the analysed rotor diam- 
eters (0.8 mm for D;=0.15 m and 1 mm for D;=0.35 m), while the increase of disks 
thickness is responsible for a wider throat section and a higher mass flow rate (higher 
power output), which entail a reduction of turbine efficiency 7 (Fig. 4.16). 
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Fig. 4. 16 Turbine efficiency a), power output b), stator—rotor losses c), kinetic energy at rotor 
outlet d) at various channel height b for different values of plates thickness s 


The change of the number of rotor channels (and consequently H;) does not mod- 
ify the behaviour of the curves and the optimising value of Б, although a reduction in 
the number of rotor channels лел can significantly improve the overall efficiency of 
the turbine (Fig. 4.17). This is directly linked to a reduction of throat section; there- 
fore, lower mass flow rates are obtained at fixed thermodynamic conditions. The use 
of just 2 channels for every module allows reaching efficiency values (7) near 60%, 
also reducing the influence of the disks thickness. On the other hand, the reduction in 
number of channels implies a drop of power output, which can be counterbalanced by 
utilizing a higher number of modules. 
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Fig. 4. 17 Turbine efficiency a), power output b), stator — rotor losses c), kinetic energy at rotor 
outlet d) at various plates thickness s for different channel heights b and three configurations 
with different total number of channels Nen 


Collecting the result of the performed sensitivity analysis guidelines for design of 
an optimized and balanced Tesla turbine can be drawn; these may be summarised in 
the following: 

e Stator, 4 convergent nozzles with a squared throat section (1x1 mm); 

e Rotor, 10 channels 0.1 mm wide each; disks diameters ratio around 0.4; the 

external radius was not fixed, in order to evaluate the effects of the turbine 
size. 


After fixing the geometry of the turbine a sensitivity analysis as function of ex- 
pander operating conditions was carried out. Originally, the superheating temperature 
and the pressure drop AP inside the stator were fixed: the minimum value of total 
pressure Poo and total temperature Тоо were set, in order to avoid sonic condition at 
the stator outlet. An optimal value of P; was found, as the turbine efficiency is affected 
both by the mass flow rate and the enthalpy of the fluid: the former decreases with 
reducing Po (leading to a positive effect on the turbine efficiency 7), whereas the 
latter has the same effect of Тоо (positive when the temperature increases). The same 
considerations can be extended to the power output, which shows an optimizing value 
at lower total inlet pressure, due to the increasing density and mass flow rate (Fig. 
4.182). 


0.226 +» x 19* 
" ,—9. . 
a) pray” Sua 7000 rpm ' 3 
озм! * Р . . 9.5 b) as 7000rpm + * 
t * * = ө . .* 
Peel . . . F1 as; „е 
m а . » = PE 
ya . е я = Fas 
= 022 E . ns E 2 
Lj 
ч Zg е Nd 
0.218 d Е «Li re 
E m as aat *- AP 
i a .* : 
0.216 ws 02 „-* o- Einon 


2.5 3 AS 4 & as 6 64 7 2 25 5 3.5 4 ss ^ 65 1 


45 3 3 45 5 
Р; [Ра] x 10° P, [Pa] x 10° 


Fig. 4. 18 Turbine efficiency and power output a) and losses b) versus static pressure at stator 
outlet (D;= 0.15 m) 
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The sensitivity analysis to the rotational velocity showed that the expander effi- 
ciency 7 and power increase at higher speeds, with a slight peak placed at lower pres- 
sure (the limited extension of the curve at higher rpm is due to the high Mach number 
reached at the rotor outlet, Fig. 4.19). 
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Fig. 4. 19 Turbine efficiency a), power output b), and losses c) versus static pressure at 
stator outlet for variable rotational speed (Dı 0.15 m) 


In Fig. 4.20, the expander efficiency and the power output are presented as func- 
tion of rotational speed and stator inlet total pressure Poo. For all different Poo, an 
optimised value of rpm is always present, lower at higher total inlet pressure; effi- 
ciency is enhanced at lower Poo, due to higher rotational speed and lower absolute 
tangential velocity at rotor outlet, which imply higher work production. The presence 
of a maximum may be explained through the momentum balance: higher velocity in- 
creases the work output, but, at the same time, also the expansion ratio and the en- 
thalpy drop through the rotor are enhanced. In this way, from the definition of total to 
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Fig. 4. 20 Turbine efficiency and power output a), losses b) versus rotational speed at variable 
stator static pressure drop (D = 0.15 m) 
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The superheating temperature level does not have a strong influence on Tesla tur- 
bine efficiency n when compared to the effects provided by pressure and rotational 
speed. 

A further analysis, which needs to be carried out, is the one to assess the influence 
of up-scaling the Tesla turbine dimensions (for example doubling the diameter), while 
keeping the outlet/inlet rotor diameters ratio and throat section width fixed. The per- 
formance curves obtained show a very close behaviour to the one, which resulted, 
from the previous results (Figs. 4.19 and 4.20) but with improved values (Figs. 4.21 
and 4.22). It was found that the up-scaled expander holds an optimized efficiency 
value for a static outlet pressure P; of 3.3 bar and a rotational speed of 10,000 rpm. 


at. 


ATS 


D 


150 
o-o-o. i b) 140 
per i O 130 
* . 
M E 
. B К 
. o-de? - 
ы oe" . — 110 
. s: 2 
"O . — 100 
1 ha = 
Е = "e EL 
. E = so 
. LJ 
16000 rpm . 79 
13000 rpm 7000 rpm e 
10000 rpm 5000 rpm 
so 
м 4 AS S SS 6 6S 31 
P, [Pa] s №? 
iw 
' 
9 as 
T as 
= [hj 
as! 19000 гри — 7000 rpm 
= 13000 rpm 5009 rpm Р 
as . 
* 10000 rpm у 
a u .. .* 
— D Lo 
po ia к 
$ a? s a” е 
3.9.5 -— 
- oe * раче 
Ot ssp fee 
° 
2 24 3 34 4 44 4 
P, [Pa] 


^ 


és 
x1 


LI 


34 


16000 rpen 

13000 rpen 
* 1000rpm 
eo 

ee e 
inn Д 
è 0000000 
9-9-9-9-—-06-0—9—9 
а 48% 4 44 
Р, [Ра] 


7000 rpm 
5000 гре 


Fig. 4. 21 Turbine efficiency a), power output Б) and losses с) versus static pressure at stator 
output, at variable rotational speed for the up-scaled expander (D; = 0.3 m) 


A very small stator static pressure drop AP, (in this case 0.5 bar) allows the 
reduction of inlet total pressure, optimising rotational speed and, on the whole, an 
increase of the turbine efficiency 7 up to 51%. 
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Fig. 4. 22 Turbine a) efficiency and power output, b) losses versus rotational speed rpm at 
various total inlet pressure, for the up—scaled (D;= 0.3 m) expander 


138 


Lorenzo Talluri 
Finally, taking into account the developed sensitivity analyses, it was possible to 
point out some useful guidelines for design and optimization of a Tesla turbine utiliz- 
ing n-Hexane as working fluid: 

e The expander efficiency, power output, mass flow rate and expansion ratio 
are in close relation: low mass flow rates m are connected to high efficiency 
and lower power output and vice versa, in agreement to the statement claimed 
in (Armstrong, 1952, Rice, 1965). 

e The length of the camber line L; and the number of the nozzles Zs are funda- 
mental geometric parameters, as they directly influence the mass flow rate, 
while the rotor variables (b and D3/D2) influence the variation of tangential 
velocity Av, and therefore the power production and the performance of the 
expander. They can be optimized in order to achieve optimal performances. 

e Generally, the best performance of the Tesla expander is achieved with low 
inlet pressure and limited mass flow rates. Through all conditions, an opti- 
mised value of rotational speed is present (the one which allows inlet tangen- 
tial relative velocity close to 0). The performance is not significantly affected 
by inlet temperature. 


Full design procedure example 

In this Section, an example of a Tesla turbine design procedure is displayed. The 
first matter to take into account, in order to achieve a complete and proper design, are 
the boundary conditions of the geometric parameters; particularly: 

e Balanced stator size (Do/Dix1.5); 

e Stator output angle a<85°; 

e Throat section length > 1mm; 

e Rotor channel width > 0.1 mm and disks thickness > 0.5 mm. 


Even the fluid conditions have constraints: 
е Ма;<1 in the throat section and Mas < Map. 


Applying the considered boundary conditions and selecting disks with an external 
diameter D5;-0.30 m, best efficiencies were obtained with a reduced throat section 
length (1 mm) and with a limited number of statoric nozzles (Z=2). Nonetheless, this 
choice leads to low power output values and a not uniform flow distribution at rotor 
inlet. Due to the applied conditions, the shape and the length of the statoric channels 
require an inlet/outlet diameter ratio Do/D, equal to 1.48, while the need of a tangential 
flow is favoured by a relatively large stator outlet angle. The parametric analysis sug- 
gested b=0.12 mm and s=0.5 mm as the optimum rotoric channel width and disks 
thickness respectively, in order to maximize turbine efficiency. In this way, if the rotor 
is made of 5 disks, the throat section height is 2.6 mm, while the diameter ratio is set 
to 0.2. 

A parametric study was thus carried out varying total inlet pressure Poo, static 
stator outlet pressure P; and rotational speed. Efficiency was maximised at very low 
inlet pressure (4 bar) and stator pressure drop (0.5 bar) and with a rotational speed of 
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10,000 rpm. These operating conditions are close to the incipient flow reversal at rotor 
inlet. This occurrence sets a limit on the values of rpm and APsia:, which, therefore, 
actually reduce the operating range. Nonetheless, very close values of efficiency can 
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Fig. 4. 23 Operating map of the Tesla expander: a) efficiency, b) power output, c) stator-rotor 
losses, d) kinetic energy at rotor outlet versus rotational speed at variable static pressure drop 
in the stator and total inlet pressure 


The curves show that the maximum achievable value of efficiency is close to 64%, 
achievable at various combinations of total inlet pressures and rotational speeds. In 
particular, raising the total inlet pressure allows lower rotational velocities in order to 
achieve the same efficiency levels. The effects of the inlet temperature are not very 
relevant and efficiency is roughly constant versus temperature (Fig. 4.24). The modest 


variation is only due to the higher-pressure losses in the gap betwe 
0.66 о) 
а) ah Po,73*10* [Pa] b) 
a . Pog? 4" 10" [Pa] è 
0.65 à Poo S*10° [Pa] so 
a A 
& . . . 
v 0.64 4 - 40 
e . ^ = 
Fa [1] "Uwe ~es а = Lai ы Ф e 
LI ^ 2 
4 LI LI = ал” . Z ш 
: *-Ma; = 1 x ы в *-Ma; *1 
OGT a- DP = 0.5 bar 20) a DP aia "0.5 bar 
a 
0.61 A 10 a “ау Aia 
100 1100 120 13M 140 150 160 170 189 190 100 110 120 1м 1 
To [°C] Fo 


140 


en stator and rotor. 


Pa73* ША [Pa] 
Рота 10° [Pa] 
Рь=5*10° [Pa] 


LIT 


1% 
[PC] 


160 150 199 


Lorenzo Talluri 


" 4 
xw x 10 
13 008 


9 d) A Poo"3* 10° [Pal 
i o—e E E E . 0.97} e Ma; =] od : гы а Pal 
* DP ae = 0.5 bar 2760 * [Pa] 
0.06 $ t 
EFT. . . E . . E u 5 А 
a = = 
= " otel mm. = ооз e н * " 
fad seee o . Pod" 10° [Pu] - Hs А 
«È Po 4°10" [Pa] йам . = 
= * Mayet Р==5"10* [Ра] E * A А 
04 * DPua 705 bar é : a & 
а a 5. 9 „** ate Q5 da aes . e а . 
a: on 
190 10 120 1390 140 150 160 170 (1596 190 100 310 3120 1596 14 1% 1640 1% 150 1x 
Too [°C] Teo [°C] 


Fig. 4. 24 Tesla expander operating map: a) efficiency, b) power output, c) losses stator—rotor, 
d) kinetic energy at rotor outlet versus inlet total temperature at different total inlet pressure 


In the present assessed case, an inlet 5 bar total pressure and the related saturation 
temperature are optimised conditions, while the static stator pressure drop AP; 
should be set at the value allowing the sonic condition in the throat section. The rota- 
tional speed showed an optimised value at 11,700 rpm. Under these design conditions, 
the Tesla expander achieves 54 W power output per channel and 61% isentropic effi- 
ciency. 

Down-scaling or up-scaling the Tesla turbine by modifying only the outer stator 
diameter D;, while fixing all the other geometric design parameters, changed slightly 
the turbine performance, as shown in Tab. 4.1. Particularly, efficiency increased due 
to the increase of power production and the decrease of the non—dimensional channel 
width ratio (b/D2). The most important feature to remark is the behaviour of optimised 
rotational speed when changing the diameter of the turbine. Particularly, high turbine 
diameters allow a reduced value of rotational speed. This matter is really important 
for some specific applications, (such as distributed power generation), as where the 
velocity could be limited for various reason (direct coupling to network frequency, 
noise, ecc.). The increase in diameter allows also a higher power production per chan- 
nel, requiring therefore a limited axial extension of the expander in order to provide 
relatively high power. 


Table 4. 1 Performance and geometric parameters of the assessed turbines 


Di [m] 0.15 0.3 0.5 
Lis [mm] 1 1 1 
H; [mm] 1.2 1.2 1:2 

Zs H] 2 2 2 
b [mm] 0.1 0.1 0.1 
D3/D: [-] 0.25 0.25 0.25 

RPM 23400 11700 7000 

o 0.03 0.02 0.01 
y 0.86 0.92 0.93 
Ns 0.029 0.013 0.008 
Ds 9.07 80.08 34.16 
n [%] 52 61 64 
Power per channel [W] 41 54 58 
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1.2.2 Turbine geometric assessment 


Once the influence of each parameter, and an optimization procedure for the de- 
sign and optimization of the turbine were developed, it was questioned if the main 
parameters could be linked together, in order to obtain geometric scaling laws, which 
are necessary when executing a performance assessment which aims to be independ- 
ent from the size of the turbine. All the main geometric parameters have therefore 
been analysed in order to obtain maximum turbine total to total efficiency. The scaling 
laws have been developed in order to link every geometric parameter to external rotor 
radius. 

The main considered parameters were found to be: 

e Stator inlet/outlet diameter ratio (Dg/D4); 

e Gap between stator and rotor (G); 

e Rotor channel width (b); 


e Rotor outlet/inlet diameter ratio (К = D2/D3); 


e Throat width ratio (TWR = Meistit 
2xT*r2*b*ngisk 
Stator inlet/outlet diameter ratio 
The stator inlet/outlet diameter ratio was fixed at 1.25. This value was assumed 
as suggested in (Glassman, 1976), where a consolidated practice relying on geometric 
definitions was carried out. 


Thus, the outlet diameter of the stator comes out: 
р; =D, +2°G (4.4) 


Gap between stator and rotor 
The gap was chosen as small as possible but limited by the thermal expansion of 
the disks: 


G = 1.5: (r2 — r3) A (Too — Tamb) (4.5) 


Where 4 is the thermal expansion coefficient, Тоо the total temperature at stator 
inlet (conservative assumption) and Tamb is the ambient temperature (this is a con- 
servative assumption, actually in/out rotor temperature should have been used). 
Therefore, the calculated value represents the linear thermal expansion of the disk, 
actually augmented by a 1.5 safety factor. 

Rotor channel width 


The rotor channel width was defined as a function of the external rotor diameter, 
by the means of an extensive parametric analysis. Fig. 4.25a shows the total to total 
efficiency of the turbine vs. channel width, at fixed 100?C total inlet temperature and 
total inlet pressure corresponding to a 10 K super heating level (or, in other words, at 
90°C saturation pressure; in order to compare all the different investigated fluids at 
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the same low temperature level). For the sake of clarity, Fig. 4.25a shows the results 
for a fixed 0.2 throat Mach number and a 0.4 rotor outlet/inlet ratio. The throat width 
ratio is fixed at 0.02 at 100°C temperature. Nonetheless, the analysis was performed 
at various throat Mach numbers (0.2, 0.3, 0.6, 0.9), various rotor diameter ratios (0.2, 
0.4, 0.6) and throat width ratios (0.02, 0.04). The effect of throat Mach number, rotor 
diameter ratio and throat width ratio do not remarkably influence the position of the 
best efficiency, but only the value (Fig. 4.25a). Therefore, Fig. 4.25b shows the loci 
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Fig. 4. 25 a) Total to total efficiency against channel width. b) linear interpolation of channel 
width against rotor external diameter at highest efficiency value 


The linear interpolated equation for the determination of the rotor channel width 
b, which allows achieving the highest efficiency vs. rotor outer diameter D» is reported 
in Eq. (4.6) for r1233zd(E) (equations for different fluids are presented in the follow- 
ing Sections). 
b = 0.0002: D; + 3: 10795 (4.6) 
Furthermore, the pressure effects are also assessed, in order to investigate the in- 
fluence of pressure on efficiency and suitable rotor channel width. Specifically, the 


super heating level is fixed at 10K and the total inlet temperature is changed from 75 
to 175?C. At low temperature, and therefore at low pressure, the rotor channel width 
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allowing best efficiencies is larger compared to what is required at high temperature 
(and therefore high pressures). This is directly related to the mass flow rate. Indeed, 
high pressure corresponds to high mass flow rates, which requires narrower rotor 
channel width in order to achieve high efficiency. Narrower rotor channels favour 
viscous entrainment, while higher mass flow rates enhance momentum balance. 
Therefore, in order to have a suitable power extraction, tight channels are required at 
high mass flow rates; otherwise, an enhancement of kinetic energy loss is present. 
Nonetheless, the difference in rotor channel width is quite small and values between 
0.00005 and 0.0001 m are recommended, in practice for all conditions. In order to 
have a comparable condition, the selected reference case is the one at 100°C for all 
the analysed fluids. 
Rotor outlet/inlet diameter ratio 

The best conditions for the rotor outlet/inlet diameter ratio (R) were evaluated 
running several parametric analyses (determining different Mach number conditions). 
It was found that when the optimal channel width correlation is applied, the best value 
for practically every turbine size is always in the range from 0.3 to 0.4, with the lower 
bound corresponding to low Mach number (0.3) and the higher bound to relatively 
high Mach number, close to 1. Fig. 4.26 displays the total to total efficiency as a 
function of R at various Ma numbers (0.3, 0.6 and 0.9) for the lower (0.08 m) and 
upper (0.44 m) diameter range bounds considered in this analysis. Smaller turbines 
can achieve higher efficiency at the price of a lower power production and higher 
rotational speeds. In the present case study, a value of R = 0.35 is selected, which 
guarantees good efficiencies at every investigated Mach number. 
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Fig. 4. 26 Total to total efficiency against rotor outlet/inlet diameter ratio 


144 


Lorenzo Talluri 


Throat/width ratio 
The throat—width ratio was determined as the ratio between stator outlet area and 


rotor inlet area, as shown in Eq. (4.7). 


TW-Hs'Zstat 
2t r2:b:ngisk 


IWE- (4.7) 

Where TWR is the throat width, Н; is the throat height, Х,а is the number of noz- 
zles; naisk is the number of rotor channels per each single stator. 

It was found that the total to total efficiency increases monotonically with reduc- 
ing TWR (Fig. 4.27). Actually, at lower TWR corresponds higher velocities at the 
throat, which are beneficial for rotor efficiency. On the other hand, the power output 
shows an opposite behaviour, steadily increasing with widening TWR. Therefore, a 
balanced solution was selected, that is, TWR=0.02. The balanced solution takes into 
account both thermodynamic matters (high efficiency at a suitable power) and manu- 
facturing issues: in fact, lowering 7WR implies dealing with manufacturing of very 
small stator channels, with a relatively complex geometry and compactness of the 
expander. Indeed, a very low TWR would certainly increase the efficiency, but would 
require a large expander: compactness factor CF — Power/Volume lower than 0.3 
W/cm? were assessed for TWR = 0.01. 
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Fig. 4. 27 Total to total efficiency and Power vs TW ratio for various Tesla turbine dimensions 


Fluid dynamics assessment 

Besides the investigation of geometric parameters, a fluid dynamics assessment 
in order to achieve best efficiencies was performed. The assessment of each thermo- 
dynamic parameter was set in terms of non—dimensional charts, which facilitate the 
methodological approach to the design of a Tesla turbine. 
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The main considered parameters were: 
e Throat Mach number at stator outlet (and consequently mass flow rate); 
e Tangential velocity over rotational speed at rotor outer diameter (or rotor in- 
let). 


Throat Mach number 

Once the geometry assessment was carried out, increasing throat Mach number 
allowed an improvement of both efficiency and power. The increase in efficiency is 
moderate, whereas the power increase is quite relevant, as shown in Fig. 4.28. 
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Fig. 4. 28 Total to total efficiency and Power vs throat Mach number for various Tesla turbine 
dimensions 


Tangential velocity ratio 


The tangential velocity ratio (с = 10) is one of the most important parameters for 
2 


Tesla turbine optimization. The right matching of rotor inlet tangential velocity and 
peripheral speed is of paramount importance to achieve a high efficiency. In practise, 
the total to total efficiency is at its highest at o = 1, or very close to 1 (Fig. 4.29). This 
is due to the right value achieved by the inlet tangential relative velocity in this con- 
dition, which is very close to zero. At higher value of o, the fluid-machine work 
transfer would not be optimal, as the velocity would drop drastically at rotor inlet, 
dissipated into heat and not usefully transmitted to the rotor by the viscous forces. On 
the other hand, if a value lower than 1 is considered, a reversal flow conditions would 
be triggered. Indeed, if the absolute tangential velocity is lower than the rotational 
speed, a negative relative tangential velocity would be obtained at rotor inlet, so that 
the turbine would behave as a compressor at least in that region. Nonetheless, values 
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only little lower than 1, but close to 1 may be considered to achieve high efficiency 
levels. Indeed, if the flow reversal region is very limited, the higher power produced 
by the remaining inner region of the rotor, operating at a higher rotational speed, while 
keeping all other parameters unchanged, counterbalances the negative effect of the 
flow reversal. Therefore, best values of the tangential velocity ratio с were found in 
the range from 0.9 to 1. 
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Fig. 4. 29 Total to total efficiency vs. tangential velocity ratio for various Tesla turbine dimen- 
sions 


Compactness and rotational speed 

As discussed in the previous Section, the right matching between inlet tangential 
velocity in the rotor and rotational speed is of paramount importance to achieve high 
turbine performances. In order to have the proper match, the rotational speed needs to 
be adapted to the rotor external diameter. Specifically, the smaller the rotor, the higher 
the rotational speed required for best efficiency (Fig. 4.30). The machine compactness 
(i.e. power output per unit volume of the turbine) is another fundamental parameter, 
depending on the specific requirements of the field of application, and is clearly re- 
lated to the rotational speed. For example — referring to expanders in the power range 
from 1 to 30 kW — for the automotive sector, compactness is a fundamental parameter, 
and therefore a small, fast-turning Tesla turbines would be preferable; on the other 
hand, for domestic application, the Tesla concept offers — with respect to other possi- 
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ble expanders, such as centripetal turbines — the attractive possibility of direct cou- 
pling with a 3000 rpm generator and a low noise emission factor; for these applica- 
tions compactness of the machine could be sacrificed. 

The compactness factor CF is the ratio between the power and the total volume 
of the turbine, expressed in W/cm? (Dumont et al., 2018); the calculated values of CF 
are shown together with the rotational speed as a function of the rotor size in Fig. 4.30, 
where Too = 100°C, Poo = 8.33 bar and mass flow rates between 0.08 and 1.1 kg/s were 
considered. 
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Fig. 4. 30 Compactness factor and rotational speed for optimized Tesla turbine geometry (Ma, 
=1;0=1;0.04 <r, < 0.22) 


Tesla turbine geometry assessment 

The obtained geometry as a function of rotor outer diameter for each fluid showed 
the same range of design configuration; that is quite similar in terms of rotor outlet/in- 
let diameter ratio (0.3 < R< 0.4), throat width ratio (TW = 0.02), throat Mach number 
(Ma, = 1) and tangential velocity ratio (с = 1); nonetheless, the most important pa- 
rameter is the rotor channel width b, which is different for each fluid (Tab. 4.2). Re- 
ferring to ORC applications of the Tesla turbine, even if an optimal value of Б results 
is present for each working fluid, values around 0.1 mm are required in order to obtain 
high efficiency. Lower rotor channel width values are beneficial for refrigerants, and 
especially for fluids with low critical temperature and high critical pressure (such as 
R1234yf). Conversely, hydrocarbon fluids such has n-Hexane, allow higher rotor 
channel width, due to their critical properties, which are opposite of refrigerants. That 
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is, coupling a high critical temperature and low critical pressure. Furthermore, as dis- 
cussed in the previous Section, high Mach numbers (in the range of 1) are beneficial 
for turbine efficiency, as they improve the power production and the efficiency. None- 
theless, the tangential velocity ratio is the most important parameter when designing 
a Tesla turbine. High values of this parameter correspond inappropriate work transfer 
between the fluid and the rotor; on the other hand, values lower than 1 imply a reverse 
flow region at inlet. This feature was verified for all the examined working fluids. 


Table 4. 2 Calculated values of rotor channels width for the investigated fluids 


Fluid Rotor channel width [m] 
R1233zd(E) 0.0002 : D, +3-1079 
R245fa 0.00015 : D, + 3:107% 
R1234yf 0.0001 : D, + 2:107% 
n-Hexane 0.0003 - D + 5:107% 


When constrained rotational velocity applications are considered, high expander 
efficiencies are directly related to the accurate selection of the rotor diameter. Indeed, 
as shown in Fig. 4.31, fixing the rotational speed implies assuming an inlet tangential 
velocity ratio (at fixed thermodynamic conditions); it can be observed that efficiency 
is higher when the inlet tangential velocity ratio ø gets close to 1. Furthermore, it 
should also be noticed that the highest expansion ratio of the machine £ is achieved 
when o approached unity (keeping the fixed condition of unit Mach number at the 
throat section). 
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Fig. 4. 31 Tangential velocity ratio, efficiency and expansion ratio at a fixed rotational speed 
of 6000 rpm (Ma; = 1; 0.04 < r; < 0.22) 
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Joint analysis: Efficiency versus Power and Expansion Ratio 

The most suitable range of Tesla turbine design expansion ratio is between 3.5 
and 5.5, depending on the working fluid. The power range, which depends on the 
number of channels of the manufactured turbine, is between few Watts and 30-35 kW 
(considering configurations in a range 2 and 100 rotor channels). At high—pressure 
ratio, the turbine is subjected to large pressure losses, thus undergoing to an efficiency 
penalty, mainly due to the stator — rotor gap and to the high kinetic energy loss at 
expander output. Fig. 4.32 displays the total to total efficiency of a 100-channels 
Tesla turbine as a function of power and expansion ratio. It is important to notice that 
the maximum efficiency level is almost the same for all the considered fluids and lies 
between 0.609 and 0.626. Also the expansion ratio in order to achieve best efficiency 
values is similar to all fluids (between 3.5 and 5.5), but shows some sensitivity to the 
different fluid characteristics. With reference to the here considered fluids, 
R1233zd(E) and R245fa (very similar for thermodynamics properties) hold the same 
optimizing range of expansion ratio, i.e. between 4 and 5. On the other hand, R1234yf 
holds higher pressure at the same temperature level, and then shows optimal condi- 
tions between 3 and 4. Conversely, n-Hexane, achieving the lowest efficiency at the 
fixed 100?C temperature level, requires higher expansion ratios, between 4.5 and 6. 
Furthermore, best efficiency conditions are achieved at low power output, especially 
in the case of hydrocarbon fluids. 
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Fig. 4. 32 Efficiency as function of power and expansion ratio for a) R1233zd(E), b) R245fa, 
с) R1234yf; d) n-Hexane (Ma; = 1; 0.9 < ø < 2.25; 0.04 < т, < 0.22,n = 100) 
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Compactness and power (optimized speed) 

Fig. 4.33 displays the behaviour of the compactness factor (Dumont et al., 2018). It 
is a fundamental indicator when selecting a micro expander. As can be noted from 
Fig. 4.33, this Tesla turbine configuration can be quite bulky (CF <1) when large rotor 
diameters are considered, especially for fluids with low pressure levels such as n— 
Hexane (Fig 4.334). On the other hand, it can become quite compact (CF >1.5) when 
high power production is achieved through the utilization of refrigerants fluids with 
high pressure conditions (such as R1234yf, Fig 4.33c). Indeed, high pressure condi- 
tions mean that the mass flow rate elaborated by the turbine is higher and, conse- 
quently, the power output increases. Therefore, a compact turbine utilizing a refriger- 
ant fluid can be suggested for applications where the compactness of the expander is 
the most important parameter, such as in automotive. 
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Fig. 4. 33 Compactness (colour legend) and power (black lines on the graphs with correspond- 
ing caption) as function of turbine dimensions, channel and rotor diameter, for a) R1233zd(E), 
b) R245fa, c) R1234yf; d) n-Hexane (Ma; = 1; с = 1) 


Fig. 4.33 shows also the power range for each considered fluid as a function of 
rotor diameter and number of channels. The Mach number was fixed at 1 in order to 
achieve the maximum possible expansion ratio. As expected, higher rotor diameters 
and higher number of channels allow higher power production at the price of a bulkier 
expander. Anyway, the right compromise between compactness and power produc- 
tion depends on the selected fluid, but certainly rotor diameters between 0.16 and 0.24 
m guarantee a compact machine with reasonable power output levels. It should be 
remarked that a high-pressure condition, as in the case of R1234yf, allows achieving 
very high power and compact expander (almost 30 kW with a CF 71.3). Furthermore, 
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as was displayed in Fig. 4.30, compactness factor can be enhanced by selecting low 
rotor diameters and high rotational speeds. 

Comparing the obtained compactness factor with those of from volumetric ex- 
panders (Dumont et al., 2018), it seems that the Tesla turbine may be in the same range 
of roots and piston, and really close to scroll expanders. Nonetheless, it would cer- 
tainly be always bulkier than screw expanders, having a very high value of CF (up to 
20). 

Inlet total pressure and degree of superheating — Effects on Efficiency and Power 

Finally, a parametric investigation of the analysed Tesla turbines (working at Тоо 
= 100°C as thermodynamic design point) was carried out for efficiency and power 
production par channel. These fundamental performance parameters were evaluated 
as functions of total pressure and superheating level at turbine inlet. 

The 4 analysed fluids show the same trend of power output (Fig. 4.34). The power 
production gets higher at high-pressures and high superheating levels. On the other 
hand, the efficiency shows a different behaviour, depending on the working fluid. 
Specifically, for R1233zd(E), R245fa and n-Hexane, high efficiency conditions are 
achieved at low pressures and moderate super heating levels. Conversely, R1234yf 
presents optimised values for high pressures and moderate super heating level. This 
behaviour is mainly due to the very different thermodynamic properties of the fluids. 
R1234yf requires relatively large pressures in order to work with high temperature 
levels. This means that when pressure is low, also the considered temperature is quite 
low (for example, for the lower pressure case of 6 bar, with a 5 K super heating level, 
the total inlet temperature of the turbine is 25.5?C). With a reduction of the tempera- 
ture, the density increases and the velocities at stator outlet decrease; therefore, the 
rotor efficiency is reduced. High superheating levels are not really fundamental, as 
they produce higher velocities at rotor outlet, thus increasing the kinetic energy losses 
at the turbine discharge. 
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Fig. 4. 34 Efficiency and power as function of pressure and super heating level for a—b) 
R1233zd(E), с-а) R245fa, e-f) R1234yf; g-h) n-Hexane (Ma; = 1; ø = 1; 0.04 < r3 < 
0.22) 


153 


Micro turbo expander design for small scale ORC 


1.2.3 Comparison with Volumetric expanders 


Power range, operating conditions, efficiency are the main parameters that need 
to be taken into account when proper selecting an expander for micro—small ORC. 
Tab. 4.3 summarises the main characteristics of the principal volumetric expanders 
compared to Tesla turbines. It is important to note that the Tesla turbine seems to be 
a direct competitor of scroll, roots and piston expanders, while for higher power 
ranges, screw expanders will still probably be the reference. According to the results 
obtained by the model predictions (deeply discussed in the next Sections), the possible 
reachable isentropic efficiency for Tesla turbines seems to be around 62%, therefore 
considering a typical value for mechanical efficiency of axial turbines of 0.95 (Cohen 
et al., 1996), the overall expander efficiency would results in the order of 59% with an 
expansion ratio in the range between 3.5 and 5.5. Furthermore, it results that the 
proper power range for this expander extends from a few Watts (for very small rotor 
diameter configurations) to some tens of kW (for the largest diameter configurations). 
These evaluations rely on a proper balance of rotational speed, maximum power and 
temperature, rather than imposing absolute limits to one specific parameter. Rota- 
tional speed is a linear function of the rotor diameter, therefore the smaller the diam- 
eter the higher the required rotational speed for a reasonable efficiency. This favours 
mechanical design, as mechanical issues could arise for high rotor diameters, but not 
for small, compact disks. The Tesla turbine is claimed to be suitable for handling two— 
phase fluids (Rice, 1965). Indeed, in (Steidel and Weiss, 1976) a Tesla turbine for geo- 
thermal application was tested in two-phase flow conditions, nonetheless, further re- 
search on this topic is still required. 


Table 4. 3 Comparison of technical data between volumetric expanders and Tesla turbine 
(Rice, 1965, Steidel and Weiss, 1976, Dumont et al., 2018) 


Parameter Scroll Piston Screw Roots Tesla 
0.005 — 2,000— 1,000- 0.005 — 
Een tw] 10,000 9901710000 599.000 30,000 30,000 
Max Rotational 10,000 3000 (swashplate: 21,000 20,000 È 
Speed [rpm] 12,000) 
Built in volume ra- 1542 2.14 8 zi E 
tio 
Maximum Pressure -40 70 Е Е Е 
[bar] 
Maximum Tempera- 250 560 Е Е = 
ture [°C] 
ee Пон Yes Low Yes Yes Yes 
handling 
Isentropic efficiency 87 70 84 47 59 


The main advantage of a Tesla turbine would be its low cost, but in order to per- 
form a sound comparison with volumetric expanders this parameter cannot really be 
taken into account, as the cost depends not only on the complexity of the machine and 
on the manufacturing process, but is also linked to the maturity of the technology and 
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market development. Concentrating on technical aspects, Tab. 4.4 compares the main 
parameters, which are of paramount importance when selecting an expander; capabil- 
ity of operating at high pressure and temperature, compactness and efficiency. There 
is not an absolute best choice when selecting and expander for micro-ORC, as it 
deeply depends on the operating conditions of the specific application. Particularly, it 
seems that scroll expanders present higher efficiency values when compared to the 
other technologies, but still they are not suited to high pressure and temperature ap- 
plications, conversely to the piston expander. Screw, scroll and roots appear to be 
adaptable to wet expansion applications, while piston would be penalized. Screw ex- 
panders are the most flexible between the compared technologies, as they allow high 
efficiency in an extended power range. Furthermore, screw expanders are very com- 
pact. The Tesla turbine does not outstand in any of the analysed categories, but on the 
other hand it has still all positive “marks”, as it is a quite efficient, reliable and flexible 
expander, suitable for a wide variety of applications. 


Table 4. 4 Comparison of advantages and drawbacks between volumetric expanders and Tesla 
turbine 


High Pressure and tem- Wet Expan- Compact- Flexibil- Effi- 


perature sion ness ity ciency 
Piston ++ - > + + 
Screw — +++ HH 
Scroll = ++ TE 
Roots = + = = 
Tesla È + + d: + 


Based on a literature review (Quoilin et al., 2012, Lemort and Legros, 2017, Lemort et 
al., 2013, Lemort et al., 2012) and on the results provided in this research, it is possible 
to compare the Tesla turbine with volumetric machines and traditional bladed tur- 
bines. From this analysis, it appears that, even if the Tesla expander is intrinsically a 
turbine, it presents several similarities with volumetric expanders, which makes it a 
“hybrid” technology between the two main categories of expanders. In terms of power 
range, pressure and temperature levels, the Tesla turbine presents characteristics sim- 
ilar to volumetric expanders. However, the very moderate dependence of the pressure 
ratio to rotational speeds makes Tesla expanders similar to traditional bladed turbines. 
In conclusion, this type of turbine is an interesting option with unique characteristics, 
which makes it suitable and attractive for application in the field of low and micro 
power (<30 kW), pressure ratios between 2 and 6, where keeping low costs is manda- 
tory. Furthermore, it is remarked that, conversely to volumetric expanders, the Tesla 
turbine does not have any requirement for internal lubrication, which represents a 
great simplification of the whole circuit. 

Fig. 4.35 highlights the main characteristics of volumetric expanders, Tesla tur- 
bines and axial and radial turbines. Of particular interest is that Tesla turbine charac- 
teristics are closer to those of volumetric expander such as rotational speed, power, 
pressure and temperature operational ranges; with the only exception of the quasi— 
independency of pressure ratio with rotational speed, proper of dynamic machines. 
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Fig. 4. 35 Comparison of volumetric expanders, Tesla turbine and classical turbines 


2. Prototypes performance maps 


Performance maps of the real prototypes were developed. Particularly, the air pro- 
totype results are displayed as function of the main Turbomachinery parameters, while 
for the ORC, a deeper assessment on various fluids has been carried out. 


2.1 Air Tesla turbine 


Performance assessment 
The performance assessment was carried out first varying the Mach number at 
: : . : t 
nozzle throat outlet and the tangential velocity ratio at rotor inlet (о = =) at fixed 
2: 


thermodynamic conditions (total inlet pressure of 1.5 bar and total inlet temperature 
of 100°C); then the thermodynamic conditions were varied, fixing the Mach number 
at stator outlet at 0.4 and the tangential velocity ratio at 1. 

The tangential velocity ratio at rotor inlet is a fundamental parameter for optimi- 
zation of Tesla turbines. Particularly, as shown in Fig. 4.36 optimal efficiency is 
reached when o = 1 or even a bit lower than one. If ø is lower than one, at rotor inlet 
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there will be a region of reverse flow, which will therefore be a loss for the turbine. 
Nonetheless, values of o between 0.8 and 1 are still recommended, as the penalty by 
the reverse flow is counterbalanced by the higher power production at higher rota- 
tional speeds. 

Tesla turbine optimal efficiency is reached for low mass flow rates, and therefore 
or small values of the flow coefficient. Fig. 4.37 displays the @ — w diagram of the 
prototype. Particularly, the suggested range of operation would be: 0.05 < ф < 
0.08 and 0.6 < у < 1. The optimal specific speed and diameter are displayed in Fig. 
4.38. Fig. 4.38 was obtained changing both Mach number and sigma, for a fixed ge- 
ometry of the turbine, this means that the variations in specific diameter is due by the 
changing of both mass flow rate and rotational speed. The optimal operation range for 
this turbine prototype would be: 10 « d, « 12 and 0.16 « n, « 0.18. Fig. 4.39 as- 
sesses the reduced mass flow rate — expansion ratio curve at various reduced speeds. 
It is important to note that the shape of the curve, as expected, is the same of a “tradi- 
tional" turbine, with expansion ratio between 1 and 3.5 and reduced mass flow rate 
between 1-Е and 5E. 
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Fig. 4. 36 Total to total efficiency of turbine prototype against o and Ma; 
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Fig. 4. 37 ф — w diagram of air Tesla turbine prototype 
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Fig. 4. 38 n—d; diagram of air Tesla turbine prototype 
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Fig. 4. 39 т-а В diagram of air Tesla turbine prototype 

The thermodynamic analysis assessment results are displayed in Fig. 4.40 (а-а). 
It is important to note the contrasting behaviour of efficiency (4.40a) and power 
(4.40c). Tesla turbine optimal efficiency is reached for very low mass flow rates 
(4.40d), which are reached for low pressures and high temperature, if the Mach num- 
ber at stator outlet and the rotational velocity ratio are fixed. The developed prototype, 
provided the right conditions, could therefore reach a power as high as 2,200 W, with 
a rotational speed of over 24,000 rpm. 

Fig. 4.40b displays the compactness factor, as defined in (Dumont et al., 2018), 
which could be a fundamental parameter for the selection of a micro expander. The 
compactness factor is the ratio between the power and the total volume of the turbine, 
expressed in W/cm». As can be noted from Fig. 4.40d, this Tesla turbine prototype is 
quite bulky (CF «1) when high efficiency range is concerned, on the other it is quite 
compact (CF 71.5) when high power production is achieved. It seems therefore, that 
in order to design properly a Tesla turbine, it is necessary to keep in mind the specific 
application requirements. 
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Fig. 4. 40 a) түтт b) Compactness factor c) Power d) Mass flow rate of air Tesla turbine proto- 
type against total inlet pressure and temperature 


2.2 ORC Tesla turbine 


Before performing experimental investigation, it was important to assess the 
power and efficiency, which could be reached by the prototype depending on possible 
working fluid selection. This assessment was carried out on the geometry outlined in 
Section 3.3.2. 

The fluids here analysed are R1233zd(E), R245fa, R134a, R1234yf, SES36, n— 
Hexane and n-Pentane. R134a and R245fa are considered as references in literature, 
because as hydrocarbon substitutes they allow to achieve the best efficiency levels. 
R1233zd(E) and R1234yf are the cleaner fluids which are going to substitute R245fa 
and R134a, respectively. N-Hexane and n—Pentane are on the other hand adopted as 
reference hydrocarbons, due to their favourable thermodynamic critical conditions 
(particularly, by the low critical pressure when compared to refrigerants). Finally, 
SES36 is chosen, as it is an azeotropic mixture, behaving like a single substance, 
showing very convenient features for the Tesla turbine, as it allows lower pressure 
levels at higher temperatures when compared to R134a and R245fa. R404a was not 
taken into account in this analysis, as the temperature levels are not appropriate for 
this fluid, nonetheless a performance analysis on this fluid is presented at the end of 
this Section. 

As low-temperature resource applications are the most likely for this expander, 
two reference total inlet temperatures, Too = 100°C and Too =150°C, were considered. 
The total pressure at inlet is different for the fluids and it was selected, for both tem- 
peratures, in order to have superheated vapour 10?C above saturation temperature; the 
only exception are the cases for R134a and R1234yf at 150?C, which means super- 
heating about 50?C as the critical temperatures of the fluid are 101 and 94.7 ?C, re- 
spectively; a subcritical inlet pressure of 3.9 and 3.0 MPa were selected in these cases. 
The main parameters evaluated for all simulations are the power per channel, as well 
as the expander total to static efficiency. The calculations were performed varying the 
specific speed and the stator outlet Mach number, which is considered a fundamental 
parameter for the Tesla turbine (Guha and Smiley, 2009). 
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Low turbine entry temperature (100°C 

The power of the Tesla turbine, such as conventional turbines, increases as the 
flow coefficient Ф is increased. The increase in flow coefficient (and hence in mass 
flow rate) is reflected by the increase in nozzle throat Mach number. Higher throat 
Mach numbers correspond to large flow coefficients and high mass flow rates. The 
increase in the flow coefficient on one hand increases power, but on the other it de- 
creases the total to static efficiency (remembering that the geometry and the specific 
speed are fixed). 

An increase in specific speed ns allows obtaining a higher power production. 
Nonetheless, due to the characteristics of the machine, and especially of the conver- 
sion of power mechanism, the specific speed cannot reach exceedingly high values. 
This is caused by the increase of the absolute velocity at rotor exit, which after a cer- 
tain rotational speed would be higher than the value at rotor inlet, with considerable 
kinetic energy losses at exhaust. This behaviour is especially relevant for refrigerant 
fluids and in particular for R134a at high specific speeds. 

Higher efficiencies and lower power are reached with hydrocarbons when com- 
pared to refrigerants. This is mainly due to the different turbine inlet pressure. Typical 
values of turbine inlet pressure for R245fa, R1233ZD(E), R134a and R1234yf are 
higher than 1 MPa; on the other hand, the inlet pressure for n-Hexane and n-Pentane 
is much lower (0.18 MPa — 0.47 MPa). Lower pressures mean lower densities; there- 
fore, sonic conditions are reached at the nozzle throat for lower mass flow rates. As 
shown in (Manfrida and Talluri, 2018), low mass flow rates are beneficial for the turbine 
efficiency, but on the other hand are adverse for power production. 

SES36 seems attractive for the Tesla turbine, as it combines relatively low—pres- 
sure levels with comparably high temperatures. This feature combines the good traits 
of refrigerants (high power to density ratio) with those of the hydrocarbons (high ef- 
ficiencies), allowing the design of a high power density machine at desirable efficien- 
cies (~50%). 

Fig. 4.41 displays the rotor and stator efficiency separately. Of particular im- 
portance is the trend of stator efficiency, which is almost flat with values around 9095. 
This is possible thanks to the vaned configuration. The rotor efficiency, on the other 
hand, is deeply influenced by the flow condition and it is the main responsible for the 
trend of total to static efficiency. It has to be also remarked that the total pressure 
losses between stator and rotor play a fundamental contribution in lowering the ex- 
pander efficiency. Specifically, the total pressure loss in the turbine is in the order of 
20-30% of the overall pressure drop, with higher values connected to high nozzle 
Mach numbers and low specific speeds. Therefore, it is important to rightly design 
the dimensions of the stator—rotor interface. 

The main performance parameters of the assessed Tesla with various workings 
fluids for a turbine inlet temperature of 100?C, are resumed in Tab. 4.5. 
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Fig. 4. 41 Efficiency and power vs. stator outlet Mach number for low turbine entry tempera- 
ture case (Too =100°C) 
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Table 4. 5 Parametric analysis of selected fluids for total inlet temperature Too =100 [°C] 


R245fa — Poo = 1.009 [MPa] R134a P = 3247 [MPa] 
n; = 0.001 
Pa- 

Maj- Maj- =; SA Ма! Ma, Ma; = 
nee ij» qe MOR Mali o e. og Male 
Ma 013 022 030 0.36 025 047 067 084 
V 362 390 436 4.87 3.18 338 393 483 
n 0.06 0.07 0:09 0.11 0.06 0.07 010 0.13 


D; 38.7 35.42 32.21 29.52 38.03 33.78 28.74 24.19 
rpm 2393 3328 4021 4420 2357 3186 3629 3811 


p»/po 0.83 0.66 0.48 0.35 0.82 0.63 0.41 0.25 
ns = 0.002 
Ma» 0.35 0.49 0.59 0.65 0.58 0.80 0.99 Not 
Y 1.07 1.46 1.89 2.29 0.91 1.33 1.82 reached 
ф 0.02 0.03 0.04 0.06 0.02 0.04 0.05 


р, 40.01 33.91 29.22 26.10 37.57 30.35 24.74 chocked 
rpm 6772 7705 8284 8582 5963 6459 6658 rotor 


p/p | 0.60 — 0.5] 0.35 0.25 0.65 — 047 028 exit 
— RIz33zdE) Pa = 0.833 [MPa] __|__RI234yf Po = 3.080 [MPa] —— 
n= 0.001 
Pa- 

Ma\= Ma\= e = Ma, Ma, Ma; I 

p A NE M Lcx E ET A ж 
Маз 013 0.21 029 0.36 024 — 046 06 082 
V 366 3.95 4.42 5.05 32] 338 387 472 
6 006 0.07 0.09 0.12 0.06 007 010 0.13 


D; 38.14 35.41 32.22 28.96 38.10 34.00 29.11 24.62 
rpm 2460 3410 4110 4580 2080 2840 3300 3530 


p»/po 0.83 0.66 0.48 0.33 0.83 0.64 0.43 0.27 
n, = 0.002 
Ma» 0.34 0.47 0.57 0.64 0.57 0.78 0..96 Not 
Y 1.10 1.48 1.91 2.40 0.90 1.30 1.77 reached 
ф 0.02 0.03 0.04 0.06 0.02 0.04 0.05 


р, 40.02 33.99 29.31 25.58 37.82 30.74 25.19 chocked 
RPM 6950 7930 8520 8850 5330 5830 6110 rotor 


p»/po 0.69 0.50 0.35 0.23 0.67 0.48 0.31 exit 
— — — SES36-Ps-04396[MPa] —— —— 
n, — 0.001 
Pa- 
rame- is Maj;-0.6 Ma;-0.8 EE 
ter ` 
Maz 0.13 0.20 0.27 0.33 
Y 3.00 3.27 3.65 4.15 
ф 0.05 0.06 0.08 0.10 
р; 38.45 35.47 32.16 28.78 
трт 2656 3647 4383 4867 
рро — 0.81 0.64 046 031 
п; = 0.002 
Ма› 0.32 0.44 0.54 0.62 
bd 0.92 1.23 1.56 1.94 
ф 0.02 0.03 0.04 0.05 


р, 40.38 33.61 28.52 24.44 
КРМ 7420 8604 9320 9665 
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n—Hexane — Poo = 0.185 [MPa] n-Pentane — Poo = 0.469 [MPa] 
n, = 0.001 

Pa- 

Ma: Ма;= Er = Mai Mai Mai m 
ЭЖ Ж ctos cecus vob Шуу 

Ma» 0.09 0.14 0.19 0.27 0.10 0.16 0.22 0.28 
bd 3.54 3.92 4.38 5.00 3.68 4.01 4.46 5.05 
ф 0.06 0.07 0.09 0.11 0.06 0.07 0.09 0.12 


р, 38.83 35.94 32.64 28.60 38.33 35.70 32.70 29.44 
rpm 3498 4755 5707 6302 3491 4818 5831 6525 


P2/Po 0.81 0.64 0.46 0.28 0.83 0.66 0.49 0.33 
ns = 0.002 

Ma» 0.21 0.27 0.34 0.43 0.25 0.34 0.42 0.48 

Y 1.16 1:55 1.94 2.43 1.17 1.55 1.94 2.39 

$ 0.02 0.03 0.04 0.06 0.02 0.03 0.04 0.06 


D; 41.00 35.31 30.50 25.52 40.57 35.00 30.53 26.63 
m 9722 11255 12225 12480 9784 11413 12447 13000 


rp 
p2/po 0.69 0.51 0.34 0.20 0.70 0.53 0.37 0.24 
Medium turbine entry temperature (150°C) 


Higher temperature levels imply a larger power production. This is true also for 
the Tesla turbine. Indeed, with an increase of 50°C (as well as with an increase of the 
pressure levels) the power produced per channel reaches almost doubled values com- 
pared to the 100°C case. 

Nonetheless, if efficiency is considered, the trend is opposite. This behaviour is 
due to the corresponding increase of pressure. As stated above, the Tesla turbine 
reaches high efficiencies at low mass flow rate and with high nozzle throat velocities. 
For a fixed geometry and mass flow rate, high velocities are reached for low density 
values (that is, at low pressures and high temperatures); when using refrigerants as 
working fluids (as an alternative to hydrocarbons) at high temperature, it is necessary 
to work at relatively high pressure and this is the reason for the larger efficiency pen- 
alty which can be noticed comparing Figs. 4.41 and 4.42. At high temperature and 
pressure, higher rotor exit velocities are also reached, which again impairs the turbine 
efficiency. 

In this range of temperature, it seems that n-Hexane would be the most suitable 
fluid for the Tesla turbine, as it can guarantee relatively high efficiency, with quite 
high power production par channel. This feature is guaranteed by the much lower total 
inlet pressure of n-Hexane when compared to all other fluids. 

The performance parameters values of a Tesla turbine are in the same range of 
those of volumetric expanders or drag turbines, corresponding to relatively high spe- 
cific diameters and low specific speeds (Balje, 1981), as resumed in Tab. 4.6. Due to 
the nature of the turbine, a low flow coefficient is needed to determine good rotor 
efficiency (  — 0.01—0.1). On the other hand, the work coefficient ‘can be very high 
(over 3-4), especially at low specific speeds (ns = 0.001). Values of "between 1 and 
1.6 are recommended for higher specific speed (n; = 0.002). 
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Fig. 4. 42 Efficiency and power vs. stator outlet Mach number for medium turbine entry tem- 
perature case (Too =150°C) 
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Table 4. 6 Parametric analysis of selected fluids for total inlet temperature Too =150 [°C] 


Ro45ta— Pao = 2.815 [MPa] R13 Po = 3.975 [MPa] 
n, = 0.001 
Pa- 
Ma;= Ma: Gs E Mai Mai Mai Ax 
Bae Qu ue MBO cala | paio ‘pie ЖУО MAT 
Ma» 0.22 0.41 0.59 0.73 0.25 0.47 0.66 0.82 
bd 3.26 3.40 3.90 4.68 3.19 3.41 4.01 4.93 
ф 0.06 0.07 0.09 0.12 0.06 0.07 0.10 0.14 
D, 38.21 34.52 30.03 25.80 38.05 33.79 28.78 24.34 
rpm 2107 2960 3492 3822 2924 3914 4401 4567 
p»/po 0.84 0.65 0.45 0.28 0.80 0.58 0.38 0.23 
n, = 0.002 
Ma; 0.54 0.73 0.90 1 0.59 0.81 1.00 
bd 0.90 1.30 1.75 2.35 0.92 1.36 1.89 Not 
ф 0.02 0.03 0.05 0.07 0.02 0.04 0.05 reached — 
D, 38.44 3147 26.07 21.84 37.46 30.17 24.57 chocked 
rpm 5587 6182 6570 6757 7375 7890 8012 rotor exit 
p»/po 0.68 0.49 0.32 0.20 0.62 0.42 0.26 
RI233zd(E) — Poo = 2.290 [MPa] RT234yf Рос = 3.080 [MPa] 
n, = 0.001 
Pa- 
Ма;= Ma;= ay è Ma; Ma; Ma; С" 
Ж CE Л см М Ж. 
Ма» 0.2 0.36 0.53 0.66 0.21 0.37 0.53 0.64 
bd 3.37 3.56 4.02 4.T] 3.36 3.59 4.11 4.89 
ф 0.06 0.06 0.09 0.12 0.06 0.07 0.09 0.13 
D, 38.16 34.75 30.50 26.41 38.14 34.59 30.33 26.39 
rpm 2230 3100 3700 4050 2800 3820 4440 4760 
p»/po 0.84 0.66 0.46 0.29 0.81 0.62 0.42 0.27 
n, = 0.002 
Ma; 0.50 0.70 0.86 1.00 0.53 0.72 0.87 0.99 
bd 0.94 1.34 1.80 2.37 0.94 1.38 1.87 2.47 
ф 0.02 0.03 0.05 0.07 0.02 0.03 0.05 0.07 
D, 38.71 31.85 26.51 22.33 38.38 31.41 26.15 22.24 
rpm 6010 6640 7030 7210 7500 8100 8400 8450 
P2/Po 0.68 0.49 0.32 0.20 0.63 0.44 0.29 0.18 
SES36 Pu = 1402 [MPa] 
n, = 0.001 
Pa- 
rame- 3 Maj;-0.6 Ma;-0.8 м 
ter i 
Ma; 0.15 0.26 0.36 0.44 
bd 3.55 3.80 4.26 4.90 
ф 0.06 0.07 0.09 0.12 
D, 38.13 35.22 31.77 28.36 
rpm 1976 2757 3331 3726 
p/p. 0.84 0.67 049 0.33 
n, = 0.002 
Ма» 0.41 0.56 0.68 0.77 
bd 1.02 1.42 1.85 2.35 
ф 0.02 0.03 0.05 0.06 
D, 39.61 33.22 28.37 24.59 
m 5550 6240 6693 6990 
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n—Hexane — Poo = 0.601 [MPa] n—Pentane — Poo = 1.329 [MPa] 
ns = 0.001 
Pa- 

Ма;= Ма;= = с Mai Mai Mai o 
uae dodo vog. OPNS CDL ре жу Жо: 
Ma; 0.11 0.17 0.23 0.29 0.14 0.24 0.33 0.40 

Р 3.67 3.99 442 4.99 3.46 3.83 4.27 4.68 
ф 0.06 0.07 0.09 0.11 0.06 0.07 0.09 0.12 
D, 3832 3575 32.80 29.64 38.27 35.49 32.26 29.02 
rpm 3248 4524 5520 6240 3320 4539 5563 6285 
р2/ро _ 0.84 0.67 0.50 0.34 0.85 0.68 0.50 0.34 
n, = 0.002 
Ma 07 0.37 0.45 0.52 0.37 0.51 0.62 0.70 
Y 1.14 1.52 1.91 2.34 1.03 1.42 1.84 2:32 
ф 0.02 0.03 0.04 0.06 0.02 0.03 0.04 0.06 
D, 40.59 34.97 30.51 26.66 40.09 33.88 29.13 25.43 
rpm 9185 10745 11775 12480 9185 10505 11415 12025 
p2/po__ 0.71 0.54 0.38 0.25 0.70 0.52 0.36 0.24 


Performance assessment 


The performance assessment of the ORC prototype was carried out for two work- 
ing fluid: R404a and R1233zd(E). The first one is the fluid, which is available at the 
test bench of University of Florence, the second one, is on the other hand, the one 
available at the test bench of University of Liége. 

As for the air Tesla prototype performance assessment, first Mach number at noz- 


zle throat outlet and the tangential velocity ratio at rotor inlet (с = =) at fixed ther- 
2 


modynamic conditions (total inlet pressure of 12 bar and total inlet temperature of 
60°C for R404a and a total inlet pressure of 8.33 bar and a total inlet temperature of 
100°C for R1233zd(E)) were changed to see their influence on the efficiency and 
power of the turbine; then the thermodynamic conditions were varied, fixing the Mach 
number at stator outlet at 0.4 and the tangential velocity ratio at 1 for both fluids. 
R404a 

Fig. 4.43 displays the change in total to total efficiency of the Tesla turbine pro- 
totype when a total inlet pressure of 12 bar and a total inlet temperature of 60°C are 
considered. As discussed in the previous Sections, the optimal efficiency is reached 
when o = 1 or close to one. High efficiency is reached when relatively low Mach 
number and proper tangential velocity ratio are achieved. Increasing the tangential 
velocity ratio will bring to a very drastic drop in efficiency. 

As was the case for the air prototype, even for the ORC Tesla turbine prototype 
optimal efficiency is achieved through low mass flow rates, and therefore of small 
values of the flow coefficient. Fig. 4.44 displays the @ — у diagram of the ORC pro- 
totype (with R404a as working fluid for the above stated fixed thermodynamic con- 
ditions). Particularly, the suggested range of operation would be the same as the air 
prototype or: 0.05 < ф < 0.08 and 0.6 < у < 1. 

Total to total efficiency against specific speed and specific diameter is displayed 
in Fig. 4.45. The optimal working range of the ORC turbine prototype is for 21 « 
d, < 23 апа 0.07 < n, < 0.1, which are achieved through the right matching of 
mass flow rate and rotational speed. 
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Finally, Fig. 4.46 presents the reduced mass flow rate — expansion ratio curves at 
various reduced speeds. Expansion ratios up to 4 are achieved when reduced mass 
flow rate are close to 8-E°. 
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Fig. 4. 43 Total to total efficiency of ORC Tesla prototype utilizing R404a as working fluid as 
function of o and Ma; 
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Fig. 4. 44 ф — w diagram of ORC Tesla prototype utilizing R404a as working fluid 
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Fig. 4. 45 n, — d, diagram of ORC Tesla prototype utilizing R404a as working fluid 
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Fig. 4. 46 M,ia — В diagram of ORC Tesla prototype utilizing R404a as working fluid 
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The thermodynamic analysis assessment results are displayed in Fig. 4.47 (а-а). 
It is important to note the still contrasting behaviour between efficiency (Fig. 4.472) 
and power (Fig. 4.47c). Indeed, Tesla turbine optimal efficiency is obtained when 
very low mass flow rates are considered (Fig. 4.47d). Low mass flow rates are reached 
for low pressures and high temperature, if the Mach number at stator outlet and the 
rotational velocity ratio are fixed. The developed prototype, provided the right condi- 
tions, could therefore reach a power as high as 2,132 W, with a rotational speed of 
about 5,500 rpm. 

The ORC Tesla prototype is not a compact expander, as can be noted from Fig. 
4.47b, which displays the compactness factor. The compactness factor values for this 
specific prototype are very small. This was due to the initial design constraint of work- 
ing with a limited rotational speed. 


a) b) 


RI 
А 
A 
к 


LES 


м к 
yu x x 5585 9 


yx yx 55 9 


c) 


Ses 855 8 


ч н 
32 VR 5» E b BM 3 


хи n 


mis 


m з м sm „ и x м „ з м “ „ te n n m m 


w 1 з 
„ты = „ҥч = 


Fig. 4. 47 a) тт b) Compactness factor c) Power d) Mass flow rate of ORC Tesla turbine 
prototype utilizing R404a as working fluid, as function of total inlet pressure and temperature 


R1233d(E) 

The same trends that were obtained for the ORC Tesla prototype working with 
R404a are found when R1233zd(E) is used instead. The o — Ma, diagram (Fig. 4.48) 
provides the same information, stressing even more that a proper value of o is funda- 
mental, as high value of this coefficient will bring to a drastic drop in efficiency. 

Even the $ — w diagram (Fig. 4.49) indicates that the optimal range of operation 
would similar, just admitting slightly higher load coefficients: 0.05 < ф < 
0.08 and 0.6 « y « 1.2. 

The n;-d; diagram displayed in Fig. 4.50 suggests that the optimal operation range 
for the ORC turbine prototype when working with R1233zd(E) would be with 23 « 
d, « 25 and 0.08 « n, « 0.12. 
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To conclude, Fig. 4.51 assesses the reduced mass flow rate — expansion ratio curve 
at various reduced speeds. It is important to note that a little higher expansion ratio 
can be reached with this working fluid for a reduced mass flow rate of about 9E. 
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Fig. 4. 48 Total to total efficiency of ORC Tesla prototype utilizing R1233zd(E) as working 
fluid as function of o and Ma, 
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Fig. 4. 49 ф — y diagram of ORC Tesla prototype utilizing R1233zd(E) as working fluid 


171 


Micro turbo expander design for small scale ORC 
25 


Poo = 8.33 bar; Too = 100°C 
24 


22 0.5876 


ds 
N 
por 


0.5456 
20 0.5036 


0.4617 


0.4197 
0.3777 
18 0,3357 


0 0.04 0.08 0.12 0.16 0.2 


ns 


Fig. 4. 50 n, — d, diagram of ORC Tesla prototype utilizing R1233zd(E) as working fluid 
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Fig. 4. 51 m, — В diagram of ORC Tesla prototype utilizing R1233zd(E) as working fluid 


The thermodynamic analysis assessment results are displayed іп Fig. 4.52 (а-а). 
As for all the other cases, the opposing behaviour between efficiency (Fig. 4.52a) and 
power (Fig. 4.52c) is still present; directly linked to low mass flow rates trend (Fig. 
4.52d). The developed ORC prototype, provided the right conditions, could therefore 
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reach a power production as high as 1,800 W, with a rotational speed of about 5,150 
rpm. As analysed for the R404a, the compactness factor (Fig. 4.52b) is definitely low. 
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Fig. 4. 52 a) тт b) Compactness factor c) Power d) Mass flow rate of ORC Tesla turbine 
prototype utilizing R1233zd(E) as working fluid, as function of total inlet pressure and tem- 
perature 


3. CFD analyses 


3.1 Air Tesla turbine 


The aim of the computational assessment developed for the air Tesla prototype 
was to provide a benchmark on the fluid behaviour inside a Tesla turbine, in order to 
strengthen the results obtained with the in-house EES code. 

Particularly, in Fig. 4.53 (a-d) the results of static pressure, absoulte radial 
velocity, absolute tangential velocity and static temperature calculated with the EES 
model are compared with those computed with the CFD analyses carried out with a 
mesh made up of 3125000 nodes. The results obtained with the k-kl-omega model 
are very close to the ones achieved with the 2D home-built model. 
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Fig. 4. 53 Comparison between the results obtained with the EES and CFD analyses; a) static 
pressure, b) radial velocity, c) tangential velocity, d) static temperature, along the radial direc- 
tion 


3.2 ORC Tesla turbine 


Rotor 


Analysis with R404a working fluid 

Fig. 4.54 resumes the main computational results (mass weighted average values) 
achieved when R404a was taken into account as working fluid. The conditions, which 
were fixed at rotor inlet for all the analysed cases, are resumed as follows: 

e Static pressure = 1.16 MPa; 

e Static temperature = 60°C; 

e Mass flow rate = 0.004 kg/s; 

e Inlet angle = 88°. 


Three different rotational speeds (1500, 3000 and 4500 rpm) were analysed in 
order to assess its influence on the expander efficiency and power output, as well as 
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to understand the changes in the velocity profiles, pressure and temperature inside the 
rotor channels. 

The results obtained at 1500 rpm are the ones, which present highest discrepancy 
between the numerical models. Nonetheless, both the 2D EES in-house code and 
ANSYS Fluent predict clearly a very drastic reduction of tangential velocity at rotor 
inlet, which is the main responsible for power production. As expected, the radial 
velocity increases monotonically from periphery to centrum, as the wetted—wall area 
decreases and the mass flow rate is fixed by the continuity equation. 

A very close match between the results of the 2D EES model and those of the 3D 
CFD model were found. Specifically, a better match in the entrance region is reached 
when the EES model is compared to the 4—equation SST Langtry-Menter model, in 
respect of the laminar simulations. This is mainly due to the underestimation of the 
velocity drop by the laminar model, as it does not correctly predict the velocity profile 
in the entry region. Basically, the two different results (which anyway are really close 
to each other) differ due to the axial velocity distribution. 

The discrepancy of the results with rotational speeds of 3000 and 4500 rpm is 
almost non-existent in these cases. The reason for this is the lower difference between 
absolute tangential inlet velocity and peripheral speed. 

When comparing the behaviour of the flow at various rotational speeds, it is nec- 
essary to assess each parameter separately. For radial velocity only the curves at 3000 
rpm are presented, as the trend (i.e. increasing towards the centrum) is almost entirely 
coincident not depending on rotational speed, as it is determined by continuity equa- 
tion. As a consequence, the static pressure decreases monotonically from periphery to 
centrum. It must be remarked that higher rotational speeds imply higher pressure 
drops. The temperature drop is larger when pressure drop is larger, which is directly 
connected to the increment of rotational speed. It is interesting to notice the tempera- 
ture behaviour at low revolution per minute (Fig. 4.54d). In the first part (i.e. at rotor 
inlet), the fluid is partially heated, due to the abrupt reduction of tangential velocity 
in a very short distance. The turbine fails to convert all the available work, which is 
dissipatedly converted into heat. Nevertheless, the related temperature variation of the 
fluid is quite modest. 

Finally, the tangential velocity is the most interesting parameter to analyse. Ac- 
cording to Euler Equation, its change is responsible for work transfer to the rotor. The 
variation of the peripheral velocity of the expander deeply affects the tangential ve- 
locity behaviour. At low speed of revolution (1500 rpm) and fixed flow velocity at 
rotor inlet, the flow suddenly decelerates, not being able to transfer all the work to the 
turbine. At intermediate speed (3000 rpm), the decrease of tangential velocity is 
smoother and enables a more efficient work transfer to the rotor. Finally, at high speed 
of revolution (4500 rpm), there is a first part of the rotor where the tangential velocity 
increases, so that the turbine locally transfers work to the fluid as a compressor. There- 
fore, from the analysed behaviour at different speeds of revolution, it is clear that, 
when the nozzle conditions are fixed, an optimised speed allowing the maximum 
power extraction at the highest efficiency exists. 
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Fig. 4. 54 Radial behaviour at various rpm with R404a of: a) Tangential velocity; b) Radial 
velocity; c) Static pressure; d) Static temperature 


Tab. 4.7 summarises the results obtained both from 2D EES and 3D Fluent simu- 
lations. Specifically, it should be noticed that the efficiency is maximised in the con- 
dition where the relative tangential velocity at rotor inlet is close to zero, as can be 
noted in Fig. 4.54. Particularly, in Fig. 4.54a it can be noted that the for the 4500 rpm 
case absolute tangential velocity increases at inlet, corresponding to a negative rela- 
tive tangential velocity and therefore to a reduction of the expander efficiency. 


Table 4. 7 Results of the 2D EES model compared to the 3D Fluent results 


Rotor mass Static pressure 


Rotor ef- Power per . Static temperature 

ра ficiency channel [W] How rile at rotor mls at rotor inlet [°C] 
[kg/s] [Pa] 
2D EES model 
1500 0.55 2.80 0.004 1155344 59.42 
3000 0.69 5.45 0.004 1158488 59.47 
4500 0.65 7.91 0.004 1160474 59.42 
3D Fluent model 

1500 0.55 2.80 0.004 1155344 59.42 
3000 0.69 5.45 0.004 1158488 59.47 
4500 0.66 7.97 0.004 1160474 59.42 


Analysis with R134a, r245fa and R1233zd(E) working fluids 

Three further fluids were simulated in order to compare the in-house 2D code and 
the 3D CFD results. As shown in Figs. 4.55, 4.56 and 4.57, the results achieved with 
EES are in close agreement with those from Fluent, for both laminar and SST cases. 
The same flow behaviour of R404a was found for R134a, R245fa and R1233zd(E). 
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Tab. 4.8 summarizes the data of the simulations run for the three fluids at various 
rotational speeds flowing with the mass flow rate of 0.004 kg/s. As for the R404a for 
radial velocity, only the curves at 3000 rpm are presented. The three fluids present the 
same efficiency tendency when total inlet conditions are fixed. Particularly, the per- 
formance of R245fa slightly undergoes at high rotational speed (4500 rpm) compared 
to R134a due to the lower critical pressure that determines higher values of density at 
throat, thus lower velocity which results in a lower power output. Indeed, the lower 
values of velocity at stator outlet determine a more pronounced flow reverse at inlet 
for the highest analysed rotational speed (Figs. 4.56a and 4.57a), responsible for lower 
values of power and efficiency. 

As it is usual in ORC applications, the fixed variable is usually the high tempera- 
ture of the cycle (due to the heat source). Therefore, when a Tesla turbine is considered 
as an expander for ORC, it seems that the optimal fluid is the one, which guarantees 
lower densities at fixed temperature. Lower densities at fixed temperature are reached 
for lower pressure of the fluid; consequently, a suitable fluid for a Tesla turbine for 
an ORC application is the one, which hold high critical temperature and low critical 
pressure. 


Table 4. 8 Comparison of results for different fluids at various rotational speeds (from 1500 to 
4500 rpm), at the mass flow rate of 0.004 kg/s, total inlet pressure (1.2 MPa) and super heating 
level (37 °C) 


Power | Static tem- 
Static pres- 


с Rotor ef- per perature at Rotational 
Fluid nn sure at rotor i 
ficiency channel inlet [Pa] rotor inlet speed [rpm] 
[W] [°C] 
R404a 0.55 2.80 1155344 59.42 1500 
R134a 0.52 2.97 1159129 82.50 1500 
R245fa 0.56 2.78 1158426 133.50 1500 
R1233zd(E) 0.55 2.81 1160346 142.00 1500 
R404a 0.69 5.45 1158488 59.47 3000 
R134a 0.69 5.64 1160386 82.5 3000 
R245fa 0.69 5.25 1161838 133.50 3000 
R1233zd(E) 0.71 5.28 1156274 142.00 3000 
R404a 0.66 7.97 1160474 59.42 4500 
R134a 0.67 8.08 1154641 82.5 4500 
R245fa 0.64 7.44 1160462 133.5 4500 
R1233zd(E) 0.69 7.26 1114178 142 4500 
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Fig. 4. 55 Radial behaviour at various rpm with R134a of: a) Tangential velocity; b) Radial 
velocity; c) Static pressure; d) Static temperature 
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Fig. 4. 56 Radial behaviour at various rpm with R245fa of: a) Tangential velocity; b) Ra- 
dial velocity; c) Static pressure; d) Static temperature 
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Fig. 4. 57 Radial behaviour at various rpm with R1233zd(E) of: a) Tangential velocity; b) 
Radial velocity; c) Static pressure; d) Static temperature 


Coupled Stator—Rotor simulations 
The numerical analyses of the Tesla turbine with coupling of stator and rotor were 


performed with R1233zd(E) as working fluid. The boundary conditions values were 
imposed from the available data of the experimental campaign carried out at Univer- 
sité de Liége. Several thermodynamic points (resumed in Tab. 4.9) were simulated 
and the power and efficiency of the turbine computed and compared to the 2D in 
house code. 


Table 4. 9 Boundary condition for coupled stator-rotor analyses 


Total inlet pres- Total inlet temper- Static outlet pres- Static outlet tem- 
sure [Pa] ature [°C] sure [Pa] perature [°C] 
1500 473535 73.42 310895 68.39 
1750 476446 73.14 311083 67.82 
2000 479870 73.25 312114 67.67 
2250 482052 73.04 310820 67.53 
2500 486119 72.98 311721 67.35 
2750 489369 72.74 311641 67.01 
3000 493133 72.66 311778 66.72 


As it will be discussed deeply in Section 4.3.2, the simulation models do not take 
into account some fundamental sources of inefficiency, due to the real geometry of 
the turbine. Therefore, the obtained power from computational analyses cannot be 
compared with those from the experimental campaign, but the results are still of fun- 
damental value as they are representing the “ideal” maximum obtainable power if the 
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fluid does not incur in windage, partialization and pumping losses. Indeed, the simu- 
lated geometry does not take into account manufacturing issues, not correct alignment 
of stator and rotor channels, as well as secondary losses at the axial edges of the rotor. 

Fig. 4.58 displays the contours of static pressure, static temperature, tangential 
velocity and radial velocity for the 3000 rpm case. Comparing the coupled model re- 
sults with the rotor model results, it is very interesting to remark that partial admission 
effects are not negligible. Its effects can be easily highlighted when analysing the 
temperature trend in Fig. 4.58b. The temperature distribution inside the rotor is not 
uniform anymore, but it displays clearly four different temperature regions, which are 
due to the spiral trajectories of the fluid from the 4 nozzle admissions. Nonetheless, 
the difference in temperature between one stream and another is almost negligible, 
especially if power calculation is taken into account. The highest temperature drop 
occurs, as expected, at nozzle exit and the temperature drop in the rotor is very small, 
as is the pressure drop. Indeed, the pressure drop in the rotor is of about 30 kPa, which 
is a very small part of the total pressure drop of about 180 kPa (Fig. 4.58a). In this 
case, the Tesla turbine works essentially as an action turbine, where the pressure drop 
is almost entirely converted in velocity in the nozzles. Radial velocity maintains the 
same trend as the one assessed for the rotor model (i.e. increasing towards the cen- 
trum); as it is determined by the continuity equation. Although, close to nozzle exit, a 
peak of radial velocity is present. This peak does not globally influence the results of 
the simulations; nonetheless, it has to be taken into account that the number of nozzles 
affects the flow field of the turbine. 

The tangential velocity behaviour is close to the one analysed for the rotor only 
case. At first, the viscous forces effect allows a considerable reduction of tangential 
velocity, while as the flow approaches the inner radii, the angular momentum effect 
prevails, therefore increasing the velocity of the fluid. Furthermore, in this case the 
partial admission effect is not influencing the performance of the turbine in terms of 
power and efficiency when mass weighted average values are taken into account. 

Finally, a comparison in terms of power and efficiency is presented in Tab. 4.10 
between the computed numerical results and the 2D in house code results. It has to be 
remarked that the 2D code is now neglecting windage, pumping and partialization 
losses, which conversely will be taken into account when assessing experimental data. 
As can be noted from Tab. 4.10, the power prediction between 2D EES code and 3D 
computational fluid dynamics is really close, both in terms of power and efficiency. 
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Fig. 4. 58 Contours of a) static pressure b) static temperature с) tangential velocity d) radial 
velocity at 3000 rpm 


Table 4. 10 Results of the 2D EES model compared to the 3D Fluent results 


Tons Rotor effi- Turbine ef- Power per 
rpm mass flow . КИ Total power [W] 
tate [kës] ciency ficiency channel [W] 
2D EES model 

1500 0.232 0.32 0.14 4.57 274 
1750 0.234 0.36 0.17 5.38 322.7 
2000 0.235 0.40 0.19 6.18 371 
2250 0.237 0.43 0.22 7.03 421.8 
2500 0.239 0.46 0.24 7.87 472.2 
2750 0.241 0.49 0.27 8.71 522.8 
3000 0.242 0.51 0.29 9.57 574.3 
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3D Fluent model 


1500 0.232 0.30 0.15 4.55 273.2 
1750 0.234 0.34 0.17 5.36 321.6 
2000 0.235 0.38 0.19 6.16 369.6 
2250 0.237 0.41 0.22 7.00 420.1 
2500 0.239 0.44 0.24 7.84 470.1 
2750 0.241 0.47 0.26 8.67 520.2 
3000 0.242 0.50 0.28 9.52 571.3 


4. Experimental Campaigns 


4.1 Air Tesla turbine 


The experimental campaign was conducted exploring two different air mass flow 
rates, varying the rotational speed of the expander. The maximum rotational speed 
was limited by the brushless generator; the maximum achievable value was 3000 rpm. 
The thermodynamic conditions at turbine inlet for all tested points have been resumed 
in Tab. 4.11. 

The comparison between experimental data and numerical prediction are dis- 
played in Fig. 4.59 (shaft power) and 4.60 (total to total efficiency). The numerical 
predictions match very well the experimental test results for the investigated mass 
flow rate tested (0.028 and 0.030 kg/s). The limiting factor of the test bench was the 
admissible revolution per minute, as showed in Figs. 4.59 and 4.60, as this specific 
prototype could not be run at optimal speed. Indeed, Tesla turbine compactness is 
mainly dependent on rotational speed. If the turbine diameter requirements are strict 
(both for the application and mechanical issues), the expander needs to hold higher 
velocities in order to obtain high efficiency. On the other hand, a high radius design 
allows reaching high efficiency at relatively small rotational speed. The assessed pro- 
totype holds a relatively small radius, which means that optimal rotational speed 
would be in the order of some ten thousands rpm. At 3000 rpm, the maximum obtained 
net power was of 94 W. The thermodynamic power was assessed at 107 W, with a 
power loss due to bearings, generator and torque meter of about 13 W. 

The maximum assessed total to total efficiency of the turbine (at 3000 rpm) was 
11.3% when the lower mass flow rate tests were run. Indeed, as previously discussed, 
it was confirmed that efficiency is high at low mass flow rate. The mechanical effi- 
ciency of the turbine (comprehensive of torque meter and generator) was assessed as 
over 85%. 
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Table 4. 11 Thermodynamic conditions at turbine inlet 


rpm m = 0.030 [kg/s] m = 0.028 [kg/s] 
Too [°C] Poo [bar] Too [°C] Poo [bar] 
300 37.81 1.49 42.78 1.39 
600 40.59 1.48 43.82 1.39 
900 43.27 1.48 44.87 1.39 
1200 45.67 1.48 45.99 1.39 
1500 47.93 1.48 47.27 1.39 
1800 50.10 1.48 48.65 1.39 
2100 52.11 1.48 50.06 1.39 
2400 53.81 1.48 51.36 1.39 
2700 55.39 1.48 52.69 1.39 
3000 56.87 1.49 53.92 1.39 
100.0 ) 
90.0 e 
80.0 70227" » 
70.0 E ai 1 Bs a 
z 60.0 4. ee 
| p aai 2-7 
ы 50.0 = ае 
© ead p". 
È 40.0 n M 
30.0 2° 
2 
20.0 No ‘dei 
10.0 622°" 
0.0 
0 300 600 900 1200 1500 1800 2100 2400 2700 3000 
КРМ 
© Experimental Data; m = 0.030 [kg/s]  ----- EES; m = 0.030 [kg/s] 
© Experimental Data; m = 0.028 [kg/s] ^ ----- EES; m = 0.028 [kg/s] 


Fig. 4. 59 Experimental Data vs. numerical prediction, shaft power 
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Fig. 4. 60 Experimental data vs. numerical prediction, total to total efficiency 
4.2 ORC Tesla turbine 


R404a tests 

The ORC Tesla turbine was firstly arranged on the test bench of University of 
Florence described in Section 3.3.2. The electric motor, the torque meter as well as 
the servo drive were the same as the one described for the air Tesla turbine tests in the 
Section 3.5.1. The test bench was developed within the project Exp-HEAT, which was 
funded by the EU (grant no. 605923) under the 7th Framework Program (FP7 Capac- 
ities). In Exp-HEAT project, a piston expander was experimentally tested; this ex- 
pander required a very high quantitative of oil in order to achieve high efficiencies, 
therefore an auxiliary lubrication system was employed. This resulted in a lot of oil 
remaining in the test bench circuit, which unfortunately hindered the test on the Tesla 
expander. Indeed, as can be seen in Fig. 4.61, the turbine worked as a “honey dipper”, 
where the oil was trapped in the gap between stator and rotor, giving rise to very high 
pumping losses. Various pressures and mass flow rates were tested, but none was able 
to achieve a stable condition, as the hindering of the oil was high. Furthermore, at first 
it was questioned if the turbine was correctly working, therefore it was arranged on 
the same test bench of the air Tesla turbine and tested with the same inlet condition 
Poo — 1.5 bar, m — 0.030 kg/s and various rotational speeds (from 300 to 3500 rpm). 
The experimental power produced by the turbine well matched the 2D EES model 
prediction, as shown in Fig. 4.62. 
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Fig. 4. 61 ORC tesla turbine jammed with oil 
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Fig. 4. 62 ORC tesla turbine jammed with oil 


185 


Micro turbo expander design for small scale ORC 


R1233zd(E) tests 


The test bench of Université de Liége allowed for the adjustment of: 


e The pump mass flow rate (through a frequency inverter), which therefore per- 


mitted the imposition of the mass flow rate of the test bench; 


e The heat input (changing the heater settings of temperature and air mass flow 


rate), which allowed the regulation of super heating level; 


e Тһе heat output (changing condenser water mass flow rate), which provided 


a control on the lower pressure of the test bench; 


e The rotational speed of the turbine (through a frequency inverter), which al- 


lowed to change the pressure drop in the turbine. 


Therefore, the experimental campaign was conducted exploring different refrig- 
erant mass flow rates (0.25 — 0.36 kg/s), varying the rotational speed of the expander 
(1000 — 5000 rpm), as well as total inlet pressure (4.7 — 6.7 bar) and super heating 
levels (3.2 — 46 ?C). In order to have confident temperature measurements, the turbine 
and the inlet and outlet pipes were insulated. The explored thermodynamic conditions 
at turbine inlet for all tested points are resumed in Tab. 4.12. 


Table 4. 12 Thermodynamic conditions at turbine inlet 


; Rotational Speed Mass flow rate Too SH Poo 
Nomenclature Points [rpm] [kg/s] [°С] [°С] [Ра] 
1 2000 0.3568 86.84 8.23 635633 
2 2500 0.3566 87.13 7.94 644701 
D1 3 3000 0.3564 87.08 7.36 653135 
4 3500 0.3565 87.19 6.79 664151 
5 4000 0.3567 87.11 6.11 673920 
6 1500 0.2541 73.42 6.34 473535 
7 1750 0.2541 73.14 5.83 476446 
8 2000 0.2541 73.25 5.67 479870 
D2 9 2250 0.2540 73.04 5.29 482052 
10 2500 0.2538 72.98 4.91 486119 
11 2750 0.2539 7274 442 489369 
12 3000 0.2539 72.66 4.04 493133 
13 3500 0.2532 108.10 37.52 518962 
14 3250 0.2531 108.56 38.26 515830 
15 3000 0.2530 108.39 38.32 512314 
m 16 2750 0.2530 108.32 38.48 508854 
17 2500 0.2529 107.60 38.00 505954 
18 2250 0.2530 107.01 37.76 501261 


186 


Lorenzo Talluri 


19 1500 0.3009 76.38 11.33 518794 
20 1750 0.3009 76.32 11.21 521596 
21 2000 0.3009 74.98 9.73 523252 
22 2250 0.3009 75.31 10.00 527202 
D4 23 2500 0.3008 74.84 9.55 530978 
24 2750 0.3007 75.44 10.06 536251 
25 3000 0.3008 75.00 9.19 539307 
26 3250 0.3007 74.96 9.28 543242 
27 3500 0.3006 75.42 9.43 547545 
28 3750 0.3008 74.59 8.22 550591 
29 5000 0.2993 120.26 44.09 598959 
30 4500 0.2994 120.13 44.60 588417 
DS 31 4000 0.2993 120.20 45.45 577387 
32 3500 0.2992 119.92 45.73 569026 
33 3000 0.2992 119.76 46.14 561617 
34 1000 0.3651 78.59 3.73 578933 
35 1500 0.3646 79.28 4.06 584299 
36 2000 0.3642 79.88 4.32 589228 
D6 37 2500 0.3640 79.93 3.90 596259 
38 3000 0.3641 79.81 3.13 605988 
39 3500 0.3635 82.08 4.65 616757 
40 4000 0.3640 81.26 3.35 624814 


Fig. 4.63 resumes the thermodynamic conditions tested as function of turbine ex- 
pansion ratio. The explored range of expansion ratio is not very wide, but it still allows 
a very significant data analysis. 

As can be depicted from Fig. 4.63 (a) and (d), higher mass flow rates are linked 
to higher total inlet pressure, as well as at higher expansion ratio. Nonetheless, higher 
expansion ratios are also obtained increasing the super heating (Fig. 4.63 (c)). The 
expansion ratio is indeed directly linked to the thermodynamic condition at nozzle 
throat. High inlet pressure implies higher mass flow rate, and therefore higher pressure 
drop. Moreover, higher super heating level, at same mass flow rate (data D1 and D3 
or D4 and DS) is connected to a higher velocity at throat section and therefore to a 
higher pressure drop. High temperature implies lower densities if all other conditions 
are fixed. Therefore, due to mass balance, higher velocity is reached at throat and 
consequently at rotor inlet. 

As expected, increasing rotational speed results in an increase of expansion ratio, 
as displayed in Fig. 4.63 (d). The slope of the curves is almost the same for all inves- 
tigated conditions. The highest obtained expansion ratio (1.87) was obtained with a 
mass flow rate of 0.36 kg/s with a super heating level of 3.35 °C at 4000 rpm. 
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Fig. 4. 63 Experimental data: mass flow rate a), rotational speed b), super heating level c) and 
total turbine inlet pressure d) as function of expansion ratio 


The maximum achieved experimental thermodynamic power was 906 W, with a 
mass flow rate of 0.299 kg/s, 44 °C super heating level and 5000 rpm rotational speed. 
The highest power conditions were indeed obtained for high super heating levels (data 
D3 and D5), followed by high mass flow rate conditions (data D1 and D6). Further- 
more, as expected, thermodynamic power increases as the expansion ratio augments. 
Fig. 4.64 displays the behaviour of thermodynamic power as function of the expansion 
ratio. 

It is very interesting to notice that the same expansion ratio can be reached either 
with a high super heating level or with a higher mass flow rate, but the power produc- 
tion is higher in the former case. Indeed, if we take as an example, an expansion ratio 
of about 1.8, it can be seen from Fig. 4.64 that the turbine produces nearly 18% more 
power with a mass flow rate of 0.299 kg/s and a super heating level of 44.6 °C com- 
pared to the case with 0.364 kg/s of mass flow rate and a 3.2 °C super heating level. 
This is due to the better match between rotational speed and tangential velocity, al- 
lowing therefore a better power conversion. 

Fig. 4.65 displays the thermodynamic power as function of expansion ratio and 
rotational speed. High expansion ratios and rotational speeds favour power produc- 
tion, while low rotational speeds and expansion ratios hinder the power production. 
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. 4. 64 Experimental thermodynamic power output vs. expansion ratio 
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. 4. 65 Experimental interpolated thermodynamic power as function of expansion ratio and 


rotational speed 
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Differently from experimental thermodynamic power, the maximum achieved 
shaft power output was 371 W for two thermodynamic states. In the first case, it was 
reached with a mass flow rate of 0.299 kg/s, a super heating level of 45.5 °C and a 
rotational speed of 4000 rpm; the same power output was also achieved for a mass 
flow rate of 0.365 kg/s, a super heating level of 4.65 °C and a rotational speed of 3500 
rpm. 
As can be easily noticed comparing Figs. 4.64 and 4.66, the shaft power output 
shows a different behaviour compared to the thermodynamic one. Particularly, the 
behaviour of the thermodynamic power output is monotonically increasing with ex- 
pansion ratio (and rotational speed), while the shaft power output presents a maximi- 
zation point. 

The trend of shaft power can be better understood when analysing Fig. 4.67. In- 
deed, in Fig. 4.67 the shaft power is displayed as a function of expansion ratio and 
rotational speed. Expansion ratio increases when augmenting the rotational speed of 
the turbine, but it also directly raises the mechanical losses due to the bearings, as well 
as the friction losses due to the electromagnetic coupling. 

On average, a 50% organic efficiency was achieved that is a really low value, 
mostly attributable to the improper alignment of the magnetic coupling, which is re- 
sponsible for a really high increase of the mechanical losses, due to the friction. 
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Fig. 4. 66 Experimental shaft power as function of expansion ratio 
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Fig. 4. 67 Experimental interpolated shaft power as function of expansion ratio and rotational 
speed 


The maximum achieved experimental thermodynamic total to total efficiency was 
30% with a mass flow rate of 0.253 kg/s, a super heating level of 37.5 °C and a rota- 
tional speed of 3500 rpm. 

Highest efficiency values are directly related to the thermodynamic power pro- 
duction; nonetheless, as expected for a Tesla turbine, higher efficiencies are achieved 
at low mass flow rates. Indeed, Figs. 4.64 and 4.68 are examined together. It clearly 
appears that there is a direct relationship between power production and efficiency. 
However, the low mass flow rates conditions (data groups D2 and D3) achieve higher 
efficiencies compared to the high mass flow rates (data groups D1 and D6). On the 
whole, an average 17% thermodynamic efficiency was obtained for this expander, 
with higher values for low mass flow rate conditions. 

The shaft efficiency (Fig. 4.69), on the other hand, is directly related to the ob- 
tained shaft power. Therefore, the highest efficiency condition is achieved at the max- 
imum power output condition, which is at a mass flow rate of 0.299 kg/s, a super 
heating level of 45.5 °C and a rotational speed of 4000 rpm; and the achieved value 
of efficiency is of 9.62%. Nonetheless, there is still an influence of the mass flow rate; 
at low mass flow rates, the efficiency still is still relatively high, even at lower power 
production. 

On the whole, an average shaft efficiency of 8.2% was obtained for this expander, 
with higher values at low mass flow rate conditions and higher power production con- 
ditions. 
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Fig. 4. 68 Experimental thermodynamic efficiency of the Tesla turbine prototype vs. expansion 
ratio 
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Fig. 4. 69 Experimental shaft efficiency as function of expansion ratio 
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As can be noted from the very big difference between thermodynamic and shaft 
efficiency, as well as from the values of organic efficiency, the mechanical power 
losses of the turbine are quite high. Fig. 4.70 resumes the experimental power losses 
of the turbine, compared to the predicted bearing losses and the sum of bearings losses 
and friction losses due to the contact of the electromagnetic coupling. For the estima- 
tion of the bearings mechanical losses, the model provided by the manufacturer was 
applied (SKF, 20184). The estimated bearing losses are expressed by Eq. (4.8). 


Pross = €x Mies (4.8) 


Where Miss is the bearing frictional torque, which takes into account the rolling 
and sliding contributions 

For the estimation of the losses derived by the contact of the electromagnetic cou- 
pling, the model of friction losses on a carrier pin was applied. The estimated pin 
friction losses are expressed by Eq. (4.9). 


Poss,pin = wM, (4.9) 


Where M, is the frictional torque, which is the friction force time the radius of 
electromagnetic coupling, with a coefficient f= 0.15 for steel over steel materials. 
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Fig. 4. 70 Mechanical losses of ORC Tesla turbine 


Comparison with 2D EES code 

As anticipated in Section 4.2.2, the 2D in-house EES code did not initially con- 
sider partialization, windage and pumping losses. However, if these are not taken into 
account, the power predicted by the model is much higher, than the experimental 
value. Therefore, in order to improve the reliability of the Tesla turbine calculation 
code, a model including the effects of partial admission, windage and pumping losses 
was added (Traupel, 1977). 

The estimated windage losses are expressed by Eq. (4.10) 

p, = 0.1 THE. p. цз 


эр (4.10) 
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Where Н is the total thickness of the rotor disks (0.0008*30) and e is the partial- 
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ization degree defined as = = 1 — 
Ainrotor 


The estimated partialization losses are expressed by Eq. (4.11) 


= Vas үрүү. б) ut 
Pot = 0.15 .m à, z (4.11) 
Where 74, is the isentropic absolute velocity at stator outlet. 

The estimated pumping losses are expressed by Eq. (4.12) 


Pop = 4: Cy: p: 02:03 (4.12) 


Where Cy is a coefficient function of Reynolds number, Cy = 0.003 · Re^? 
(Traupel, 1977). 

This model was developed for partial admission steam turbine, and therefore it 
takes into account experimental coefficients (0.1 for windage losses and 0.15 for par- 
tialization losses), which do not proper match for an organic fluid. 

Fig. 4.71 displays the effect of partial admission, windage and pumping losses in 
the thermodynamic diagram. Particularly, these losses increase the temperature at ro- 
tor inlet, besides increasing the entropy, which results in a reduction of the available 
work. 
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Fig. 4. 71 Thermodynamic representation of losses incurring in stator—rotor gap 


It was then decided to try interpolating the ratio between actual experimental 
power and the power calculated by the EES 2D model, including partial admission, 
windage and pumping losses. The experimental fit coefficient was interpolated as a 
function of Mach number and temperature at the stator outlet, after having included 
the enlargement losses. 

The interpolated equation is expressed by Eq. (4.13). In spite of the not very high 
amount of available experimental data, which do not allow an accurate experimental 


194 


Lorenzo Talluri 


fitting of the whole dataset, the procedure and the interpolated coefficient still pro- 
duces a reliable prediction. Fig. 4.72 shows the expression of the experimental inter- 
polation surface, as a function of Mach number and temperature. 


Сехр = —31.08 — 0.023 — 31.08 T + 100.8 * Ma + 0.0199 - T: Ma 
— 77.04. Ma? 


Fig. 4. 


(4.13) 


72 Interpolation surface of experimental constant 


Finally, the modified 2D EES in-house code, including the above models to ac- 
count for the partial admission, windage and pumping losses, with the tuned experi- 
mental coefficient was applied and the results are displayed in Figs. 4.73 and 4.74. 

In this way, the 2D EES code, properly describes the behaviour of both thermo- 
dynamic power output and efficiency, with a satisfactory agreement level between 
calculated and measured data. 
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73 Experimental data and numerical prediction of thermodynamic power vs. expansion 
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Fig. 4. 74 Experimental data and numerical prediction of thermodynamic efficiency vs. expan- 
sion ratio 


Comparison with 2D EES code and CFD 
To conclude, 2D model (both including and neglecting partialization, windage 


and pumping losses), 3D CFD and experimental data (D2) are compared and dis- 
played in Tab. 4.13. The 2D model results not including losses are close to the one 
obtained from 3D CFD simulations, but they are not well matching the experimental 
data. As discussed in the previous Section, this is due to the ideal conditions of the 
simulations, which assume that all the mass flow rate which exits from the stator en- 
ters the rotor. Unfortunately, it was found experimentally that part of the mass flow 
rate does not enter the rotor, therefore, partialization, windage and pumping losses 
cannot be neglected. 

Finally, Tab. 4.14 shows the relative deviations between the results predicted by 
the EES 2D code, the SST-CFD code and the experimental data (D2). The relative 
error is defined as the ratio between the absolute value of the local difference between 
the experimental data and the numerical results, and the value of the experimental data 
(the relative error is referred to the experimental data). It is interesting to notice that 
the introduction of partialization, windage and pumping losses, allows to strongly de- 
crease the relative deviation between numerical and experimental data. 
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Table 4. 13 Results of 2D EES model (with and without partialization, windage and pumping 
losses) compared to 3D Fluent and experimental data (Dataset D2) 


Efficiency Total power [W] 
2D 2D 
р PES 3D Experi- сы EES 3D Experi- 
rpm EES model EES model 
; CFD mental ‘ CFD mental 
model with model with 
losses losses 
1500 0.14 0.14 0.15 0.15 274 298 273 304 
1750 0.17 0.16 0.17 0.17 323 331 322 353 
2000 0.19 0.17 0.19 0.19 371 357 370 400 
2250 0.22 0.17 0.22 0.17 422 371 420 363 
2500 0.24 0.17 0.24 0.17 472 382 470 372 
2750 | 027 0.17 0.26 0.17 523 368 520 374 
3000 0.29 0.15 0.28 0.18 574 339 571 401 


Table 4. 14 Relative deviation between experimental data and numerical models 


Efficiency Total power [W] 
2D EES 2D EES 
rpm eu model with 3р CFD Li model with 3р CFD 
losses losses 
1500 4.2% 2.1% 2.7% 9.9% 2.0% 10.1% 
1750 1.4% 6.1% 1.4% 8.6% 6.2% 8.9% 
2000 0.7% 10.7% 0.7% 7.3% 10.8% 7.6% 
2250 31.0% 2.2% 31.0% 16.2% 2.2% 15.7% 
2500 40.9% 2.6% 40.9% 26.9% 2.7% 26.4% 
2750 60.0% 1.4% 54.0% 39.8% 1.6% 39.1% 
3000 62.6% 15.5% 57.0% 43.296 15.596 42.596 
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Chapter 5 
Conclusions and Recommendations 


The present work collects a wide part of research activities performed during the 
Ph.D. course, dealing with the design of a “new—old” expander technology, for small 
scale organic Rankine cycle applications. Conclusions of the developed research are 
now drawn, as well as recommendations for future works, which could be a spark for 
researchers to move forward in a field, which still presents very engaging topics to be 
unveiled. 


1. Conclusions 


Over the last years, energy systems research grown a strong attention towards the 
study and development of proper solutions for small and micro distributed systems 
for heat and power generation, to be applied in domestic buildings or industrial facil- 
ities. The utilization of ORCs is a sound solution in this field, especially when the 
cycle is made of efficient, reliable and low cost components. A critical part of micro 
ORCs is the expander, as it often does not combine low cost and reliability require- 
ments. In this domain, Tesla turbine seems to be a valuable candidate to tackle these 
issues, as its simple structure guarantees a low cost, reliable and quite efficient ex- 
pander. 

Open literature review showed that Tesla turbine has gained a renewed interest in 
the last few years, mainly due to the higher attention that micro power generation 
gained on the energy market. Four principal lines of research stirred up the interest on 
the assessment of Tesla turbine. The first is the one carried forward by W. Rice, which 
developed a throughout assessment from analytical to experimental campaigns on air 
Tesla turbines. The second is the one conducted by A. Guha and S. Sengupta, which 
realized both analytical and computational assessment of the Tesla turbine, deeply 
assessing the role of each force inside the rotor, as well as conducting several inter- 
esting analyses to make clear the causes of inefficiency of a Tesla turbine. The third 
line of research is the one proposed by VP. Carey and his team, which, improved the 
model firstly realized by W. Rice and applied it to Watt and sub—Watt applications. 
The last line of research is the one brought forward by C. Schosser and M. Pfizner, 
which numerically and experimentally investigated the flow path behaviour inside an 
air Tesla turbine. Outside these four main lines of research, very interesting spot work 
were developed, such as those of P. Lampart or J. Song, which started to analytically 
and numerically analyse the Tesla turbine for ORC applications. 
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As highlighted by the state of the art, several analytical and numerical models 
were realized, and many experimental campaigns were carried out, nevertheless, the 
majority of these investigations were developed on Tesla turbines utilizing air as 
working fluid, while, in the present study, the analysis is carried out from a slightly 
different perspective, which is setting the focus on Tesla turbine for ORC applications. 
For this reason, a clear and complete design and optimization procedure for ORC 
Tesla turbines with a model including real gas equations and concentrated pressure 
losses was developed. Furthermore, a straightforward methodology for geometry as- 
sessment as well as computational analyses on the flow behaviour inside the disks of 
a ORC Tesla turbine, and experimental investigations utilizing air and R1233zd(E) as 
working fluids, were realized. 

In the first part of the work a 2D in-house code was developed in Engineering 
Equation Software ambient. Initially, it was tested and validated against available lit- 
erature data utilizing air as working fluid. After the validation, the model was ex- 
tended to organic working fluids and a full design and optimization procedure for 
ORC Tesla turbines was realized. A pivotal point of this first part of research was the 
innovative design procedure of the expander compared to the previous literature: it 
did not only take into account the rotor or the stator separately, but it also included all 
the assembled components together, from the plenum chamber to the diffuser, passing 
through the losses in the gap between the stator and the rotor. Furthermore, a revised 
conceptual modular design of the turbine was proposed and the performance of sev- 
eral working fluids with specific prototype size assessed. Finally, a scheme for geom- 
etry assessment was also featured in order to perform a qualitative comparison with 
small-size expanders, which are the direct competitor of the Tesla turbine. 

The key outcomes of this part of the work may be summarised as follows: 

e Anovel methodology for the complete design ofa Tesla turbine working both 
with air and organic fluids was proposed and assessed. Each component was 
designed taking into account the mutual relationships between the different 
parts of the expander. 

e An innovative model for the solution of the rotor flow field was developed. 
Starting from an existing literature approach, the new one was generalized, 
considering real fluid behaviour and influence of Reynolds number on the 
velocity profile inside the rotor channels. 

e A sensitivity analysis to each geometric and thermo-fluid dynamic parameter 
was carried out. It was found that performance, mass flow rate and expansion 
ratio are strictly connected: low mass flow rates allow better efficiency and 
lower power output. Low mass flow rates, for a fixed geometry of the nozzle 
and fixed velocity at the throat, are obtained for low density at nozzle exit 
(from continuity equation); therefore, high temperatures and low pressures 
are necessary for a proper design of the Tesla turbine rotor. The right choice 
of channel height and in/out rotor diameter ratio are of primary importance 
in optimization of the expander efficiency. 

e An improved design concept was introduced simplifying the shaft/rotor as- 
sembly with a modular, robust construction principle and possibility of sealed 
operation. 
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e Due to its intrinsic working principle — work transmitted by friction — the 
Tesla turbine results to be competitive with conventional expanders only for 
low-power application, but it is not a good candidate for medium to high 
power as the several losses involved, such as the high kinetic energy at ex- 
haust, the rotor efficiency drop due to the increase of mass flow rate, and the 
higher pressure drop into the stator/rotor gap. 

e The numerical assessment indicated that the Tesla turbine appears potentially 
competitive with other expanders for low n; (0.001—0.005) and high а, (20— 
50) (typical range for volumetric expanders or drag turbines) with special 
reference to efficiency. The rotational speed has a strong influence on the 
expander power and efficiency, but generally, the turbine can be sized to 
work properly within 4000—8000 rpm. 

e The most critical parameters for achieving good turbine performance were 
found to be the rotor inlet tangential velocity ratio, the stator throat Mach 
number, the rotor channel width and the rotor outlet/inlet diameter ratio. From 
the analyses of several computations on different working fluids, it was found 
that, for all fluids the throat Mach number and the rotor inlet tangential ve- 
locity ratio should be close to 1, in order to achieve high efficiencies. Fur- 
thermore, the proper value for rotor channel width was found to be a linear 
function of rotor inlet diameter and optimal values for rotor outlet/inlet diam- 
eter ratio were found to be between 0.3 and 0.4. 

e Suitable design expansion ratios for the Tesla turbine were determined be- 
tween 3.5 and 5.5. This range of expansion ratios is quite common in low 
temperature applications, which may be considered, therefore, to be the opti- 
mal field of application of this turbine. 

e Two prototypes (one utilizing air as working fluid, the other organic fluids) 
were finally designed, taking into account thermodynamic considerations, 
mechanical restriction, as well as test bench operating ranges. 


Once the fundamental principles of the turbine were unveiled, and the design of 
two prototypes was carried out, a computational assessment on the performance of 
the Tesla turbine rotor, as well as on the stator rotor interaction, was performed. The 
analysis was conducted through the evaluation of the flow behaviour and the conse- 
quent performance parameters of the turbine, whose size was referred to the actual 
built prototype. A pivotal point of this part of the work was the comparison of the 
results achieved by the 3D CFD model realized with the commercial code ANSYS 
Fluent and the 2D home built EES model. The realization of the stator rotor interaction 
simulations allowed understanding that partial admission plays a fundamental role in 
the behaviour of the fluid inside a Tesla turbine rotor. 

The key achieved results of this part of research may be summarized in the fol- 
lowing: 

e Тһе results of the rotor simulations in ANSYS Fluent model and 2D EES in- 

house code were compared showing an excellent matching. 
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The only rotor simulations were performed both with laminar and SST com- 
putational models. The employment of Langtry—Menter turbulence model al- 
lowed to determine the presence of transitional effects in the entrance region, 
which were not identified with the laminar model. 

High rotor efficiencies were predicted through the rotor only simulations for 
the small-scale prototype, as high as 69% at 3000 rpm for all assessed fluids. 
Stator—rotor interaction simulations allowed the understanding of the flow 
behaviour due to partial admission. Especially temperature is influenced by 
the streamlines, developing different independent bands of temperature. 
Globally the rotor only and stator—rotor simulations do not affect the predic- 
tion of work and efficiency, nonetheless, the latter simulations allow to in- 
vestigate the critical points of the expander, such as the nozzle outlet opera- 
tion and the influence of the gap. 

The results of the stator-rotor simulations in ANSYS Fluent model and 2D 
EES in-house code were compared showing a good matching of the perfor- 
mance prediction. 


Finally, in order to increase the Technology Readiness Level (TRL) of the re- 
search, experimental activities were necessary. The TRLs are a type of estimation 
system defined in order to evaluate the maturity level of a specific technology, firstly 
developed by NASA in 1974. There are nine technology readiness levels. TRL 1 is 
the lowest and TRL 9 is the highest (NASA, 2018). At the beginning of this research 
the TRL of air Tesla turbine was around 3, while ORC Tesla turbine TRL was between 
1 and 2. At the end of this research work, it can estimated that the TRL for ORC Tesla 
turbine is increased at TRL at least 3, or between 3 and 4. Indeed, as clearly displayed 
in Fig. 5.1, a TRL 3 is considered when active research and technology design are 
performed and when there is the combination between analytical and laboratory stud- 
ies, which means that first preliminary experimental results are required in order to 
for the technology to be considered at “proof of concept” level. 
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Experimental investigations on both air Tesla turbine and ORC Tesla turbine were 
carried out. The experimental assessment of the air Tesla turbine allowed confirming 
the validity of the 2D in-house code, while the experimental campaign on the ORC 
Tesla turbine was the first documented with organic fluids, and therefore allowed a 
proof of concept of the technology. A pivotal point of this last part of research was 
the comparison between experimental and numerical results, which resulted in fair 


agreement. 


The key achievements of this part of research are summarised in the following: 


An experimental test campaign was conducted with air Tesla turbine, with 
two different air mass flow rates at various rotational speeds. The maxi- 
mum achieved efficiency was of 11.2% at 3000 rpm with a total mass flow 
rate of 0.028 kg/s; the maximum power output (again at 3000 rpm, but 
with 0.030 kg/s mass flow rate) was 94 W; 

The experimental results of the air Tesla campaign revealed a very good 
agreement between numerical and experimental results, thus demonstrat- 
ing the soundness of the developed numerical model of the expander; 
The shaft efficiency of the air Tesla prototype (including generator and 
torque meter losses) was estimated over 85%; 

The maximum achievable thermodynamic efficiency of the air Tesla pro- 
totype was about 58%, with power output of 500 W. 

An experimental test campaign was conducted with ORC Tesla turbine, 
with several mass flow rates of R1233zd(E) at various rotational speeds 
and total inlet conditions. The maximum achieved thermodynamic effi- 
ciency was of 30% with a mass flow rate of 0.253 kg/s, a super heating 
level of 37.5 °C and a rotational speed of 3500 rpm; while the maximum 
experimental thermodynamic power obtained was of 906 W, with a mass 
flow rate of 0.299 kg/s, a super heating level of 44 °C and a rotational 
speed of 5000 rpm; 

The maximum achieved shaft efficiency of the ORC Tesla turbine was of 
9.62 % with a mass flow rate of 0.299 kg/s, a super heating level of 45.5 
°C and a rotational speed of 4000 rpm; while the maximum experimental 
shaft power obtained was of 371 W for two thermodynamic states: for a 
mass flow rate of 0.299 kg/s, a super heating level of 45.5 °C and a rota- 
tional speed of 4000 rpm; and for a mass flow rate of 0.365 kg/s, a super 
heating level of 4.65 °C and a rotational speed of 3500 rpm; 

The shaft efficiency of the ORC Tesla prototype (including generator and 
torque meter losses) was estimated at 50%, which is a really low value, 
but which was mostly due to the not right alignment of the magnetic cou- 
pling, which brought a really high increase of mechanical losses, due to 
contact friction losses; 

The modified 2D EES in-house code, taking into account the partial ad- 
mission, windage and pumping losses, as well as the experimental constant 
coefficient properly described the trend of both thermodynamic power and 
efficiency, as the numerical prediction was always within the uncertainty 
of the measured data. 
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As a final remark, it can be affirmed that the here proposed analysis showed how 
the Tesla turbine might be a competitive solution when applied to micro and small 
power applications, from recovery of low pressure waste air flows to micro ORC ap- 
plications. The research on this expander technology is worth of further develop- 
ments, given the encouraging experimental results which demonstrated the reliability 
of the numerical simulations and that predicted very interesting power output and ef- 
ficiency potential for a very wide range of possible micro power applications in the 
kW scale. 


2. Recommendations 


The most difficult task of a project is to decide its boundaries and to properly 
select and assess only the fundamental parts, while leaving some interesting aspects 
to be developed in proper future works. Hereafter a list of possible lines of research 
is presented, which could rise from the leveraging of the acquired know-how from 
this research. 


e Perform further experimental tests on air Tesla turbine with different mass 
flow rates, rotational speeds and total inlet conditions, in order to validate the 
model for a wide range of conditions; 

e Re-design the ORC Tesla turbine, leveraging the experimental data obtained, 
improving the stator—rotor coupling and substituting the electromagnetic cou- 
pling; 

e Test the ORC Tesla turbine in a wider range of conditions, with higher ex- 
pansion ratios; 

e Test the ORC Tesla turbine with different working fluids, in order to further 
prove the suitability of this technology for ORC applications; 

e Design a Tesla turbine for СО» applications, as СО» expansion ratios could 
properly match the optimal range of Tesla turbine expansion ratios; 

e Develop a two-phase flow model for the Tesla turbine, combining two-phase 
flow and Tesla turbine theories; 

e Compute numerical analyses on two—phase flow in a Tesla turbine rotor; 

e Design a Tesla turbine for two—phase flows; 

e Perform an experimental campaign on two-phase Tesla turbines. 
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ORC Organic Rankine Cycle 

P Power [W] 

р Pressure [Ра] 

Pr Prandtl number [-] 

Q Heat [W] 

r Radius [m] 

R Rotor Inlet/Outlet diameter ratio 
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Re Reynolds number [—] 

rpm Revolution per minute 
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Se Length of viscous entry region 
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SH Super heating level [?C] 
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В Relative flow angle [°] 
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Subscripts and superscripts 
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Micro turbo expander design for small scale ORC 


The Tesla expander was first developed by N. Tesla at the beginning of the 20" 
century. In recent years, due to the increasing appeal towards micro power genera- 
tion and energy recovery from wasted flows, this cost effective expander technol- 
ogy rose a renovated interest. In the present study, a 2D numerical model is real- 
ized and a design procedure of a Tesla turbine for ORC applications is proposed. 
A throughout optimization method is developed by evaluating the losses of each 
component. The 2D model results are further exploited through the development 
of 3D computational investigation, which allows an accurate comprehension of the 
flow characteristics. Finally, two prototypes are designed, realized and tested. The 
former one is designed to work with air as working fluid. The second prototype is 
designed to work with organic fluids. The achieved experimental results confirmed 
the validity and the large potential applicative chances of this emerging technology 
in the field of micro sizes, low inlet temperature and low expansion ratios. 
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